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The  increased  complexity  of  the  very  large  scale  integrated  cir¬ 
cuits  (VLSI)  has  greatly  impacted  the  field  of  computer-aided  design 
(CAD).  One  of  the  problems  brought  about  is  the  interconnect ion 
problem. 

In  this  research,  the  goal  is  twofold.  Firat  of  all,  a  more 
accurate  numerical  method  to  evaluate  the  interconnect  capacitance, 
including  the  coupling  capacitance  between  interconnects  and  the 
fringing  field  capacitance,  waa  investigated,  and  the  integral  method 
was  employed.  Two  FORTRAN  programs  "CAP2D"  and  "CAP3D*  based  on  this 
method  were  developed.  Second,  a  PASCAL  extraction  program 
eaiphasizing  the  extraction  of  interconnect  parameters  was  developed. 
It  employs  the  "cylindrical  approximation  formula"  for  the  self¬ 
capacitance  of  a  single  interconnect  and  other  simple  formulas  for 
the  coupling  capacitances  derived  by  a  "least  square  method".  The 
extractor  assumes  only  Manhattan  geometry  and  NMOS  technology. 
Four— disMnaional  binary  search  trees  are  used  as  the  basic  data 
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INTRODUCTION 


The  advent  of  very-large-scale  integrated  (VLSI)  circuits  as  a 
resnlt  of  rapid  progress  lias  greatly  iapaeted  on  the  various  aspects 
in  aieroelectronics.  For  exaaple,  as  feature  diaensions  approach  the 
snbaicron  range,  the  "short-channel  effects"  doainate  the  physics  of 
the  devioes.  This  has  brought  forth  aore  research  and  studies  into 
device  physics,  processing  technologies,  aaterial  science,  etc.  On 
the  other  hand,  the  exponential  growth  of  IC  coaplexity  aay  push  to 
the  extent  that  the  future  progress  in  circuit  integration  vill  no 
longer  be  liaited  only  by  technology,  but  also,  and  to  a  auch  greater 
extent,  by  our  ability  to  reduce  the  astronoaical  design  tiae.  This 
concern  has  spawned  research  and  studies  in  coaputer-aided-design 
(CAD),  coaputer  architecture,  design  aethodology,  design  autoaation, 
etc.  (1-491.  All  these  efforts  are  aiaed  at  further  advancing  the 
progress  in  aieroelectronics  that  we  have  enjoyed  for  the  past  two 
decades.  Hopefully,  through  these  studies  and  with  the  coaputational 
power  VLSI  circuits  have  brought  us,  we  can  "bootstrap"  ourselves  to 
the  next  generation  of  technologies. 

A  eoaprehensive  CAD  systea  consisting  of  various  tools  has 
beeoae  indispensable  in  VLSI  circuit  design.  A  part  of  the  CAD  sys- 
tea  is  design  verification  whioh  shows  that  a  given  mask  geoaetry 
vill  perfora  the  desired  digital  functions  as  originally  designed. 


Extraction,  the  first  step  in  design  verification,  transforms  the 
geometric  information  on  the  masks  into  circuit  information.  After¬ 
wards  a  comparison  can  be  made  with  the  original  design,  or.  a  circuit 
simulation  can  be  performed  to  check  the  responses  of  the  circuit. 

One  of  the  problems  that  the  increased  complexity  has  brought 
about  in  CAD  is  the  interconnection  problem.  Not  only  does  the  rout¬ 
ing  among  s  myriad  of  devices  on  the  chip  become  coa^licated,  diffi¬ 
cult  and  time-consuming,  but  also  the  electrical  parameters  of  these 
very  interconnects  have  gradually  grown  into  an  integral  part  of  the 
circuits.  The  increased  complexity  has  augmented  the  lengths  of 
interconnects  on  the  chip.  Moreover,  the  reduced  dimensions  of  the 
interconnects  themselves  increase  the  resistances  of  the  lines,  snd 
the  fringing  capacitances  are  no  longer  negligible.  Consequently,  it 
is  necessary  to  include  these  interconnect  parameters  in  the  circuit 
model.  It  is  on  this  aspect  that  this  thesis  concentrates. 

First  of  all,  an  equivalent  circuit  model  for  the  interconnects 
on  the  chip  has  to  be  established.  A  T-type  lumped  circuit  model  is 
used  in  this  thesis  for  the  interconnects  assuming  silicon  technol¬ 
ogy.  In  Appendix  A,  oriteria  for  the  validity  of  the  lumped  circuit 
model  are  derived  from  transmission  line  theory.  It  is  shown  that  the 
lumped  circuit  models  are  generally  applicable  under  the  present 
technologies. 

In  the  first  part  of  the  thesis,  s  more  accurate  account  of  the 
interconnect  capacitance  by  numerical  methods  is  considered.  Both 
the  differential  method,  e.g.,  the  finite  difference  method,  snd  the 
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integral  method  (by  way  of  Green's  functions)  have  been  studied. 
Based  on  the  applioations  of  both  methods  to  solve  various  problems 
in  the  literature  [2-10],  combined  with  the  experience  in  this  study, 
it  was  found  that  the  integral  method  seems  to  offer  more  advantages 
in  the  flexibility  of  programming  and  in  computation  time.  There¬ 
fore,  it  vas  employed  for  the  numerical  computation  of  two-  and 
three-dimensional  interconnect  capacitances.  Basically  the  method 
follows  the  "method  of  moments"  proposed  by  Harrington  [11]  and  is 
similar  to  those  in  [8-10] .  Both  testing  functions  and  basis  func¬ 
tions  are  ohosen  the  same  aa  pulse-type  functions,  which  are  con¬ 
stants  over  the  subdivisions,  and  variable-length  subsections  are 
divided  according  to  the  roots  of  Chebyshev  polynomials.  It  vas 
found  that  this  approach  offered  more  accuracy  than  a  comparable 
number  of  constant  subdivisions  or  using  delta-type  testing  func¬ 
tions,  and  the  positive-definiteness  of  the  resulting  kernel  matrix 
also  guaranteed  a  solution.  Some  results  will  be  presented  in 
Chapter  2. 

Next,  with  a  view  to  easily  incorporating  these  capacitances  in 
the  extraction  program  for  design  verification,  simple  approximation 
formulas  were  sought.  An  approximation  formula  for  the  self¬ 
capacitance  of  a  single  interconnect  was  first  derived.  It  is  based 
on  the  known  expression  of  the  capacitance  of  a  cylinder  over  a 
ground  plane.  Employing  this  expression  to  approximate  the  fringing 
part  of  the  self-capacitance,  we  have  a  reasonable  approximation  for¬ 
mula  with  an  accuracy  within  10%  of  the  numerically  computed  values. 


and  it  is  named  "cylindrical  approximation  formula".  As  to  the  cou¬ 
pling  capacitances  between  interconnects,  a  "least  square  fit"  method 
on  the  data  obtained  from  the  numerical  computation  was  employed  to 
find  approximation  formulas.  A  few  formulas  have  been  obtained  for 
different  interconnect  configurations,  e.g.,  coupling  capacitances 
between  two  parallel  interconnects  or  two  pe  jndicular  cross-over 
lines . 

The  second  part  of  the  theais  is  mainly  unearned  with  the 
extraction  program.  It  was  intended,  if  possible,  to  be  part  of  a 
hierarchical  CAD  system  which  supports  a  number  of  functions  at  each 
level  in  the  design  hierarchy  with  different  tools.  Hence,  the  out¬ 
put  form  of  the  extractor  is  in  a  SPICE  input  compatible  format  such 
that  the  circuit  simulation  can  be  performed  without  the  data  conver¬ 
sion.  The  basic  data  structure  in  the  extractor  is  4-d  binary  search 
trees  of  geometrical  rectangles  on  the  masks.  The  output  includes 
not  only  the  pertinent  information  of  transistors,  e.g.,  channel 
width  and  channel  length,  but  also  the  parameters  of  interconnects. 
The  resistances  and  the  capacitances,  both  self-  and  coupling  capaci¬ 
tances,  are  reported. 

In  aummary.  Chapter  2  focuses  on  the  numerical  computation  of 
the  interconnect  capacitance,  both  in  two  and  three  dimensions.  It 
comprises  the  detailed  formulation  of  the  Green's  functions,  the 
even-odd  analysis  to  find  the  coupling  capacitances  in  multi¬ 
conductor  cases,  the  integral  method  and  the  comparison  with  the  dif¬ 
ferential  method,  and  a  few  examples.  A  user's  manual  for  the  pro- 


grass  "CAP2D",  "CAP3D"  which  execute  the  aforementioned  numerical 
computation  is  described  in  Appendix  B.  Chapter  3  covers  the  approx¬ 
imation  formulas  of  capacitances.  The  "cylindrical  approximation 
formula"  is  derived,  and  other  approximation  formulas  for  the  cou¬ 
pling  capacitances  are  also  discussed.  In  Chapter  4,  the  network 
extraction  program's  data  structure  and  the  procedure  of  extracting 
transistors  and  interconnect  parameters  are  detailed;  its  charac¬ 
teristics,  limitations  and  some  examples  will  also  be  discussed. 
Finally,  the  conclusions  are  in  Chapter  5. 


NUMERICAL  CONFUTATION  OF  INTERCONNECT  CAPACITANCE 


i.l.  Introduction 

The  present  very  large  scale  integration  circuits  lead  to  very 
large  and  complex  systems  with,  small  physical  dimensions.  Not  only 
do  the  devices  have  small  dimensions  which  require  special  attention, 
but  the  interconnects  also  become  important  and  have  to  be  reckoned 
with  in  the  circuit  simulation.  A  better  understanding  of  the 
electrical  characteristics  of  the, interconnect,  especially  the  capa¬ 
citance.  is  essential  in  the  design  and  analysis  of  today's 
integrated  circuits. 

As  a  general  rule  of  thumb,  the  highest  frequency  components 
contained  in  the  signals  in  most  present  large  scale  integrated  cir¬ 
cuits  correspond  to  wavelengths  well  exceeding  the  physical  dimen¬ 
sions  of  the  circuit  elements.  Hence,  leaped  circuit  models  can  be 
used  for  the  electrical  analysis  of  the  circuits.  6.  Bilardi  et 
al.  [12]  evaluated  various  circuit  models  for  interconects  in  VLSI 
circuits  and  it  was  concluded  that  both  current  and  the  projected 
silicon  technologies  in  the  near  future  fall  within  the  realm  of  the 
capacitance  model.  That  is  to  say  that  a  dispersive  line  esn  be 
replaced  by  a  capacitance  proportional  to  its  length.  From  another 
point  of  view,  the  lumped  circuit  model  can  be  considered  as  a 


first-order  approximation  of  the  transmission  line  model  as  discussed 
in  Appendix  A. 

Another  important  aspect  in  the  eleotrical  analysis  of  the 
integrated  oircuit  is  the  "general  impedance  level".  It  is  defined 
in  [13]  as  the  lossless  characteristic  impedance  of  the  average  con¬ 
nection  in  the  system.  Since  most  digital  integrated  circuits  are 
based  on  JtOSFET  devices  which  are  typically  of  higher  impedance  than 
the  general  impedance  level,  voltage  and  charge  become  the  "state" 
variables  in  the  analysis.  Consequently,  capacitance  emerges  as  the 
dominant  circuit  element.  The  deciding  factor  of  the  performance  of 
most  digital  circuits  is  the  time  delay  through  the  system,  and 
interconnect  capacitance  certainly  plays  an  important  role. 

A  cross-sectional  view  of  a  typical  integrated  circuit  is  shown 
in  Fig.  2.1.  Since  the  channel-stopping  ion  implantation  is  usually 
situated  beneath  the  field  oxide  where  the  interconnect  is  generally 
placed,  it  can  be  considered  as  the  ground  plane  as  shown  in  the  com¬ 
putation  model  in  Pig.  2.2.  In  general  cases,  there  can  be  three 
different  dielectric  media,  and  interconnects  can  reside  inside  the 
oxide  layer  or  the  passivation  layer  as  shown  in  Fig.  2.2.  The  top¬ 
most  layer  is  assumed  to  be  always  the  free  space.  In  order  to 
reduce  the  computational  effort,  we  can  also  consider  the  intercon¬ 
nect  capacitance  in  a  two  layer  or  a  homogeneous  medium  if  the  layers 
in  question  are  thick  enough.  The  following  discussion  of  the  numeri¬ 
cal  computational  method  of  the  interconnect  capacitance  is  based  on 
this  model  as  shown  in  Fig.  2.2. 


£.£.  Numerical  Coi 


The  computation  of  capacitance*  has  appeared  in  the  past  for 
various  applications,  e.g.,  the  calculation  of  the  characteristic 
impedance  of  microstrip  transmission  lines  is  one  prime  example  [42- 
45].  Thongh  the  emphasis  and  the  assumptions  are  different  in  dif¬ 
ferent  applications,  the  methods  are  still  applicable  in  the  computa¬ 
tion  of  interconnect  capacitance.  In  summary,  the  available  computa¬ 
tion  methods  of  interconnect  capacitance  can  be  categorized  into 
three  types,  namely:  analytic,  differential  and  integral  methods. 

The  analytic  method,  in  effect,  can  calculate  the  capacitance  in 
closed  form  and  no  numerical  analysis  or  discretization  is  required. 
This  will  be  the  subject  of  the  next  chapter  so  its  discussion  vill 
be  deferred.  In  this  ohapter  we  vill  concentrate  on  the  numerical 
methods,  i.e.,  the  differential  and  integral  methods. 

Let  na  first  consider  the  differential  method.  The  finite 
difference  method  for  solving  boundary  value  problems  is  the  prime 
example  of  a  differential  method.  It  starts  vith  Poisson's  differen¬ 
tial  equation,  defines  the  boundary  of  the  problem  in  question,  and 
then  sets  up  the  boundary  conditions.  After  solving  all  the  poten¬ 
tials  on  the  grid  points  vithin  the  boundaries,  ve  can  calculate  the 
charges  on  the  surface  of  the  interconnect  and  then  the  capacitances 
accordingly.  For  example,  the  method  employed  in  [2,3]  can  be 
categorized  as  the  differential  method.  It  is  generally  easy  and 
straightforward  to  formulate.  Multiple  layers  of  media  or  irregular 


boundaries  can  be  accommodated  without  too  much  inc; 


in 


computation.  However,  if  the  region  in  consideration  is  an  open 
boundary  aa  shown  in  Fig.  2.2,  a  suitable  boundary  condition  has  to 
be  defined.  In  [2]  a  "ref lective"-  or  Neumann-type  boundary  condi¬ 
tion  at  a  reasonable  distance  from  the  source  was  defined  to  simulate 
the  open  boundary.  In  [3]  an  ad  hoc  arrangement  on  the  boundary 
grids  was  set  up  in  the  numerical  process  to  take  care  of  the  open 
boundary.  Though  there  are  other  ways  to  reasonably  simulate  the 
open  boundary,  they  may  compromise  the  accuracy  of  the  computation  or 
increase  the  computation  time.  This  is  one  disadvantage  of  the  dif¬ 
ferential  method.  Other  disadvantages  may  be  attributed  to  the 
rather  large  system  of  linear  algebraic  equations  of  grid  potentials 
and  the  difficulty  of  extending  it  to  three-dimensional  cases. 
Depending  on  the  grid  and  the  regularity  of  the  boundaries  and  the 
source  regions,  this  matrix  may  not  lead  itself  to  some  special  tech¬ 
niques  in  numerical  linear  algebra.  Iterative  methods  like  Gauss- 
Seidel  may  be  too  time-consuming  and  slowly  convergent.  In  three- 
dimensional  cases,  the  number  of  discretixed  grids  may  be  too  large 
to  handle  in  the  finite  difference  method. 

On  the  other  hand,  the  integral  method  starts  with  the  Green's 
function  of  the  Poisson’ s  equation,  and  the  integral  equation  of  the 
convolution  of  the  Green’s  funetion  with  the  charge  density  function 
is  solved.  Then  the  capacitances  are  calculated  after  the  charge 
denaities  on  the  interconnect  have  been  found.  In  [4]  the  Green's 
function  for  two-dimensional  interconnects  in  two  layers  of  media  was 
derived,  and  it  will  be  shown  in  a  latter  section.  In  [6]  the 


three— dimensional  Green's  function  wee  formulated  for  the  computation 
of  irregular  conductors  over  dielectrics. 

In  comparing  these  two  methods,  the  differential  method  can  be 
thought  of  as  solving  the  problem  "microscopically",  and  the  integral 
aMthod  as  "maeroscopically".  Because  the  system  of  linear  algebraic 
equations  generated  in  the  finite  difference  scheme  relates  only 
local  neighboring  points  in  the  grid,  the  resulting  matrix  tends  to 
be  large  but  sparse.  The  effects  of  the  boundaries  and  the  sources 
on  the  grid  potentials  are  reflected  only  after  the  system  of  equa¬ 
tions  is  solved  by,  for  exaaple,  the  Ganss-Seidel  method.  Vhile,  on 
the  other  hand,  the  effects  of  the  boundaries  are  already 
"integrated"  in  the  Green's  function  in  the  integral  swthod;  thus, 
the  region  in  question  only  concerns  the  source  region,  and  the 
resulting  kernel  matrix  is  small  but  full  in  this  case.  There  are 
other  differences  between  these  two  methods,  and  they  are  summarized 
in  Table  2.1. 

• 

One  of  the  important  characteristics  of  these  two  methods  is 
that  the  capacitance  computed  by  the  integral  method  will  be  a  lower 
bound  to  the  true  value  and  that  by  the  differential  method  will  be 
an  upper  bound.  This  is  parallel  to  the  computation  of  the  charac¬ 
teristic  impedance  of  a  waveguide  discussed  in  [41],  where  the  varia¬ 
tional  properties  of  the  two  approaches  were  also  discussed.  As  in 
[41],  it  would  be  ideal  to  compute  the  capacitance  twice,  once  by 
each  method  and  take  the  average.  However,  there  may  be  too  much 
computation  involved  this  way.  Besides,  the  accuracy  of  the  results 


TABLE  2.1  Comparison  of  the  Differential  and  Integral  Method*. 


Itea  of  Coaparison 

Differential  Method 

Integral  Method 

Basie  equation 

Poisson's  differen¬ 
tial  equation 

Green's  function  and 
its  convolution  in¬ 
tegral 

"State"  variable* 

Potentials 

Charge  densities 

Diseretised  region 

Large  region  enclos¬ 
ing  the  source  and 
the  aedia  around  it 

Only  the  surface  of 
the  source  region 

Multiple  layers  of 
aedia 

Easy  to  acccaaodate 
without  too  auch  in¬ 
crease  in  coaputation 

Infinite  snaaation  in 
Green's  function 
causes  a  substantial 

tion 

Multiple  conductors 

Cannot  take  advantage 
of  the  even-odd  node 
excitation,  substan¬ 
tial  increase  in  coa¬ 
putation;  hard  to  su- 
toaate  in  progranaing 

Easy  to  _  take  advan¬ 
tage  of  the  even-odd 

node  excitation,  and 
easy  to  automate  the 
input  of  the  diaes- 
sions  and  positions 
of  the  conductors 

Autoaatic  aesh  gen¬ 
eration 

Harder  to  ia^pleaeat 
in  a  prograa 

Easy  to  iaqpleaent 

Three-diaens ional 

ease 

Much  aore  complicated 
than  two-dimensional 
counterpart 

Not  too  auch  aore 
coaplicated 

Solution  aethod 

Solution  of  a  large 
linear  algebraic 
equation,  i.e., 
Oauss-Seidel 

Moaent  aethod 

Boundaries  shape 

Can  take  care  of  aore 
irregular  boundaries 
provided  the  grid  is 
refined  enough 

Rigid,  regular  boun¬ 
daries,  at  aost  slant 
side  walls  and  plane 
interfaces 

Open  boundaries 

Soae  approxiaation  or 
other  aechaaisa  has 
to  be  devised 

No  problea  with  open 
boundaries 

Characteristics  of 

Upper  bound  on  the 
true  vain* 

Lower  bound  on  the 
true  value 

*  State  variable  seen a  the  unknown  in  the  ba*ie  equation  first  was 
solved  bp  either  the  differential  or  integral  aethod. 

**  See  Ref.  [41]  by  Collin  for  a  complete  discnssioa. 


can  be  increased  by  refining  the  grid  in  either  method.  In  comparing 
the  pros  and  cons  between  these  two  methods,  the  integral  method 
seems  to  be  the  better  choice,  especially  in  the  case  of  homogeneous 
medium  and  multiple  conductors.  It  is  easy  to  program  for  multiple 
conductors  with  diverse  geometries  and  positions;  the  discretized 
region  only  concerns  the  conductor  surfaces,  which  reduces  the  number 
of  linear  algebraic  equations.  More  important,  we  can  take  advantage 
of  the  LU  factorization  of  the  kernel  matrix  in  solving  charge  densi¬ 
ties  for  the  multiple  conductor  case.  However,  two  major  disadvan¬ 
tages  of  the  integral  method  are  in  the  cases  of  multiple  layers  of 
media  and  irregular  shapes  of  the  interconnects  and  the  boundaries. 
If  accurate  numerical  computation  is  imperative  for  those  cases,  the 
finite  element  method,  one  of  the  differential  methods,  may  be  the 
only  viable  solution. 


2.£.  Formulation  of  Multi-Conductor  Capacitance  Computation 

In  this  section,  a  general  formulation  for  computing  the  self- 
and  the  coupling  capacitances  among  multiple  conductors  is  described. 
Either  the  differential  or  the  integral  method  has  to  be  associated 
in  this  frame  of  formulation  to  compute  the  capacitances.  First  the 


general  case  of  N  conductors  is  derived,  then  two  examples  of  N  -  3 
and  N  ■  4  are  shown.  Special  consideration  is  required  for  the  case 


V 


1.2.1.  General  formal at ion  for  N  conductors 


Conaidc 


lystem  of  N  conductors  over  a  ground  plane  as  shown 


in  Fig.  2.3  where  N  ■  3  is  depicted.  All  the  capacitances,  both 
self-  and  coupling,  are  shown.  The  algebraic  relationship  between 
the  charges  Qi#  i«l,2,*..,N  and  the  potentials  0i#  i-l,2,...,N  on  the 
conductors  can  be  expressed  as  [14] 


. + 


0*  ■  C*i,(0a“^x)  +  + 


+  C2n(®*“®n> 


(2.1) 


Qn  "  Cul<V«*>  +  *  •••  +  CM0 


n2  n 


Note  that  the  potentials  0's  are  referenced  with  respect  to  the 
ground  plane,  and  the  capacitances  are  symmetric,  i.e.,  -  Cji# 
for  1  <  i,j  £  N,  because  of  reciprocity. 

An  analysis  naaed  "even-odd  node  analysis"  is  eaqployed  to  solve 


the  capacitances  C^'s  in  Eq.  (2.1).  Different  potential  patterns 
are  first  assigned  on  the  conductors,  even  node  or  odd  modes,  then 
the  corresponding  charges  can  be  expressed  algebraically  in  terms  of 
the  capacitances  as  will  be  shown  in  the  N  -  3  and  N  -  4  examples. 
Solving  for  the  charges  numerically  in  the  respective  modes  by  either 
the  differential  or  the  integral  method,  then  we  can  find  easily  all 
the  capacitances. 
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3 


s’. 


fj 


( i)  Even  mode  :  la  the  even  mode  the  potentials  are  assigned  as 
•a  *  0,  0^  ■  1.  From  Eq.  (2.1)  the  self-capacitances  can  be 
written  as  ,  i-l,2,...,N.  The  superscript  on  the  charges 
indicates  that  those  charges  are  evaluated  under  the  even-node  poten¬ 
tial  pattern.  All  the  self-capacitances  are  determined  by  this 
even^ande  excitation. 

Since  there  are  N(N-l)/2  coupling  capacitances  left  to  be 
determined,  it  is  necessary  to  generate  the  sane  number  of  indepen¬ 
dent  algebraic  equations  from  the  odd-mode  potential  patterns.  It  is 
a  matter  of  choice  to  assign  these  odd-mode  potential  patterns.  A 
systematic  procedure  is  used  here.  Except  for  the  case  N  *  4,  the 
procedure  employed  here  produces  a  system  of  N(N-l)/2  linear  indepen¬ 
dent  algebraic  equations  which  can  be  programmed  easily  and  then 
solved  uniquely  for  the  coupling  capacitances. 

( ii)  Odd  mode  1  :  the  potentials  are  assigned  as  01  *  0,  -  i, 
and  all  the  other  potentials  equal  -1.  Solving  Eq.  (2.1),  we  have 

®a  “  ^ii  +  2(C1J+C14+» • »+C^a) 

“  Caa  +  2(C„+C>4+...+C2ft) 

The  sum  of  the  above  two  equations  is  chosen  to  be  one  of  the  N( N- 
l)/2  independent  equations.  Since  Cxx  and  Caa  have  been  found  in  the 
even-mode  analysis,  the  sum  of  the  two  equations  can  be  expressed  as 


TV-.- 


^ 
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cx»  +  cX4  +  ...  +  cla  +  ca,  +  c14  +  ...  +  c2a 

-  0.5(Q*X)  +  <4X>  -  q£#)  -  ai#)>  (2.2) 

The  superscript  1  on  the  charges  here  indicates  the  first  odd-mode 
ezoitation. 

Other  odd-aK>de  patterns  can  be  obtained  as  in  (ii)  by  assigning 
two  of  the  N  potentials  1  volt,  and  all  the  other  potentials  -1  volt. 
Since  there  are  totally  C(N,2)  *  N(N-l)/2  comb  ins  t  ions  of  this  type, 
a  total  of  N(N-l)/2  linear  algebraic  equations  as  Eq.  (2.2)  can  be 
generated.  It  can  be  shown  that  the  determinants  of  these  systems  of 
algebraic  equations  are  nonzero  for  N  *  3  to  10  except  N  ■  4.  So  the 
linear  independence  is  achieved  and  a  unique  solution  is  guaranteed. 

Example  for  H  "  i 

First  Eq.  (2.1)  is  rewritten  in  the  following  form  for  N  *  3. 

Qi  *  CXX9X  +  CXX(0X-9X)  +  Ci,(0x-0,)  (2.3a) 

®s  «  CXX(9X-0X)  +  Cxx9»  +  C,,(0,-0,)  (2.3b) 

*  CXj(0|—0i)  +  Caj(0j— 0j)  +  C||0|  (2.3c) 

The  evenr-  and  odd-mode  excitations  are  listed  as  follows  : 
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( i)  Even  mode  :  let  a  0,  3  0,  *  1  then  we  have 

C  m  •  C  —  a(#)  *  C  m 

'*11  Qj  '  '"ll 

( il)  Odd  node  1  :  let  0X  ■  0*  “  1*  and  0S  ■  -1,  then  ve  have 

cn  +  Clt  -  0.5( q£,1)  +  Q^1)  “  ®ie)  "  0i#)) 

( iii)  Odd  node  2  :  let  0X  »  0,  =  1,  and  0X  *»  -1,  then  we  have 

Cxa  +  C„  -  0.5(q12)  +  oj**  -  Qie)  ~  Qie)) 

( iv)  Odd  node  3  :  let  0X  »  0,  *  1,  and  0X  ■  -1,  then  we  have 

Cxx  +  Cxl  -  0.3(Qi,)  +  Q<’>  ~  Q$e)  "  Q(o)) 

It  can  be  easily  shown  that  the  determinant  of  the  system  of  the 
above  three  linear  algebraic  equations  from  (ii),(iii)  and  (iv)  is  2. 
Hence,  the  coupling  capacitances  Cxs,  Cx,  and  CX1  are  solved  in  terms 
of  the  corresponding  charges  as 


(2.4.) 
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C,,  -  0.25(Ql[1)-Qi,)+aS1)-Qi,>+Qi1)+Qi1)-2Qie>) 


(2.4c) 


Note  that  the  subscript  i  on  the  charge  indicates  the  nnaber  of 
the  conductor  the  charge  is  on,  while  the  superscript  j  shows  the 
nnaber  of  the  odd-aode  excitation  this  charge  is  extracted  from. 
Since  there  are  other  ways  to  define  the  odd-mode  potential  patterns, 
the  coupling  capacitance  formula  in.  Eq.  (2.4)  is  not  unique.  For 
exasiple,  the  capacitance  formula  in  [5]  is  'the  result  of  another  way 
to  specify  the  odd-mode  excitation  and  not  the  same  as  that  in 
Eq.  (2.4). 


2.2..i.  Special  consideration  for  N  »  4 

The  procedure  described  in  Sec.  2.3.1  would  result  in  a  system 
of  linearly  dependent  equations  for  the  case  N  *  4.  Another  way  to 
specify  the  even-odd  mode  excitation  has  to  be  devised  here.  The 
even-mode  excitation  can  still  be  the  same,  while  some  other  odd-mode 
excitations  have  to  be  defined.  Instead  of  snouting  the  corresponding 
equations  resulting  from  the  odd-mode  excitations  as  mentioned  above, 
they  can  be  treated  as  separate  equations.  The  resulting  system  of 
linear  algebraic  equations  can  be  shown  to  be  independent.  So  the 
even-odd  mode  analysis  for  the  N  -  4  case  can  be  summarized  as  fol¬ 
lows  : 


even  mode  :  0X  >  0a  ■  0,  -  04 


1 


Equations  (2.5)  to  (2.10)  constitute  s  system  of  six  algebraic  equa¬ 
tions  with  determinant  4.  Thus,  the  six  coupling  capacitances  are 


uniquely  determined.  Note  that  if  the  general  procedure  in 
Sec.  2.3.1  ia  employed,  then  the  next  odd-mode  pattern  will  be 
t%  ■  0|  *  1.  0X  ■  04  ■  -1.  This  is  the  exact  complement  to  odd-mode 
3  listed  above.  Consequently,  the  summing  of  the  corresponding  equa¬ 
tions  vill  result  in  an  identical  left-hand  side  from  both  excita¬ 
tions  and  a  dependent  system  of  equations.  For  cases  other  than 
N  ■  4,  it  is  not  possible  to  find  pairs  of  "complementary  excita¬ 
tions";  hence,  the  equations  are  independent. 


£.±.  Integral  Method 

In  this  section,  a  detailed  discussion  of  the  integral  method  is 
given.  The  method  of  moments  [11]  is  formulated,  *2ong  with  the 
incorporation  of  the  "even-odd  mode  analysis"  discussed  in  Sec.  2.3. 
The  two-  and  three-  dimensional  Green's  functions  in  multiple  layers 
of  media  vill  be  derived  in  the  next  two  sections,  and  the  detailed 
expressions  used  to  compute  the  kernel  matrix  vill  be  given  in 
Appendices  C  and  0. 


Assume  that  the  Green's  function  to  the  Poisson's  equation  is 
expressed  as  G(r,r').  In  the  derivation  of  this  Green's  function, 
the  boundary  conditions  have  already  been  taken  into  account.  In 
analogy  to  system  analysis,  the  Green's  function  can  be  considered  as 
the  impulse  response,  and  the  potential  in  Poisson's  equation  can 
then  be  expressed  as  a  convolution  integral  : 


0(r)  ■  Js  dr# 


(2.11) 


la  Eq.  (2.11)  0  is  the  poteatisl  sad  o  is  the  charge  dsasity  oa  the 
scares  sarfacs  whose  area  is  designated  as  S.  Note  that  ths  duaay 
variable  ia  ths  iategral  r'  indicates  ths  coordiaats  of  ths  soarcs 
region,  aad  r  is  ths  observation  poiat.  Since  potentials  oa  ths 
source  rsgioa  are  decided  by  ths  evenrodd  aods  excitation  as  dis¬ 
cussed  in  ths  previous  section,  aad  ths  Green's  function  is  also 
known,  the  only  unknown  variable  ia  Eq.  (2.11)  is  the  charge  density 
function  o(r').  Therefore,  it  is  an  iategral  equation,  aad  the 
aethod  of  aoaeats  can  be  eaployed  to  solve  it. 

First,  assaae  that  the  unkanwu  charge  density  function  o(r') 
can  be  expressed  in  terns  of  a  set  of  basis  functions  w  i-l,2,***,N 


a 

«(r')  •  J  °i«i(*') 


(2.12) 


Thus,  the  unknown  function  o(r»)  ia  trsnforaed  into  N  unknown  con¬ 
stant  ooef feeieats  c^.  Substituting  Eq.  (2.12)  into  Eq.  (2.11),  we 


an 

•C^ls-5  Jg  wl(r')G(^,7') 


(2.13) 


Note  that  the  observation  point  is  chosen  at  the  source  surface  S 
over  which  the  convolution  integral  is  integrating,  because  in  even- 
odd  node  excitation,  the  potentials  are  specified  over  S. 
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Next,  another  set  of  testing  functions  t  ,  j«l,2,...,N  is 

defined.  If  the  testing  functions  ere  chosen  exactly  the  same  as  the 
basis  functions,  the  aethod  is  named  Galerkin's  aethod  [11],  Defining 
an  inner  product  <f.g>  ■  Jgf(r)g(r)dr  for  any  two  real  functions  f 
and  g.  we  then  take  the  inner  products  on  both  sides  of  Eq.  (2.13) 
with  all  the  testing  functions  t ^ ,  sad  have 


f  a 

J  0(r)  I  t.UJd?  -  J  «J  «£j  <£'  -iU'JGU.J'HjU) 

S  S  c  s 


s  Js 


j  -  1,2, *  *  * ,N 


(2.14) 


Let  us  further  aiaplify  Eq.(2.14)  by  defining 


bj  ■  J  0(r) I  tj(r)dr 
14  S  S 


and  A 


ij 


i, j  ■  1 ,2, . . .  ,N 


Since  b^  and  Ajj  are  scalars  which  can  be  calculated  froa  the  known 
functions  0(r)|g,  tj(r),  a^r'),  and  G(r.r'),  Eq.  (2.14)  is 
essentially  a  systea  of  siaultsneous  linear  algebraic  equations.  It 
can  be  written  in  the  aatrix  fora  Ax  *  b.  The  vector  x  is  the  un¬ 
known  coefficients  ai  to  be  solved.  Note  that  the  sixe  of  the  "ker¬ 
nel"  aatrix  A  depends  on  the  nuaber  of  basis  functions.  Usually  it 
is  considerably  saaller  than  that  generated  in  the  finite  difference 
scheae,  so  a  direct  aethod  like  LU  decoaposition  or  Gaussian  elimina¬ 
tion  can  be  eaployed.  It  should  also  be  pointed  out  that  tiie 


"kernel"  matrix  is  usnally  full  and  tkose  obtained  by  Galerkin's 


method  should  be  symmetric. 

One  of  the  main  tasks  in  the  above  procedure  is  the  ohoioe  of 
the  basis  functions.  As  mentioned  in  [11] .  they  should  be  linearly 
independent*  and  the  superposition  of  them  as  in  Eq.  (2.12)  should 
approximate  the  unknown  function  b(r)  as  olosely  as  possible.  Some 
additional  factors  have  to  be  considered  in  the  choice  of  both  test¬ 
ing  functions  and  basis  functions  :  (1)  the  accuracy  of  the  solution; 
(2)  the  ease  of  evaluation  of  the  "kernel"  matrix's  elements,  A^; 
and  (3)  the  size  of  the  matrix  A  and  its  behavior,  whether  it  is 
well-conditioned  or  not. 

For  numerical  computation  of  capacitances  using  the  capacitance 
model  in  See.  2.1,  the  basis  functions  and  the  testing  functions  are 
chosen  the  same  as  "pulse"  functions  which  are  constant  over  the 
subareas  on  the  conductors.  In  other  words,  the  surfaces  of  all  the 
conductors  are  divided  into  subareas  and  the  charge  densities  over 
these  subareas  are  assumed  to  be  constant.  Then  the  above  procedure 
is  adopted  to  solve  for  the  unknown  coefficients.  The  reasons  for 
using  pulse  functions  are  as  follows  :  (1)  they  are  linearly  indepen¬ 
dent  of  each  other;  (2)  the  charge  density  function  can  be  reason¬ 
ably  well  expressed  in  terms  of  them,  provided  that  the  subdivisions 
are  refined  enough;  (3)  a  closed-form  expression  can  be  derived  for 
the  double  integral  in  the  evaluation  of  the  kernel  matrix  elements; 
(4)  if  the  basis  functions  and  the  testing  functions  are  *he  same, 
the  kernel  matrix  is  not  only  symmetric  but  also  positive  definite. 


and  this  guarantees  a  well-behaved  kernel  matrix.  In  order  to  more 
cloaely  approximate  the  charge  density  function  with  the  choice  of 


pnlse-type  basis  functions,  variable  width  pulses  are  also  used.  It 
improves  the  accuracy  compared  with  that  of  constant  pulses  with  the 
same  number  of  subdivisions.  In  practice,  those  variable-width 
subareas  are  chosen  according  to  the  roots  of  Chebyshev  polynomi¬ 
als  [16] .  because  the  charges  tend  to  congregate  at  the  edges  of  the 
conductors,  and  Chebyshev  polynomials  have  smaller  widths  between 
roots  at  both  ends  of  the  interval  [-1,1].  Additionally,  an  explana¬ 
tion  [10]  based  on  expressing  the  kernel  matrix  resulting  from 
Galerkin's  method  as  a  quadratic  form  of  energy  proves  that  it  is  a 
positive  definite  matrix. 

1.1.1.  Integral  method  apd  even-odd  mode  analysis 

In  this  section,  we  show  the  incorporation  of  the  integral 
method  with  the  even-odd  mode  analysis  in  the  case  of  a  multi¬ 
conductor.  Let  us  continue  with  the  three-conductor  example  given  in 
See.  2.3.2.  Aasume  that  the  number  of  subareas  on  these  three  con¬ 
ductors  are  n1#  n»,  uj ,  so  that  the  total  number  of  subareas 

N  “  ai+n*+n, ,  then  Eq.  (2.12)  can  be  rewritten  as 

nx  ni+ua  N 

J  °i*»i(r’ )  +  J  «*•*(?)  +  J  aiWi(?) 

i-1  i-nx+l  i-Oi+Uj+l 

Since  Galerkin's  method  is  esqiloyed,  the  testing  functions  are  the 

same  as  the  basis  functions,  t.(^)  •  m.(r),  j-l,2,***,N,  and 

J  J 

Eq.  (2.14)  becomes 


s  •*-. 


Nr)  I  »j(r)dr  -  2  drJ  dr'<-  ^r'JGCr.r'U  U) 

S*  i«l  si  si  J 


tti+ni 


+  5  °ij  drJ  d?-i(^')G(J. r*)«  (P) 

Si  S»  J 


inx+l 


N 

+  I  drf  dr,*»l(r,)G(r,rf  )m4(r) 

i-ux-4+l  S*  S»  J 

j  ■  1 #2 » • • • #n1 


9(r)l  Mj(r)dr  «  J  “if  drf  dr'»4(r’ )G(r,r* )« . (J) 
S-  J  .  .  lJS4  JSX  J 


ai+n* 


+  $  aiJ  drJ  dr'“i(r' )G(r,r* )» .  i 

Sa  Sa  J 


i«ax+l 


+  5  °ij  drJ  dr'o^r' )G(r,r*  )«•  (r) 

i«n4+ii4+l  S*  S* 
j  m  nx+l» • • • ,nx+a4 


><r)l  (dj(r)dr  ■  J  oj  drJ  d^'w^r’ >GU.r' )«  (J) 

S*  i«l  si  Sx  J 


nx+n» 

+  ^  aiJ  drf  dr'«i(r’)G(r,r')<i».(r) 

i«nx+l  S*  S* 
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*  f 

i-ni+0,+1 
j  “  “x+ni+l, • 


aiJ  ^rJ  dr,*»|(r' )G(r*r#)»j(r) 

S,  S, 

..N 


Sx,  s*  and  S«  arc  respectively  the  areaa  of  the  three-conductor 
surfaces.  Note  that  for  different  potential  patterns,  only  the 
potentials  at  the  conductors  are  different,  i.e.,  only  0(r)|c 

0<r)ls>.  and  0(r)l§  in  the  above  expressions  are  changed.  This 
corresponds  to  only  changes  in  the  right-hand  side  vector  b  ot  the 
aatrix  equation  Ax  -  b.  As  an  example,  for  the  even  and  odd  modes 
specified  in  Sec.  2.3.2,  the  b  vectors  are 


b  * 


yS/nx+i 

^ Sj^nj+uj+l 


for  the  even  mode. 


and  for  the  three  odd  modes. 


ysx“* 

• 

Js** 

• 

Zfs*x 

• 

*>X  - 

• 

V»i*i 

• 

b,  - 

• 

I/Sx-nx+i 

• 

b,  - 

• 

ySx“nx+l 

e 

• 

J  S,“nx+nx+i 
• 

• 

y S,“nx+nx+i 
• 

e 

^•““x+ns+l 

e 

*  e 

• 

• 

On  the  other  hand,  the  kernel  matrix  A  only  concerns  the  Green's 
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fraction  and  the  basis  fractions,  and  it  stays  the  sane  for  different 
even-odd  node  excitations.  This  is  one  of  the  inportant  advantages 
of  the  integral  nethod.  Because  the  evaluation  of  the  kernel  natrix 

elenents  always  is  the  nost  computationally  intensive,  we  nave  to 
compute  it  only  once,  and  UJ  factorize  it  only  once,  then  use  back¬ 
ward  substitution  to  find  the  charge  densities  corresponding  to  bx, 
b3,  and  b,  as  listed  above.  If  the  differential  method  is  used,  we 
have  to  essentially  solve  four  boundary— value  problems,  and  cannot 
take  advantage  of  the  even-odd  mode  excitations  as  above.  From  this 
point  of  view,  the  integral  method  is  better  in  the  case  of  multiple 
conductors. 

Derivation  of  Two-Dimens  tonal  Green*  s  Fraction 

Two-dimensional  Green's  functions  to  Poisson's  equation  in  s 
homogeneous  medium  or  in  two  layers  of  media  have  been  derived  in 
[4,14],  and  they  are  given  here  for  reference;  on  the  other  hand,  we 
will  concentrate  on  the  derivation  of  the  three-layer  case  in  this 
section. 

asdiaa 

Consider  a  unit  line  charge  at  (x',y')  in  a  homogeneous  medium 
with  dielectric  constant  s,  and  the  ground  plane  is  at  y  *  0.  Then 
the  Green's  function  is 


\  «*,  •  . 
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-1  (x-x')  +(y-y') 

G(x,y|x'  ,y' )  -  -  la  - r - —  (2.15) 

2ne  (x-x')  +(y+y') 


Not*  that  the  logarithmic  function  in  the  above  expression  character¬ 
ises  a  two-dimensional  line  source,  and  the  complementary  terms  of 
y— y'  and  y+y'  reflect  the  consequence  of  the  source  and  the  image 
across  the  ground  plane. 


£.£.£.  Two  layers  of  media 


Let  there  be  two  layers  of  media  over  the  ground  plane  as  shown 
in  Fig.  2.4.  The  first  layer  has  a  thickness  h  over  the  ground  plane 
and  a  dielectric  constant  sx.  The  second  layer  on  top  of  the. first 
one  has  a  dielectric  constant  e3.  There  are  four  subfunctions  asso¬ 
ciated  with  the  Green's  function  depending  on  the  relative  positions 
of  source  point  (x',y')  and  observation  point  (x,y)  namely: 
Gixtx.ylx* ,y' ) ,  Gax(x,y|x’ ,y' ) ,  Gsa(x,y|x* ,y* )  and  Gxa(x,y|x' ,y' ) . 
The  first  digit  of  the  subscript  specifies  the  region  where  the 
observation  point  resides,  and  the  second  digit  indicates  the  posi¬ 
tion  of  the  source  point.  This  convention  will  also  be  used  for 
other  subfunctions  of  other  Green's  functions  discussed  later.  These 
four  subfunctions  have  been  derived  by  the  method  of  multiple  reflec¬ 
tions  across  the  ground  plane  and  the  interface  between  the  two  media 
[4].  They  are  listed  as  : 


7  1 

LAYER  1 

f  (x,y)  1 

' 

€x 

Gn<x.yU'*y') 


l 


«»“8l 


4JI8J 


■5. 

k-0 


*i+«i 


-)*  In 


(x-x')  +(y+y'+2kh) 

{ x-x ' ) *  +( y-y' +2kh) * 


1  ?  «»-«i  (x-x')a+(y+y'-2(k+l)h)a 

+ - 2  < - )k  1  .ln  - - — 

<*«x  t«o  **+#*  (x-x')*+(y-y*-2(k+l)k)* 

0  <  y  i  h  and  0  i  y'ik  (2.16e) 


G»(x.yli'.y') 


1  v  «i-«x  V  (x-x')1+(y+y'+2kJi)* 

-  2  ( - >k  In  - 

4n#*  k-0  **+,x  (x-x')a+(y-y'-2kb)* 


1  £  «a-«i  k+1  (x-x')a+(y+y'+2(k-l)k)* 

- 2  < - )k  1  ln - — - 

k-0  *»+,i  (x-x')a+(yy'-2(k-l)k)a 

0  <  h  <  y  «d  0  <  li  i  y'  (2.16b) 


Giatx.ylx*  »y’) 


l 


2ji(«,+«t) 


■  2  < - 

]t»0  *j+*i 


in 


0  i  y  i  i>  i  y' 


( x-x ' )  *+(  ry’  +2kb) a 

(x-x')*+( ly-y* I +2kh) 
(2.16c) 


G»x<x.y|x' ,y»)  -  G13(x,y|x' ,y* )  except 


0  1  y’  i  k  £  y 


(2.16d) 


Compared  with,  the  formula  in  Eq.  (2.15),  the  above  formulas  basically 
are  still  logarithmic-type  functions  with  some  modification  due  to 
the  two  media.  The  infinite  summation  in  the  formula  is  a  direct 


result  of  the  infinite  reflections  of  images  between  the  two  infinite 
parallel  planes. 

£*5.*sL*  Three  levers  of  media 

The  method  of  multiple  reflections  employed  in  [4]  to  derive  the 
Green's  function  for  two  layers  of  media  is  too  complicated  to  formu¬ 
late  the  three-layer  case.  Hence,  another  method  of  boundary  match¬ 
ing  in  the  spectral  domain  [19]  is  employed  instead. 

As  shown  in  Fig.  2.5,  the  dielectric  constants  of  the  three 
layers  of  media  are  aa,  sa,  and  s0.  The  outermost  layer  is  assumed 
to  be  free  space  always,  and  the  interconnect  will  not  reside  in  this 
region.  Thus,  there  are  four  subfunctions,  Gxx,  Gax,  and  Gxa,  Gaa, 
to  be  considered.  The  convention  on  the  subscript  of  those  subfunc¬ 
tions  is  the  same  as  that  discussed  in  the  previous  section.  Let  us 
consider  the  case  when  source  point  (x',y')  is  in  the  layer  ex,  and 
derive  the  subfunctions  Gla(x,y|x'  ,y» )  and  Gia(x,y|x' ,y’ )  .  The  other 
two  subfunctions  can  be  derived  similarly. 

As  depicted  in  Fig.  2.5,  an  artificial  interface  is  drawn  past 
the  source  point  (x* ,y')  in  layer  es,  which  increases  the  number  of 
regions  to  four.  In  analogy  to  Fourier  analysis,  assume  that  the 


potentials  in  these  four  regious  are  : 
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Not®  that  the  a^k)'*  and  Pi(k)'i  in  Eq.  (2.17)  may  be  considered  ss 
the  unknown  spectral  components  of  the  corresponding  potentials. 
Because  of  the  ground  plane  in  region  1,  potential  0x(x,y)  should  be 
zero  at  y  ■  0  and  does  not  contain  a  cosh(ky)  term.  For  the  outer¬ 
most  region  4,  e  **  represents  a  decaying  wave  in  the  open  boundary. 

At  the  interfaces  between  different  regions,  the  potentials  and 
the  corresponding  "displacement  vectors"  D's  have  to  be  continuous. 
Hence,  the  boundary  conditions  can  be  written  as  follows  : 


(2.18b) 


d0i 


6<y-y' ) 


(2.18d) 


90i  90, 

**  9y  '  #*  9y  W 


**»  lyb  “  ly-b 


90,  t  904  , 

«.  'y-b  “  ••^“ly-b 


ay 


ay 


(2.18o) 


( 2 . 1 8  £ ) 


Note  that  these  six  boundsry  conditions  uniquely  determine  the  six 
unknown  spectral  components,  a^kj's  snd  fl^k^s  in  Eq.  (2.17). 
After  some  tedious  algebraic  manipulation,  the  spectral  components 
can  be  solved  : 


I 

°i(k) 


2 


«(«,+«,) 


-1—  .-k=,,<i-«,.-2l(b-’r')> 

kD(k) 


(2.19a) 


°,(k) 


— - - -  e"ky ' ( l-$ ,co  sh ( 2ka ) ) ( l-$ , e"2k (b”y ' * ) 

ns,  kD(k)  (2.19b) 


P,(k) 


— - -  e_ky’  ((,sinh(2ka) ) (l-5,e“2k(b"y' ) ) 

ns. 


kO(k) 


(2.19c) 
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«,(k)  -  — - - - - <8inli(ky')+51*iah(k(2*-y,)))(l+{1e_2kb) 

kO(k)  ( 2 .19d) 


1  1 

P»  (k)  - - ( sinMky* )  +$aainh(k(2a-y* ) ) )  (l-?ae-Zkb) 

n«a  kD(k)  (2.1»e) 


2  1 

a4(k)  - - (sinMky' )+$xsinh(k(2a-y* ) ) )  (2.19f) 

*<«!+*»)  kD(k) 


«s“«i  «•-«» 

where  5X  -  -  sad  5a  -  - 

•s+«i  «e+«s 


«»d  D(k)  -  1  -  ?xe-2k*  '  «ae_2kb  +  ?x5ar2k(b_a) 


Substituting  Eq.  (2.19)  back  into  Eq.  (2.17),  vo  can  obtain  the 
potentials  in  spectral  integral  forai  analogous  to  Fourier  integrals. 
Sinee  the  source  in  the  above  derivation  is  in  the  layer  ea,  the 
potentials  0a  and  0, ,  both  in  layer  es  as  a  result  of  the  artificial 
interface  introduced  in  the  derivation,  can  be  coabined  to  fora  the 
subfunction  Gaa.  0a  is  naturally  the  subfunction  Gaa. 

Vritten  in  spectral  integral  fora,  0a,  as  an  ezaaple,  can  be 
expressed  as 


•»(x.y)  "  - f 

2nsa 


cos(k(x-x* ) ) 


■(i-?,e-2k(b-y'>) 


•-  -vv:-.vv 


-  *•  -  -  •  •.•aAn.1 1^“  a.'  O  JVJk  JVA  , 


kD(k) 
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( 0"k( 7' -y) _#-k( y' +y) i#-k(y' +y-2a) 

b  >  7*  >  y  2  » 


-«1e”k(y'-y+2*>)dk 

(2.20) 


la  order  to  "reverse  Fourier  transform"  Eq.  (2.20) ,  or  integrate  ont 
the  dummy  variable  k,  first  the  factor  [D(k)]-1  is  expanded  in  bino¬ 
mial  form  as  : 


[DU)]"*  -  5  5  1  (-l)k<i)(k)  «*“  *“2kbn  e”2k*(i-»“k) 

j-0  n-0  k-0 

Then  sobstitnting  this  binomial  expansion  into  Eq.  (2.20)  utiliz¬ 
ing  the  following  definite  integral  from  [17]  : 


o-7*_e-0x  „V 

I  - -  1  oos(bx)  dx  -  (1/2) In 

x 


.»  * 

b  +t r 


Rep>0,  Rer>0 


we  can  integrate  out  the  spectral  variable  k  and  express  f>a(x,y)  in 
the  space  domain  as  : 


1  r  J  “ 

•»(x.y)  - Y  Y  Y  (-Dk(i)(j)«xj"n+k?»n 

4rt,»  j-o  n-o  k-0 

_  (x-x')*+(y+y,+2bn+2a( j-k-n))* 

[  ln - ; - ; — 

( x-x ' )  +(-y+y'+2bn+2a( j-k-n) ) 


+  ?!  In 


(x-x')  +(-y+7*  +2bn+2a( j-k-n+1) ) 


(x-x' ) *+(y+y' +2bn+2 a (j-k-n-1) ) * 


(x-x')4+(-y“y*  +2b(n+l)+2a(  j-k-n)  )* 


(x-x' ) *+(y-y' +2b(n+l) +2a( j-k-n) ) * 


( x-x ' )  * +( ry’  +2b  ( n+1) +2a ( j-k-n-1) ) 4 

+  In  - ; - - 

<  x-x' )  +( -y-y' +2b (n+1) +2a( j-k-n+1) ) 


0  <  ‘  i  7  i  7'  i  b  (2.21) 


Because  tbs  constraint  on  0,(x.y)  i*  »  j  y'  i  y  <  b,  and  0,  baa  a 
similar  formula  as  that  in  Eq.  (2.21),  tbs  combination  of  0,  and 
results  in  substituting  absolute  values  for  y-y'  in  tbe  above 
Eq.  (2.21) .  Tbe  other  subfuactions  can  be  obtained  in  a  similar 
manner.  In  summary,  ve  have  tbe  two-dimensional  Green's  function  in 
three  layers  of  media  as  follows  : 


l  *  J  a 

®a*(x,ylx',y»)  - J  }  (-l)k(i)(J)?1j"n+k?,n 

j-0  n-0  k-0 

.  (x-x')*+( ly+y' l+2bn+2a( j-k-n) ) * 

[  in - - - 

(x-x')  +( lyy'  I  +2bn+2a(  j-k-n) ) * 
(x-x')4+(“ly+y' l+2bn+2a( j-k-n+1) )* 

+  in  - - - 

(x-x')  +(-|y-y' l+2bn+2a( j-k-n+1) ) 


Jl 
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+  In 


( x-x ')*+(- I y+y' I +2b (n+1) +2a( j-k-n) ) * 
(x-x' ) *+(-|y-y' |+2b(n+l)+2*( j-k-n) )* 


( x-x' ) *♦< ly+y' I +2b(n+l) +2»<  j-k-n-1) ) * 

+  1®  ““ "  3  '  "  a  ] 

( x-x ' )  +( I y-y' I +2b (n+1) +2»( j-k-n-1) ) 

0  <  7.7#  i  *  <  b  (2.22.) 


®»x(n*yl*' .y')  *  Gia(x,y|x' »y' )  except 


°<y'  1  *lyib 


(2.22b) 


G»a(n.yl*' .y') 


—  255 

4**»  j-0  a-0  k-0 


[  lu 


(x-x')*+( ly+y’ l+2bn+2»( j-k-n) )* 
(x-x')*+( I y-y' l+2bn+2a( j-k-n))* 


+  ?i  In 


♦  In 


(x— *')*+( ly-y' |+2bn+2*( j-k-n+1) )* 
( x-x' )*+( ly+y' I +2bn+2e( j-k-n-1) ) * 


(x-x')‘+(-ly+y' l+2b(n+l) +2*( j-k-n) ) * 
(x-x' ) *+(-|y-y' l+2b(n+l)+2*( j-k-n) ) * 


V 


-1 


.1 
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E 


V 


(x-x')*+(“ly-y’  I  +2b ( n+1)  +2* ( j-k-n-1) ) * 

+  5l5i  ln  - - - - - 

(x-x')  +(-Iy+y'  I +2b(n+l) +2a( j-k— n+1) ) 


0  <  *  <  7#y'  i  b 


(2.22c) 


j  n 


Gis(x.yU'  ,y') 


2*(«x+s,)  j«0  n-0  k-0 


III 


[  1, 


(x-x')*+( ly+y' l+2bn+2a( j-k-n) ) * 
( x-x' )  *  +(  | y-y'  |  +2bn+2a ( j-k-n) )  * 


(x-x')a+("ly+y# l+2b(n+l)+2a( j-k-n))* 

+  In  .  '  1 

(x-x')  +(-|y-y* |+2b(n+l)+2s( j-k-n)) 

0  1  y  i  «  1  y'  i  b  (2.22d) 


V. 


Note  tbet  these  expressions  still  retain  the  basic  logarithmic 
function  as  that  in  Eq.  (2.15)  .  Tith  the  pnlse  basis  function,  the 
double  integral  in  the  evaluation  of  the  kernel  elements  will  result 
in  a  closed  form.  The  detailed  expressions  of  these  closed-form 
evaluations  are  included  in  Appendix  C.  Furthermore,  the  complex 
triple  summation  in  the  above  formulas  can  be  interpreted  as  the  mul¬ 
tiple  reflections  among  the  three  parallel  planes,  i.e.,  the  ground 
plane  and  two  dielectric  interfaces.  The  program  which  executes  the 
aforementioned  procedure  of  the  integral  method  is  named  "CAP2D*  and 
a  user's  manual  for  it  is  included  in  Appendix  B. 


F  » 
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*  -■  -*  •'*„ 

Jm.  -h  m  haJfc*  a 


■  ad 


.  %  "w  * 

■V 
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2,.£.  Derivation  of  Three-D imens tonal  Green's,  Fnnction 

There  are  soae  circumstances  in  which  the  two-dimensional 
analysis  cannot  satisfactorily  compute  the  interconnect  capacitances. 
For  example,  the  cross-over  of  two  perpendicular  interconnects  on 
different  levels  is  one  such  example.  Hence,  the  three-dimensional 
Green's  functions  in  homogeneous,  two-layer  and  three-layer  media  are 
examined  in  this  section. 

£.£. i,.  Homogeneous  medium 

Consider  a  unit  point  charge  at  r'  ■  (x',y',z')  in  a  homogene¬ 
ous  medium  with  dielectric  constant  a  and  the  ground  plane  at  z  ■  0. 
The  Green's  function  to  the  Poisson's  equation  satisfying  this  boun¬ 
dary  condition  has  been  derived  16]  and  is  written  as  : 


GU,7')  UA(*-*')V1/2  -  [p1+(z+*')1]”1/2}  (2.23) 

4na 

where  p*  ■  (x-x' ) 1+(y-y* ) * 

Note  that  r  *  (x,y,z)  and  r'  ■  (x',y',z')  are  the  observation  and 
the  source  point,  respectively.  The  core  of  the  expression  is  a 
function  proportional  to  the  reciprocal  of  the  distance  which  is  dif¬ 
ferent  from  the  logarithmic  function  in  the  two-dimensional  case. 
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£.£.£.  Two  layers  of  media 

The  basic  principle  need  in  the  derivation  of  the  three- 
dimensional  Green's  function  in  multiple  layers  of  media  is  similar 
to  boundary  matching  in  the  spectral  domain  as  discussed  in  Sec.  2.5 
for  the  two-dimensional  case.  The  formulation  had  been  introduced  in 
[fij .  For  two  layers  of  media,  let  the  dielectric  constant  of  the 
lower  layer  be  Sj,  and  the  thickness  be  h.  Assume  that  the  ground 
plane  is  at  z  *  0,  and  the  dielectric  constant  of  the  top  layer  is 
•a.  Utilizing  the  same  convention  as  the  two-dimensional  Green's 
function,  the  three-dimensional  Green's  function  for  two  layers  of 
media  is 


-  - - - - -  J  Sn[lP*+<*+*'+2ak>a]"1/2 

4naxr  4ns t  B<p,0 

-  ?tp*+(z-z'+2(n+l)h)*r1/2  -  $[p*+(z-z'-2(n+l)h)ar1/2 

+  ?[p1+(z+z'-2(n+l)h)*]"1/2] 

0  <  z.z'  <  h  (2.24a) 


Gas(***’) 


1 


4ns#r 


— -  5  Sa[-?[pl+<*+*'+2(n-l>10  V1/2 

4R,a  a-o 


1 


+  [p’+U+z^nh)’]-1/2] 


0  <  h.  <  *,*• 


(2.24b) 


1  2t%  m 

,(r.r*)  -  - ( - — )  J  5*[[pa+(*-*»-2ah)*]“1/2 

4ff*»  «a+«x  n„0 

-  [p*+(z+z'+2nh)‘]"1/2] 

0  <  z  <  h  i  z' 


(2.24c) 


1  2«  ™ 

*x<*»*'>  -  - ( - )  3  5nf[p*+(z-z'+2nh)*]_1/2 

Am  m  -  -l.  ~  *• 


4n*x  *a+*i 


-  [p1+(z+z'+2nh)*]“1/2] 

0<»-ai: 


(2.24d) 


where  p  -  (x-x' )  *+(y-y' )  *  ,  r*  -  p*+(z-z')*  and  $  «  -*  *3 


Three  layers  o£  media 

Consider  the  configuration  in  Fig.  2.$.  The  thicknesses  of  the 
layers  are  a  and  b-a  and  the  dielectric  constants  are  e#,  e,,  and  ex. 
These  parameters  are  the  same  as  those  in  the  two-dimensional  case. 
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Utilizing  one  well-known  Fourier  integral  [15,6] : 

r"1  -  p  J#(k p)  •  dk  (2.25) 

where  r  *  p  +  (z-z*)  ,  and  J#(x)  is  the  Bessel* s  function  of  first 
kind  order  0,  we  can  derive  the  Green's  function. 

First,  consider  the  case  in  which  the  source  point  resides  in 
layer  t,  and  the  subfunctions  Gia  and  Gaa  are  derived.  As  before, 
the  potentials  in  layers  s#,  ea,  ea,  are  represented  as  Fourier 
integrals: 


0,  .  -1 -  f"  Q(k)J0(kp)e_k(x"z,>  dk 

4n«, 

0  <  a  <  z*  ibiz  (2.26a) 


*  Gaa(r,r*) 


«• 

— — )  J0(kp)e 


-klz-z'  I 


dk 


+ 


f  Ox(k) J0 (kp) e”k<  z~z> 


dk 


+  J^8a(k)J.(kp)ek(z_z')dk] 

0  <  *  <  z,  z'  £  b  (2.26b) 

1  ® 

“  Gia(r,r* )  -  — - - ff  T1(k)J,(kp)e_k(z"z')dk 

4ne,  LJ* 
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+  J#Tl(k)J#(kp)ekU”z')<ik] 

0  <  z  £  a  and  i  1  i'  £,  b  (2„26c) 

%  2  % 
where  p  *  (x-x*)  +  (yy1) 

The  functions  Q(k) ,  <J,.(k) ,  03(k),  Ti(k),  and  ti(k )  are  the  unknown 
spectral  components  to  be  determined  by  boundary  matching.  As  in  the 
two-dimensional  case,  the  exponential  term  e~fc(z~z')  in  0o  dictates  a 
decaying  wave  towards  the  open  boundary  in  layer  s0;  also,  the  abso¬ 
lute  value  on  z-z*  is  removed  because  z  is  always  greater  than  z'  in 
this  region.  Note  that  the  first  term  in  0,  accounts  for  the  singu¬ 
larity  of  the  source  point  in  the  same  region.  Because  the  relative 
magnitudes  of  z  and  z*  are  not  fixed  in  this  region,  the  absolute 
value  is  needed  in  the  exponent. 

As  in  the  two-dimensional  case,  the  continnity  of  the  potentials 
and  the  displacement  vectors  dictates  the  following  necessary  boun¬ 
dary  conditions  : 


“•'.-o  ■  0 


30i  d0s 

**~T“U-a  “  8*”T  U-a 

3z  3z 


(2.27a) 


(2.27b) 


(2.27c) 


(2.27d) 


1  z«b  *  z=b 


30»  30. 

z-*>  “  z«b 

3z  dz 


(2.27 t) 


Substituting  Eq.  (2.26}  into  Eq.  (2.27).  we  can  solve  the  five  unk¬ 
nown  spectral  components  in  Eq.  (2.26).  After  some  algebraic  manipu¬ 
lation.  they  are  listed  as  follows  : 


*a 


(2.28a) 


0, 


(k)  -  "-?,-~A_1[e"2k(*“z')(l-51e-2k*)  +  (5x-e”2k*)]  (2.28b) 


2*. 


xx(k)  -  --^L.A-le-2k.(5i_e2k(b-z') 


(2.28c) 


2s. 


*XU>  ~-A-1.-2k*(e2kM1e2kz': 

•x  ■*’*s 


(2.28d) 


where  A  -  -  {l,2Mb-2.)  +  ,2‘<‘-‘) 
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Note  that  the  function  Q(k)  is  not  associated  vith  the  Green's  func¬ 
tion  in  question,  so  it  is  not  listed  above.  Utilizing  the  binomial 
expansion  of  A-2  and  the  following  definite  integrals  [17]  : 

f*  •”oxJ#<0x)dx  -  (a*  +  Re(a+i0)  >  0 

J  0 

r  e_axcoshOx)J0(rx)  -  O.SUyVp-a)1]"*7*  +  tr'+O-o)1]"17*}, 

Rea  >  lRep|,  y  >  0 

we  can  transform  the  Fourier  integrals  in  Eq.  (2.26)  back  into  the 
space  domain  in  order  to  obtain  the  two  subfunctions  of  the  three- 
dimensional  Green's  function  in  three  layers  of  media  as  follows  : 


G*» (r ,r' ) 


1 

4ns,r 


— - — }  y  j  (-i>n  k(i)(j)5xn5»j'k- 

4r«i  j ®0  n»0  k*G 


[[p*+(2((j-k)b+(2k-u)a)+z+z')*]~1/2 
-  ?1tp‘+<2((j-k)b+(2k-n-l)a)+i+z'),]_1/2 
"  ?,{[p1+<2((j-k+l)b+(2k-n)a)-z+z'))*]_1/2 
+  [p*+(2( ( j-k+l)b+( 2k-n) a)+z-z' ) ) *]”1/2 
-  [p*+(2( ( j-k+1) b+( 2k-n) a)-z-z' ) )*]_1/2} 

+  515,{[p*+(2((j-k+l)b+(2k-n-l)a)-z+z')),]'1/2 


+  [p*+(2<(j-k+l)b+(2k-n-l)a)+z-z'))*r1/2 

-  [p*+(2( ( j-k+l)b+(2k-n+l) a)-z-z' ) ) *]”1/2} ] 

0  <  »  <  :  i  b  «d  0  <  a  1  z'  1  b  (2.29a) 


Gu(r,r') 


1  l  5 


2”(*i+*a)  j«0  a-0  k-0 


[{[p*+(2((j-k)b+(2k-n)a)+|z-z* i))*]“1/2 
-  [p*+(2( ( j-k)b+(2k-a) a)+| z+z' |))*]"1/2J 


-  ?,{[p*+(2( (j-k+l)b+(2k-n)a)-lz+z' I))*]"172 


+  [p*+(2( ( j-k+l)b+(2k-a)a)-|z“z' I))*]"172)] 


0  <  z  £  a  <  b  and  0  <  »  i  z'  1  b 


(2.29b) 


The  other  two  subfunctions  of  the  Green's  function  can  be  similarly 
formulated  and  they  are 


4ns 


—  I  $  \  <-i)i<i>(i)s»*t.J**fk- 


ir  *n*i  j»0  n-0  k«0 


[lp*+(2(-( j-k) z' +( j-n+k)b+(n-2k)a)+z+z' ) ) *]”1/2 
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+  ?i5»[p*+(2(-(j-lt)i'+(  j-a+k+l)b+(n-2k-l)a)+z+z'))a]"1/2 

-  ?! { [ p*+< 2< -< j-k) z' +( j-n+k)b+< n-2k+l)  a)  +z-z' ) )  * ]~1/: 2 

-  [p1+(2(-(j-k)z'+(j-n+k)b+(n-2k+l)a)-z-z'))1]“1/2 

+  [p1+(2(-(j-k)z'+(j-n+k)b+(a~2k+l)a)-z+z'))*]“1/2} 

-  ?a{[p,+(2(-(j-k)z'+(j-a+k+l)b+(a-2k)a)+z-z'))1]-1/2 

-  [p*+(2(-(j-k)z'+(j-n+k+l)b+(n-2k)a)-z-z'))*]"1/2 

+  [ p*+<  2< -< j-k) z' +( j-n+k+1) b+( a- 2k) a ) -z+z ' ) ) * ] ~1  ( ' 2) ] 

0  <  z  <  a  <  b  a ad  0  £  z'  £  a  <  b  (2.29c) 

Gii(r»r')  ■  Gxa(r,r')  except 

0  <  z'  £  a  <  b  and  0  <  a  £  z  £  b  (2.29d) 

It  should  be  poiated  oat  that  the  "core"  function  in  the  above 
expressions  is  still  the  inverse  square  root  function  mentioned  ear¬ 
lier.  Associated  with  the  pulse  basis  functions  in  Appendix  D,  the 
inverse  square  root  function  results  in  closed-form  integration  for 
the  evaluation  of  the  kernel  matrix  elements.  A  detailed  derivation 
of  these  kernel  integrals  is  included  in  Appendix  D.  The  program 
vhich  implements  the  numerical  computation  of  the  three-dimensional 
interconnect  capacitances  described  above  is  named  "CAP3D",  and  the 


user's  manual  for  it  is  the  same  as  the  one  for  "CAP2D",  and  is 
included  in  Appendix  B  with  soae  examples. 


2.2*  Some  Results  of  Two-  and  Three-Dimens ional  Numerical  Computa¬ 
tion 

2*2 *i*  Comparison  of  two-  and  three-dimens ional  capacitance 

In  two-dimensional  capacitance  computation,  the  length  of  the 
interconnect  is  assumed  to  be  infinite;  in  other  words,  the  calcu¬ 
lated  capacitance  is  the  capacitance  per  unit  length  of  the  intercon¬ 
nect.  In  comparison  with  the  three-dimensional  capacitance,  as  shown 
in  Fig.  2.7,  the  capacitances  are  plotted  versus  the  "normalized” 
length.  Since  the  two-dimensional  capacitances  and  the  parallel- 
plate  capacitance  are  proportional  to  the  length,  they  are  linear 
lines  in  Fig.  2.7.  As  expected,  the  discrepancies  between  the 
corresponding  three-dimensional  and  two-dimensional  capacitances  are 
diminishing  as  the  length  increases.  Some  of  the  difference  may  be 
attributed  to  the  "end  effect"  or  the  fringing-f ield  capacitances  at 
the  ends  of  the  interconnect  which  are  included  in  the  three- 
dimensional  computation.  Also,  the  difference  between  two- 
dimensional  capacitance  and  the  parallel-plate  capacitance  is  due  to 
the  fringing-f ield  capacitance  at  both  sides  of  the  interconnect 
which  had  been  taken  into  account  in  the  two-dimensional  computation. 


Two  cases  of  thickness  are  also  considered  in  Fig.  2.7.  In  the 
case  of  t/h  -  0,  the  interconnect  is  considered  infinitesimally  thin. 


NORMALIZED  LENGTH  (  l/h*Icm  ) 


A  :  3-d  and  t/h  ■  0.8 
B  :  2-d  and  t/h  *0.8 
C  :  3-d  and  t/h  -  0 
D  :  2-d  and  t/h  "  0 
E  :  parallel-plate  capacitance 

*  Capacitancea  are  in  nnits  of  8.854xl0~14e  Farad,  if 
the  length  is  in  units  of  centimeter,  where  e  is  the 
dielectric  constant  of  the  homogeneous  medium. 


Figure  2.7  Comparison  of  two-  and  three-dimensional  capacitance 
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end,  naturally,  the  fringing-f ield  capacitances  from  the  side  walls 
are  reduced  and  all  the  corresponding  capacitances  are  smaller  com¬ 
pared  with  those  of  thickness  t/h  >0.8.  From  this  comparison,  it 
can  be  concluded  that  as  the  dimensions  of  interconnects  are  reduced, 
the  fringing-f ield  capacitance  does  not  decrease  in  the  same  propor¬ 
tion  as  the  parallel-plate  field  capacitance.  Thus,  the  fringing 
capacitance  gradually  emerges  as  the  dominating  component  in  the 
evaluation  of  the  capacitances.  Hence,  more  accurate  computation  is 
needed  to  evaluate  them  in  order  to  take  all  the  capacitance  com¬ 
ponents  into  account. 

2.2.2.  Coupling  capacitance  and  slant  side  walls 

Shown  in  Fig.  2.8  are  the  cross-sections  of  two  interconnects. 
The  self-  and  coupling  capacitances  with  respect  to  the  separation 
"s"  and  the  pitch  angle  of  the  side  walls  "a"  are  given  in  Table  2.2. 
First  of  all,  it  can  be  deduced  from  the  table  that  the  coupling 
capacitance  is  inversely  proportional  to  the  separation,  and  the 
self-capacitance  decreases  as  the  coupling  capacitance  increases. 
This  effect  may  be  attributed  to  the  fact  that  the  fringing  fields  on 
the  side  wall  are  directed  more  towards  the  adjacent  interconnect 
inatead  of  the  ground  plane  as  the  separation  decreases.  Conse¬ 
quently,  some  fringing-f ield  capacitance  in  the  self-capacitance 
shifts  to  the  coupling  capacitance.  However,  as  compared  with  the 
differences  in  capacitances  shown  in  Table  2.2,  the  coupling  capaci¬ 
tance  increases  more  than  the  reduction  in  the  self-capacitance. 


Us- 
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TABLE  2.2  Self-  and  Coupling  Capacitance  vs  Pitch  Angle  and  Separation, 


*  pitch  angle  (a) 


a/h 

o 

o 

ye 

70° 

75° 

80° 

90° 

1.0 

22.332 
6.729  „ 
(4.266) 

22.667 

7.984 

(4.626) 

22.839 

8.801 

(4.813) 

23.018 

9.811 

(5.006) 

23.407 

13.236 

(5.426) 

2.0 

23.469 

4.145 

(3.129) 

23 .844 
4.830 
(3.449) 

24.036 

5.236 

(3.616) 

24.234 

5.706 

(3.79) 

24.660 

6.964 

(4.173) 

5.0 

25.111 

1.666 

(1.487) 

25.613 

1.922 

(1.68) 

25.870 

2.060 

(1.782) 

26.136 

2.209 

(1.888) 

26.710 

2.555 

(2.123) 

10.0 

25.968 

0.659 

(0.63) 

26.573 

0.759 

(0.72) 

26.885 

0.812 

(0.767) 

27.208 

0.866 

(0.816) 

27.911 

0.987 

(0.922) 

*  The  capacitances  are  in  units  of  8.854xl0~i4  Farad/cm, 

and  the  dielectric  constant  of  the  homogeneous  medium  is  4.0. 


••  The  decrease  in  the  self-capacitance  as  a  result  of  the 

adjacency  of  another  interconnect.  The  reference  capacitance 
is  the  self-capacitance  of  a  single  interconnect. 


Next  the  pitch  angle  seems  to  have  a  greater  effect  when  the 
separation  is  small.  For  example,  at  s/h  *  1,  the  coupling  capaci¬ 
tance  with  a  90°  pitch  angle  more  than  doubles  that  with  a  60°  pitch 
angle;  whereas  at  s/h  “  10,  the  difference  in  coupling  capacitances 
is  about  50%.  This  is  due  to  the  fact  that  the  closer  the  two  inter¬ 
connects,  the  more  strongly  coupled  they  are,  and  the  more  evident  is 


the  change  in  pitch  angle.  Since  the  design  rules  of  the  separations 
between  adjacent  interconnects  are  usually  dictated  by  the  considera¬ 
tion  for  aetal  migration,  it  may  not  be  necessary  to  consider  the 
pitch  angle  effect.  However,  this  ability  to  take  into  account  the 
pitch  of  the  side  walls  of  interconnects  is  helpful  if  detailed 
information  is  needed.  Subsequently,  in  all  further  capacitance  com¬ 
putation  of  interconnects,  a  90°  pitch  angle  will  be  assumed. 


In  the  channel  routing  problem,  there  is  some  concern  about 
whether  the  knock-knee  or  the  cross-over  configuration  will  have 
smaller  coupling  capacitance.  In  Fig.  2.9  and  Fig.  2.10,  both  the 
knock-knee  and  the  cross-over  configurations  of  two  interconnects  in 
different  levels  are  shown.  For  the  sake  of  simple  computation,  the 
thickness  of  interconnects  is  assumed  to  be  infinitesimally  thin.  As 
discussed  previously,  if  the  thickness  is  included,  the  coupling 
capacitance  will  be  increased  due  to  the  increased  side  wall  areas. 
The  calculated  capacitances  for  these  two  cases  are  listed  in 
Table  2.3  to  Table  2.6.  For  simplicity,  both  interconnects  are 
assumed  to  have  the  same  width  and  length.  The  parameter  "l"  is 
defined  comparably  in  both  knock-knee  and  cross-over  cases  as  shown 
in  Fig.  2.9  and  Fig.  2.10.  Caa  is  the  self-capacitance  of  the  lower 
interconnect,  and  Caa  is  that  of  the  upper  one.  More  important,  the 
coupling  capacitance  between  the  two  interconnects  is  Cla,  Comparing 
Table  2.3  with  Table  2.4,  or  Table  2.5  with  Table  2.6,  it  is  apparent 
that  either  the  knock-knee  or  the  cross-over  configuration  has  very 
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TABLE  2.3  Knock-Knee  Capacitance. 


(Wx/^l 

■  W,/llx 

-  1.  k. 

/kx  “  2) 

1/hx 

2.0 

3.0 

5.0 

10.0 

20.0 

Cxx 

6.65 

9.03 

14.20 

27.87 

56.00 

5.59 

7.32 

11.08 

21.22 

42.47 

Css 

2.06 

2.68 

3.42 

4.13 

4.52 

TABLE  2.4  Cross-Over  Capacitance. 
<wl/hi  »  Wj/hx  »  1,  h,/hx  ■  2) 


TABLE  2.5  Knock-Knee  Capacitance 


(Wj/ha  = 

Wj/hj  = 

1,  1/ba 

=  10.0) 

hj/ha-1  0.5 

1.0 

2.0 

5.0 

10.0 

Ctl  27.59 

27.87 

28.38 

29.41 

30.14 

Cj j  22.98 

21.22 

19.38 

18.49 

18.20 

Cla  5 .66 

4.13 

2.88 

1.41 

0.59 

TABLE  2.6 

Cross-Over  Capacitance. 

(Wj/ha  * 

wa/hx  * 

1.  1/ba 

-  10.0) 

h,/ha-l  0.5 

1.0 

2.0 

5.0 

10.0 

CX1  28.22 

28.61 

29.40 

30.75 

31.61 

C,a  23.56 

21.78 

20.41 

19.4-8 

19.27 

Cx,  6.71 

5.26 

3.70 

1.74 

0.69 

siailar  values  of  coupling  capacitance  between  the  two  interconnects 
in  different  levels.  Therefore,  it  does  not  matter  whether  the  chan¬ 
nel  routing  problem  is  solved  by  the  knock-knee  or  the  cross-over 
configuration  from  the  viewpoint  of  the  capacitance  penalty. 
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In  summary,  the  two-  and  three-dimensional  capacitance  computa¬ 
tions  can  be  applied  to  various  interconnect  configurations .  As  the 
dimensions  of  the  interconnects  and  the  proximity  between  them  become 
smaller,  not  only  the  self-capacitances  need  more  accurate  computa¬ 
tion,  but  also  the  coupling  capacitances  become  more  important.  In 
this  chapter,  the  numerical  computation  of  the  interconnect  capaci¬ 
tances  was  discussed.  Fortran  programs  "CAP2D”  and  "CAP3D"  were  also 
developed  to  implement  the  integral  method  discussed  here.  In  order 
to  easily  incorporate  these  capacitances  into  an  extraction  program, 
some  simpler  closed-form  formulas  are  needed  instead.  These  will  be 
discussed  in  the  next  chapter. 
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CHAPTER  3 

APPROXIMATION  FORMULAS  FOR  INTERCONNECT  CAPACITANCES 

2.-1.  Introdoct  ion 

As  the  complexity  of  the  present  VLSI  circuits  increases,  the 
interconnect  parameters  have  become  more  important  in  both  circuit 
design  and  circuit  verification.  In  Chapter  2,  numerical  approaches 
to  the  computation  of  interconnect  capacitances  were  discussed.  How¬ 
ever,  those  computations  are  usually  quite  involved  and  require  a 
large  amount  of  computation  time.  In  extraction  programs  for  design 
verification,  it  is  impractical  to  go  through  those  computations  for 
every  interconnect  in  the  layout.  Simple  formulas,  though  uot  as 
accurate  as  the  numerical  '.-mputation,  have  to  be  derived  and 
employed  in  the  extraction  programs.  From  another  point  of  view, 
simple  or  approximate  formulas  can  also  provide  some  physical  insight 
or  "guidelines"  for  the  designers.  Hence,  designers  can  take  the 
electrical  parameters  of  the  interconnects  into  account  in  the  design 
process  and  reduce  the  number  of  design  iterations.  In  general, 
these  approximation  formulas  can  serve  as  useful  tools  in  both  design 
and  design  verification. 


3_.2.  Se  1  f-Cao  ac  i  t  anc  e 


In  this  section,  the  approximation  formula*  for  the  self- 
cspscitsnce  are  discnssed.  The  self-capacitance,  by  definition,  is 
the  capacitance  of  the  interconnect  with  respect  to  the  ground  plane. 
It  should  be  pointed  ont  that  the  presence  of  other  interconnects  in 
the  vicinity  will  reduce  this  capacitance  with  respect  to  the  ground 
plane  as  discnssed  in  Sec.  2.7,  and  shown  in  Table  2.2. 

5. .2.1.  A  review  of  simple  formulas  for  self-capacitance 

The  simple  analytic  formulas  for  the  two-dimensional  self¬ 
capacitance  have  been  discnssed  in  the  literature  [20-23] .  For  the 
sake  of  comparision  and  discussion,  a  general  review  of  chose  simple 
formulas  is  first  given. 


Parallel-plate  formula 

The  parallel  plate  formula  is  a  fundamental  expression  for  the 
computation  of  capacitance.  It  can  be  written  as 

C  *  «w/h  (3.1) 


where  e  is  the  dielectric  constant,  w  is  the  width  and  h  is  the  dis¬ 
tance  from  the  ground  plane  to  the  plate.  Note  that  the  formula  m 
Eq.  (3.1)  is  two-dimensional,  and  the  capacitance  calculated  by  it  is 
in  units  of  Farads  per— unit-length.  Generally  speaking,  this  formula 
is  applicable  nnder  the  assumption  that  the  thickness  of  the  plate  is 
negligibly  small  and  the  width  w  is  considerably  larger  compared  with 
the  distance  h.  It  can  be  considered  as  a  first-order  approximation 
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formula  and  it  was  used  successfully  for  the  calculation  of  intercon¬ 
nect  and  gate  capacitance  before  the  dimensions  became  so  small  that 
fringing  capacitance  was  no  longer  negligible. 

Thin— plate  formula 

It  is  known  that  the  charges  on  a  conductor  surface  will  congre¬ 
gate  at  places  where  there  is  a  sharp  curvature  change,  e.g.,  edges 
and  points.  Hence,  the  electric  fields  at  the  edges  are  not  constant 
as  in  the  middle  of  the  plate  as  assumed  in  the  parallel-plate  for¬ 
mula.  Thus,  it  gives  rise  to  the  "edge  effect"  or  fringing-f ield 
capacitance  of  the  plate.  In  [20]  a  general  formula  was  derived  from 
the  conformal  mapping  method  to  include  the  edge  effect  of  infini¬ 
tesimally  thin-plates.  An  approximate  formula  assuming  that  width  w 
is  larger  than  the  distance  h  can  be  expressed  as 

w  2h  jtw 

C  -  e - [  1  +  - (  1  +  In - )],  w>>h  (3.2) 

h  irw  h 

The  second  term  in  the  expression  can  be  considered  to  account  for 
the  fringing-f ield  capacitance.  As  w  is  increased  larger  than  h, 
this  term  becomes  negligible  and  the  thin-plate  formula  Eq.  (3.2),  as 
expected,  reduces  to  the  parallel-plate  formula  of  Eq.  (3.1). 

i*l*i*i*  Chang's  formula 

Advancing  one  step  forward,  Chang's  formula  [21]  furthermore 
takes  into  account  the  edge  effect  and  the  finite  thickness  of  the 


plate.  Utilizing  the  Schwartz-Chr istof fel  transformation  of  the  con¬ 


formal  mapping  method  and  an  ingeniously  contrived  transformation 
function,  Chang  was  able  to  derive  an  analytic  formula  for  two- 
dimensional  self  capacitance.  It  can  be  summarized  in  the  following: 

C  -  e(2/n)ln(2Jab/Ri) ,  w/h  2  1  (3.3) 

lnR^  «  -i  -  (mr/2h)  -  (p+l)p“1/*tanh“1(p-1/1)  -  In  ((p-l)/4p) 

R^1)-[(R^1)-l)(R^1)-p)]X/1+(p+l)tanh"1[(R^1)-p)/(R^1)-l)]1/l 
-2pl/,*tanh“X[(R^1^-p)/p(R^1^-l))1^*+(nw/2h)p1^*,  5  >  w/h  >  l 

otherwise 

8^X^  *  q  +  (p+l)/2*lnA 

q  ■  p  ^*(jrw/2h)  +  (p+l)/2[l  +  ln(4/(p-l))]  -  2pX^*tanh-Xp-X^X 
A  -  max(q,p) 

p  -  2B*  -  1  +  [(2B*-1)*-1]X/X 

B  -  1  +  t/h 


The  parameters  w,h,e  are  as  defined  in  Eq.  (3.1),  and  "t 


is  the 
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thickness  of  the  piste.  It  should  be  pointed  out  that  this  formula, 
though  complicated,  is  rather  accurate.  In  comparison  with  the  two- 
dimensional  numerical  computation  [4],  its  percentage  errors  gre  all 
within  ±2%  [21]  for  different  combinations  of  w/h  and  t/h.  Still 
this  formula  is  considered  too  complicated  to  be  employed  in  the 
extraction  programs. 

i*l*I*l»  Elmasry*  s  formula 

In  line  with  the  idea  to  provide  an  easy,  simple  formula  in  the 
preliminary  design  phase  and  in  the  CAD  programs,  another  simple  for¬ 
mula  [22]  was  proposed.  It  can  be  expressed  as 

w  h  t  t  w/  2 

C  -  s - [  1  +  2 - ln(l+  - )  +  2 - ln(l+  - )]  (3.4) 

h  w  h  w  h+t 

Note  that  the  parameters  are  defined  as  before.  Physically  the 
second  and  the  third  terms  in  the  above  expression  correspond  to  the 
capacitances  resulting  from  the  charges  on  the  side  walls  and  the  top 
wall  of  the  plate,  respectively.  Since  this  is  an  empirical  formula, 
the  factor  2  in  front  of  the  second  and  the  third  terms  can  be 
adjusted  empirically  according  to  the  experimental  data  or  the  data 
from  numerical  computation  for  a  better  correlation  with  the  "true" 
capacitance  value. 

All  in  all,  the  above  formulas  can  be  applied  in  the  calculation 
of  the  self-capacitance  in  a  homogeneous  medium.  Chang's  formula  is 


accurate  enough  but  still  too  complicated  to  be  incorporated  in  the 
extraction  programs.  Elmasry's  formula  is  not  sufficiently  accurate, 
thus,  another  approximate  formula,  the  "cylindrical  approximation 
formula",  is  proposed. 

Cylindrical  approximation  formula 

1*  Formulation 

As  a  preliminary,  let  us  consider  the  capacitance  of  an  inter¬ 
connect  whose  cross-section  is  a  square  with  side  t  and  a  distance  h 
over  the  ground  plane.  It  is  known  that  the  capacitance  is  deter¬ 
mined  by  the  geometric  shape,  the  perimeter  and  the  distance  from  the 
ground  plane.  Hence,  the  capacitance  per  unit  length  of  the  square, 
Cs,  will  be  a  function  of  t  and  h.  Since  the  capacitance  per  unit 
length  of  a  conductor  over  a  ground  plane  with  a  circular  cross- 
section  is  given  by  the  relatively  simple  expression  [151 


C 

c 


2ns 


S  i. 1/1 

d+(d  -a  )  ' 
ln( - ) 


(3.5) 


where  "a"  is  the  radius  of  the  circle,  and  "d"  is  distance  from  its 
center  to  the  ground  plane,  then  intuitively,  we  might  employ 
Eq.  (3.5)  as  an  estimation  for  the  capacitance  C#  with  a  proper 
choice  of  parameters.  Next  we  examine  a  few  possible  choices  for 


these  parameters 


Consider  the  inscribed  circle  to  the  square;  that  is,  a  »  t/2, 
and  d  ■  a  +  h.  Because  like  charges  repel  on  the  conductor,  they 
will  keep  as  far  apart  as  possible.  Consequently,  saooth  surfaces 
like  the  circle  will  have  less  charge  accumulated  near  the  ground 
plane  than  the  square  which  has  sharp  corners  for  the  congregation  of 
charges.  Also,  the  "inscribed**  circle  has  a  perimeter  irt  compared  to 
a  perimeter  of  4t  for  the  square.  Due  to  this  shape  factor  plus  the 
smaller  perimeter  of  the  circle,  it  is  expected  that  the  "inscribed 
circle  approximation"  will  always  underestimate  the  capacitance  per 
unit  length  of  the  square  configuration. 

Next  let  ns  consider  the  conductor  with  a  rectangular  cross- 
section  ,  as  shown  in  Fig.  3.1  with  width  w  and  thickness  t.  Also, 
shown  in  Fig.  3.1  is  a  cylindrical  conductor  whose  ends  are  circular 
with  a  radins  t/2.  For  this  conductor  the  total  "side-wall  capaci¬ 
tance"  can  be  estimated  with  Eq.  (3.3),  where  a  ”  t/2  and 
d  ■  t/2  +  h.  In  the  center  part  of  the  conductor,  a  uniform  charge 
distribution  is  assumed  throughout  the  bottom  side  of  the  conductor, 
so  that  the  capacitance  of  this  section  can  be  computed  with  the 
parallel-plate  formula  of  Eq.  (3.1).  It  was  found  by  trial  and  error 
that  the  capacitance  of  the  cylindrical  conductor  closely  approxi¬ 
mated  the  capacitance  of  the  rectangular  conductor  when  a  section  of 
width  t/4  was  taken  from  each  end  of  the  rectangle  and  replaced  with 
a  semicircular  section  or  radius  t/2,  as  shown  in  Fig.  3.1.  The 
total  perimeter  of  the  semicircular  ends  is  jtt  while  the  total  perim- 
eter  removed  from  the  rectanngle  is  3t.  Thus,  there  is  a  close 
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Figure  3.1  Conductor  with  t  rectangular  cross-section 
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match.  With  the  above  assumptions  and  for  the  case  v  >  t/2.  the 
capacitance  per  unit  length  of  the  rectangular  conductor  can  be 
estimated  with  the  sample  formula 


w-t/2  2n 

■  s  [  +  V/  »  ^ ' 

h  ln[l+2h/t+{2h/t(2h/t+2) )  '  ] 


»  >  t/2 

(3.6? 


When  w  <  t/2.  it  is  not  possible  to  remove  sections  of  width  t/4 
from  each  end.  Although  this  case  is  uncommon  in  integrated  cir¬ 
cuits,  it  was  found  that  a  reasonable  approximation  to  the  capaci¬ 
tance  could  be  obtained  by  using  the  parallel-plate  calculation  for 
the  entire  width  and  modifying  the  circular  capacitance  approximation 
in  order  to  account  for  the  fringing  capacitance  and  to  maintain  con¬ 
tinuity  with  Eq.  (3.6).  The  resulting  formula  for  estimating  the 
total  capacitance  per  unit  length  in  this  case  is 


e[- 


n(l-0.0543t/2h) 


+  1.47], 


ln[l+2h/t+(2h/t(2h/t+2))  '  ] 


<  t/2 
(3.7) 


In  both  Eq.  (3.6)  and  Eq.  (3.7),  the  second  term  attempts  to  account 
for  the  additional  capacitance  due  to  fringing  effects  at  the  edges 
of  the  conductor.  In  the  case  of  two  layers  of  media,  the  first 
term,  parallel-plate  term,  and  the  fringing  field  term  are  weighed 
with  different  dielectric  constants.  This  will  be  discussed  later. 
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2.*Z>2.2.  Results  and  comparison 

Since  the  accuracy  of  Chang's  formula  has  been  established  in 
comparison  with  the  two-dimensional  numerical  computation,  the  capa¬ 
citance  values  generated  by  it  are  used  as  the  "true”  capacitance  *  in 
subsequent  comparisons.  First  of  all,  the  normalized  capacitances  by 
the  "cylindrical  approximation  formula"  are  compared  with  those  by 
Chang's  formula  in  Fig.  3.2.  Those  capacitances  are  normalized  with 
respect  to  the  parallel-plate  capacitance,  and  they  are  plotted  in 
Fig.  3.2  for  t/h  ■  0.1,  1,  and  10.  As  w/h  increases,  all  the  normal¬ 
ized  values  approach  constant  1.  In  other  words,  the  capacitances, 
as  expected,  approach  the  parallel-plate  capacitance  when  the  width 
is  relatively  larger  than  the  distance  from  the  ground  plane.  For 
t/h  -  10.  the  cylindrical  approximation  is  within  5%  of  the  actual 
capacitance,  and  for  t/h  *1,  the  error  can  be  as  large  as  10%,  but 
only  for  w/h  <<  1.  For  t/h  *0.1,  the  error  is  within  20%. 


Furthermore,  for  a  more  detailed  comparison  of  various  approxi¬ 
mation  formulas,  the  percentage  errors  of  the  capacitance  values  com¬ 
puted  by  the  thin-plate  formula,  Elmasry' s  formula,  and  the  cylindri¬ 
cal  formula,  are  plotted  in  Fig.  3.3,  Fig.  3.4,  and  Fig.  3.5  for 
t/h  *  0.1,  1,  and  10,  respectively.  The  referenced  "true"  capaci¬ 
tance  is  obtained  by  Chang's  formula  Eq.  (3.3).  Note  that  the  thin- 
plate  formula  and  Elmasry' s  formula  are  logarithmic-type  approxima¬ 
tions,  and  the  errors  are  relatively  larger  for  smaller  w/h.  While 
the  cylindrical  formula  is  a  reciprocal  logarithmic  approximation, 
the  error  trend  is  revei. id  for  smaller  w/h.  Hence,  the  error  for 
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Figure  3.5  Comparision  of  errors  of  the  approximate  formulas 


the  cylindrical  approximation  is  bounded  within  a  smaller  range  and 


is  more  constant  over  the  range.  For  example,  at  t/h  “  1.  the 
cylindrical  approximation  has  an  error  within  +5%  for  w/h  from  1  to 
50.  As  for  the  other  two  approximate  formulas,  the  errors  are  30% 
and  40%.  However,  for  the  extremal  case,  t/h  *  0.1,  the  thin-plate 
approximation,  Eq.  (3.2),  generally  has  less  error  than  the  cylindri¬ 
cal  approximation.  It  may  be  inferred  that,  for  this  case,  the 
charges  are  distributed  more  on  the  bottom  side  of  the  interconnect, 
so  that  the  infinitesimally  thin-plate  approximation  is  more  realis¬ 
tic  than  the  cylindrical  approximation.  Consequently,  Eq.  (3.2)  is  a 
better  approximation  is  this  case.  For  general  cases,  the  cylindri¬ 
cal  approximation  formula  compares  fairly  favorably  with  the  other 
approximate  formulas. 

3. Two  layers  of  media 

If  the  interconnect  is  situated  within  a  passivation  layer, 
<s„g.,  SiN,  on  top  of  the  oxide,  SiO,,  then  the  difference  of  the 
dielectric  constants  has  to  be  tafceu  into  account  in  the  calculation 
of  the  capacitance.  This  can  be  accomplished  by  selecting  different 
weighting  factors  for  the  "parallel-platf ”  term  and  the  "fringing" 
term  in  Eq.  (3.6)  an  Eq.  (3.7).  Consequently,  for  two  layers  of 
dielectric  media,  the  total  capacitance  is 


In  order  to  examine  the  errors,  two-dimensional  numerical  compu¬ 


tation,  e.g.,  "CAP2D" ,  is  performed  to  calculate  the  actual  capaci¬ 
tance.  In  this  calculation,  it  was  assumed  that  the  passivation 
layer  is  thick  enough  so  that  the  field  lines  passing  through  the 
outer  boundary  are  negligible.  Thus,  the  thickness  of  the  passiva¬ 
tion  layer  was  considered  to  be  infinite  in  the  calculation  of  the 
actual  capacitance.  The  percentage  error  between  this  capacitance 
and  the  capacitance  obtained  from  the  simplified  formulas  is  shown  in 
Table  3.1.  The  cylindrical  approximation  means  the  formulas  of 
Eq.  (3.6)  and  Eq.  (3.7)  are  modified  according  to  Eq.  (3.8),  and  the 
thin-plate  approximation  means  Eq.  (3.2)  is  modified  according  to 

TABLE  3.1  Percentage  Error  in  Capacitance  Estimation  with 
Two  Dielectric  Layers. 


Percent  Error 
cylindrical  approximation 


thin  plate 
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Eq.  (3.8).  Note  that  the  errors  of  the  cylindrical  approximation  at 
t/h  *  0.1  are  relatively  large  compared  vith  those  of  the  thin-plate 
approximation,  except  at  v/h  «  1 .  As  discussed  previously,  modified 
Eq.  (3.2)  is  a  better  approximation  than  the  modified  formulas  of 
Eq.  (3.6)  and  Eq.  (3.7)  in  this  case.  However,  for  the  other  cases, 
t/h  -  1  and  10,  the  percentage  errors  of  the  cylindrical  approxima¬ 
tion  are  well  within  ±10%.  Despite  the  already-existing  errors  in 
Eq.  (3.6)  and  Eq.  (3.7),  the  simple  weighting  factors  in  Eq.  (3.8) 
approximate  the  case  of  two  dielectric  layers  satisfactorily.  There¬ 
fore,  the  simple  formulas,  except  Eq.  (3.3),  may  also  be  applicable 
for  the  capacitance  of  interconnect  in  VLSI  circuits  when  the  pas¬ 
sivation  layer  is  taken  into  account. 

Cylindrical  approximation  on  the  overlap  capacitance 

The  overlap  capacitance  between  the  gate  and  the  source  or  the 
gate  and  the  drain  in  an  MOS  device  is  an  important  parameter  in  cir¬ 
cuit  simulation.  It  affects  the  MOS  model  in  the  circuit  simulation 
program  and  impacts  on  the  circuit  response.  Typically  it  was 
estimated  by  the  parallel-plate  formula  as  Eq.  (3.1).  As  a  side 
topic,  it  was  found  out  that  the  cylindrical  approximation  formula 
can  also  be  applied  to  the  estimation  of  the  overlap  capacitance. 

A  simple  model  and  formula  were  derived  for  the  overlap  capaci¬ 
tance  in  an  MOS  device  in  123}  .  Assuming  the  side  walls  of  the  gate 
and  those  of  the  source  are  vertical  boundaries,  as  shown  in 
Fig.  3.6,  the  simple  formula  derived  in  [25]  was  expressed  as 


cln<1+  “T— J *  +  “ox—; -  +  "7*81  ln(1  +  , 

'ox  tox  ”  toz 


■sm0) 

(3.9) 


'ox  . 

A  m  -  - - - 

2  sinB 


Note  that  Mt  *  is  the  gate  oxide  thickness,  "d**  is  the  overlap 
width  between  the  gate  and  the  sonrce,  "x^"  and  nx.H  are  the  thick¬ 
nesses  of  the  gate  and  the  sonrce,  respectively.  The  accuracy  of  this 
forarala  is  manifested  in  [25]  by  comparison  with  experimental  data 
and  two-dimensional  device  model  simulation. 

Employing  the  basic  idea  of  the  cylindrical  approximation  for 
the  fringing  fields  of  the  finite  thickness  of  the  gate  and  the 
source  of  an  MOS  device,  we  can  obtain  the  cylindrical  approximation 
on  the  overlap  capacitance  as 


where  0 


ln(l+0+(0(p+2>)  '  ) 

•  in(i+rKr<r+2)) 

;p  ;  r  -  4tox/Xj 

(3.10 

+  «„.>/2 

82 


r-  .*» 


The  first  tern  in  Eq.  (3.10)  accounts  for  the  parallel-plate  com¬ 
ponent  of  the  capacitance*  the  second  one  results  from  the  fringing 
field  of  the  gate*  and  the  third  term  is  due  to  the  fringing  field  of 
the  source.  Because  the  source  is  embedded  inside  the  silicon  sub¬ 
strate,  the  difference  in  the  dielectric  constant  between  silicon 
oxide  and  the  substrate  is  taken  into  account  by  the  factor  e 

avg* 

A  few  selected  data  points  comparing  Eq.  (3.10)  with  Eq.  (3.9) 
are  given  in  Table  3.2.  The  errors  are  mostly  bounded  within  1%  for 

TABLE  3.2  Comparison  of  Simple  Formulas  for  Overlap  Capacitance. 


d 

t 

X^ 

XJ 

• 

c 

c*  * 

ox 

p 

ov 

ov 

0.2 

0.035 

0.4 

0.4 

4.104 

4.074 

0.4 

0.035 

0.4 

0.4 

6.128 

6.098 

0.6 

0.035 

0.4 

0.4 

8.152 

8.122 

0.8 

0.035 

0.4 

0.4 

10.176 

10.145 

1.0 

0.035 

0.4 

0.4 

12.200 

12.169 

0.5 

0.5 

1.0 

1.0 

1.295 

1.301 

0.5 

0.5 

0.5 

0.5 

1.009 

1.082 

1.0 

0.02 

0.2 

0.5 

20.234 

20.342 

0.5 

0.02 

0.5 

0.2 

11.025 

11.134 

10.0 

1.0 

10.0 

10.0 

5.519 

5.467 

*  C  is  the  capacitance  calculated  by  Eq.  (3.9),  and  C’ 
is  that  by  Eq.  (3.10) .  ov 
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various  comb iaat ions  of  paramters.  Thus,  it  can  be  said  that  the 
cylindrical  approzimation  on  the  overlap  capacitance  is  as  good  a 
siaple  formula  as  that  proposed  in  [25] .  Hence,  the  cylindrical 
approziaation  can  be  more  generally  applied  to  the  estimation  of  the 
fringing  capacitance  due  to  the  thickness  of  a  plate. 

To  sua  up,  the  "cylindrical  approziaation  formula"  seeas  capable 
of  being  generally  applied  to  the  estimation  of  the  fringing  field 
capacitance,  be  it  due  to  the  side  walls  of  an  interconnect,  or  those 
of  a  gate  or  a  source.  Besides,  the  errors  from  it  are  bounded 
within  a  reasonable  limit,  and  the  ezpressions  are  simple  enough  for 
the  application  in  the  eztraction  programs  for  design  verification. 


!•£•  Coupling 


In  VLSI  circuits,  as  the  design  rules  become  more  stringent 
interconnects  are  cloaer  together;  inevitably,  the  coupling  capaci¬ 
tances  between  the  interconnects  become  important  and  have  to  be 
reckoned  with  in  the  circuit  simulation.  By  numerical  computation 
methods,  they  can  be  evaluated  rather  accurately  along  with  the 
self-capacitances.  However,  as  mentioned  before,  it  is  essential  to 
have  a  simple  formula  for  coupling  capacitances  in  order  to  incor¬ 
porate  them  into  the  eztraction  programs  for  design  verification. 

Since  the  coupling  capacitances  are  strongly  related  to  the 
self-capacitances,  it  is  not  easy  to  treat  them  separately  and  derive 
a  simple  formula  solely  for  the  coupling  capacitance.  For  ezample, 
the  simple  formulas  derived  in  [24]  are  for  the  estimation  of  the  sum 


of  the  coupling  end  the  self-capacitances  for  two  interconnects. 
However,  one  method  for  generating  a  simple  formnla  for  the  coupling 
capacitance  uses  a  "least  square  fit"  on  the  data  obtained  from 
numerical  computation.  In  [24]  the  simple  formulas  have  been 
obtained  by  this  method.  Employing  the  programs  described  in  Chap. 2. 
i.e.,  "CAP2D"  and  "CAP3D".  we  can  collect  a  few  sets  of  data  on  dif¬ 
ferent  configurations  of  interconnects.  Then  utilizing  a  subroutine, 
we  can  fit  different  formulas  with  multiple  coefficients  to  the  data 
in  search  of  better  simple  formulas. 

1.2.4,.  Coupling  capacitance  between  two  parallel  interconnects 

Let  us  consider  the  configuration  of  two  parallel  interconnects 
in  a  Homogeneous  medium  similar  to  the  figure  shown  in  Fig.  2.8.  In 
the  most  general  case,  there  are  seven  parameters  to  consider  the 
width  V,  the  thickness  "t" ,  and  the  distance  Mhn  from  the  ground 
plane  for  each  interconnect,  plus  the  separation  "s”.  In  order  to 
reduce  the  variables  in  the  simple  formnla,  let  us  assume  that  both 
interconnects  possess  the  same  parameters,  and  one  of  the  parameters 
is  utilized  as  the  normalizing  factor;  then,  in  essence  there  are 
three  variables.  In  [24]  a  simple  formula  for  two  interconnects  with 
identical  parameters  was  derived.  It  estimated  the  sum  of  the  self- 
and  the  coupling  capacitances  for  one  interconnect.  It  is 


c  -  «0XU.15(w/h)  +  2 . 80 ( t/h)  +  (s/h) 

[0.03( w/h)  +  0.83( t/h)  -  0.07(t/h)# **il] }  (3.11) 

The  Accuracy  of  the  formula  haa  been  shown  in  [24]  to  be  within  +10% 
of  the  two-dimensional  numerical  computation.  But  the  coupling  and 
the  self-capacitances  are  not  separated. 

In  order  to  obtain  a  simple  formula  solely  for  the  coupling 
capacitance,  the  same  "least  square  fit"  method  is  applied  on  the 
data  calculated  by  "CAP2D",  and  a  similar  formula  is  obtained  : 


Cia  *  eox[0.'004«(w/t)1*‘1  +  0.0271(h/t)X**<  +  (s/t)-1*#* ' . 

[O^dCm/t)***4*  +  0.914(  h/t)#*2#1]}  (3.12) 

The  errors  of  this  approximate  formula,  Eq.  (3.12),  are  given  in 
Table  3.3.  Hence,  Eq.  (3.12)  can  be  nsed  in  conjunction  with 
Eq.  (3.11)  to  evaluate  the  self-capacitance  and  the  coupling  capaci¬ 
tance  individually. 

Three-dimens ional  cross-over  geometry 

In  this  section  the  configuration  of  perpendicular  cross-over  of 
two  interconnects  in  Fig.  3.7  is  considered.  As  in  the  previous  sec¬ 
tion,  to  reduce  the  variables  considered,  both  interconnects  are 
assumed  to  have  the  same  width  and  thickness.  Since  the  coupling 


between  the  two  over  the  cross-over  region  is  of  interest,  the 
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TABLE  3.3  Percentage  Error  of  Eq.  (3.12). 


s/t 

w/t 

h/t 

.500 

1.000  2.000 

5.000 

1.00 

.50 

-3.7245 

-9.5482  -3.5060 

42.2329 

1.00 

1.00 

-2.2482 

-12.6091  -12.9382 

20.2085 

1.00 

2.00 

1.6377 

-12.7294  -18.5078 

2.8354 

1.00 

5.00 

15.0499 

-.5205  -6.5252 

14.9383 

•  * 

2.00 

.50 

-1.0380 

-9.9699  -8.9054 

21.9106 

2.00 

1.00 

.2764 

-12.7097  -16.4034 

6.5781 

ft 

2.00 

2.00 

3.2419 

-13.3229  -21.3423 

-6.1157 

2.00 

5.00 

14.4955 

-3.5693  -11.8742 

4.3947 

5.00 

.50 

4.4536 

-8.4302  -12.6018 

5.5028 

5.00 

1.00 

5.6243 

-10.7517  -18.1644 

-4.0452 

5.00 

2.00 

8.1620 

-11.3784  -22.0076 

-12.5019 

3 

5.00 

5.00 

17.2850 

-3.9135  -14.8168 

-3.9239 

10.00 

.50 

11.4921 

-2.1553  -6.1729 

14.9096 

10.00 

1.00 

12.1799 

-4.7770  -11.6876 

5.6213 

10.00 

2.00 

14.3737 

-5.7158  -15.7960 

-3.4468 

10.00 

5.00 

22.6563 

.6820  -10.2202 

1.2573 

:« 

lengths  of  the  interconnects  are  assumed  to  be  five  times  the  widths 
to  reduce  further  the  number  of  variables.  Using  "t"  as  the  normal¬ 
izing  factor,  we  have  a  total  of  three  variables,  w/t,  h/t,  and  s/t. 
In  this  case,  the  parameter  ”h”  designates  the  distance  of  the  lower 
interconnect  from  the  ground  plane.  The  program  "CAP3D"  was  employed 
to  evaluate  the  coupling  capacitances  in  this  configuration,  and  the 
data  are  shown  in  Table  3.4.  A  few  approximate  formulas  composed  of 
the  variables  w/t,  h/t,  and  s/t  were  tested.  The  following  formula 
seems  to  have  the  smallest  error: 


TABLE  3.4  Coupling  Capacitance  of  the  Cross-Over  Geometry 


w/t 

h/t 

0.500 

s/t 

1.000  2.000 

5.000 

1.000 

.500 

39.323 

27.732 

16.808 

6.132 

1.000 

1.000 

41.811 

29.811 

18.386 

6.986 

1.000 

2.000 

44.749 

32.364 

20.450 

8.243 

1.000 

5.000 

48.032 

35.659 

23.406 

10.385 

2.000 

.500 

78.091 

60.677 

41.846 

19.473 

2.000 

1.000 

84.008 

65.220 

45.190 

21.395 

2.000 

2.000 

89.104 

70.923 

50.281 

24.652 

2.000 

5.000 

98.349 

79.420 

57.593 

29.873 

5.000 

.500 

204.420 

168.070 

126.430 

74.454 

5.000 

1.000 

209.290 

174.900 

133.640 

79.888 

5.000 

2.000 

207.090 

179.510 

142.740 

88.719 

5.000 

5.000 

224.790 

199.590 

163  .530 

105.730 

10.000 

.500 

500.330 

412.840 

313  .010 

192.930 

10.000 

1.000 

480.290 

407.360 

318.000 

201.730 

10.000 

2.000 

436.340 

386.570 

317.970 

214.240 

10.000 

5.000 

419.600 

387  .050 

336.120 

243 .620 

xa  ”  «ox{0* 

0034(w/t) 

i  .»«• 

5 ,624( h/t) 

+ 

_0  .  »  : 

(s/t) 

[6 

.74( v/t) J 

i.xis 

+  6 

,25( h/ t) * * 

171n 

Note  that  the  dielectric  constant  eQx  in  Eq.  (3.13)  is  in  units  of 
Farads.  The  percentage  errors  of  the  approximate  formula  (3.13)  com¬ 
pared  vith  the  coupling  capacitances  in  Table  3.4  are  listed  in 


Table  3.5 
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Note  that  Eq.  (3.13)  has  a  similar  structure  to  that  of 
Eq.  (3.11).  Also  the  maximum  percentage  errors  are  +20%  and  -18%, 
and  the  nomial  errors  are  around  10%  as  seen  in  Table  3.5.  Since 
more  levels  of  interconnection  are  needed  in  increasingly  complicated 
VLSI  circuits,  Eq.  (3.12)  would  be  helpful  to  estimate  those  coupling 
capacitances  between  levels  of  interconnects. 
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The  approximation  formulae  discussed  iu  this  chapter  can  be 
▼ieved  as  a  means  to  incorporate  the  capacitance  information  computed 


in  Chapter  2  into  the  extraction  programs.  Another  method,  e.g., 
table  look-up,  can  also  be  used  instead  of  using  the  "least  square 
fit"  on  the  data  to  find  an  approximation  formula.  However,  these 
approximation  formulas  may  also  provide  a  physical  interpretation  and 
are  simpler  to  handle.  Thus,  they  are  also  useful  in  the  circuit 
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CHAPTER  4 


NETWORK  EXTRACTOR 


Tie  design  process  for  LSI  or  VLSI  circuits  still  is  not  fully 
autoaated  todsy,  but  contsins  aore  or  less  aanuel  design  steps. 
These  annuel  end  tie  interactive  layout  techniques  which  are  used  to 
optiaize  the  layout  design  inevitably  introduce  huaan  errors  and 
inaccuracy  into  the  design.  Consequently,  it  is  imperative  to  fully 
cheek  the  layout  and  verify  the  design.  Tith  the  ever-growing  cir¬ 
cuit  coaplexity.  soae  autoaated  tools  for  design  verification  are 
necessary.  Figure  4.1  shows  a  process  flow  of  design  and  the 
interaction  of  the  design  tools  at  various  stages.  The  first  step  in 
the  design  verification  is  the  artwork  (or  topological  or  aask  or 
layout)  analysis.  This  step  basically  transforas  the  layout  inforaa- 
tion  back  to  the  circuit  info rant ion  for  recheeking  the  circuit  per- 
foraance .  It  can  generally  be  classified  into  three  distinct,  but 
highly  interrelated  categories  [33]  :  ( i)  connectivity  checks,  (ii) 
design  rule  checks,  and  (iii)  the  calculation  of  electrical  paraae- 
ters.  The  connectivity  oheeks  ezaaine  the  physical  connection  of  the 
layout  eleaents  to  see  if  there  are  short  circuits  or  open  circuits. 
Design  rule  checks  basically  ezaaine  the  tolerances  between  geoaetric 
eleaents  to  see  whether  they  are  in  accordance  with  the  process 
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requirement*.  The  calculation  of  electrical  parameters  should 
include  both  the  transistors  and  the  interconnect  parameters.  For 
exasqple,  the  ehannel  width,  the  channel  length,  the  areaa  of  source 
and  drain,  and  the  overlap  capacitances,  etc.,  are  important  electri¬ 
cal  parameters  for  transistors.  The  capacitance  and  the  resistance 
are  the  electrical  parameters  for  the  interconnect.  In  the  litera¬ 
ture,  most  network  extractors  that  have  been  developed  concentrated 
on  the  connectivity  cheeks  and  the  design  rule  checks  [27,36],  or  the 
efficiency  of  the  algorithms  [32,34,33],  or  the  hierarchical  combina¬ 
tion  of  the  coaputer-aided  tools  [32,37] .  Not  too  much  effort  has 
been  spent  on  the  extraction  of  the  electrical  parameters.  Hence,  in 
the  network  extractor  developed  here,  our  main  goal  is  to  concentrate 
on  the  third  aspect  of  the  layout  analysis,  i.e.,  the  calculation  of 
electrical  parameters,  especially  the  interconnect  capacitances. 

In  the  face  of  the  growing  circuit  complexity,  the  lengths  of 
the  interconnects  on  the  chip  have  to  increase  as  pointed  out  in 
[18].  Due  to  the  smaller  sizes  of  the  transistors  and  this  increase 
of  interconnects,  the  electrical  parameters  of  the  interconnects 
become  essential.  Furthermore,  those  parameters  usually  are  not 
taken  into  account  in  the  design  phase  by  the  interactive  layout 
technique,  so  it  is  important  to  include  them  in  the  design  verifica¬ 
tion  to  ensure  satisfactory  circuit  performance.  In  the  previous 
chapter,  approximate  formulas  for  the  interconnect  capacitances  have 
been  discussed.  They  are  used  in  the  extractor  discussed  in  this 
chapter  to  calculate  the  interconnect  capacitances. 


4,.£.  General  Pronertiee  of  the  Extractor 

The  extractor*  ia  coded  in  PASCAL,  and  some  important  propertiea 
of  it  can  be  summarized  as  follows  : 

(1)  CIF  format  :  The  input  of  the  layont  information  to  the  extrac- 
tor  should  be  in  CIF  (Caltech  Intermediate  Form)  format  aa 
defined  in  [38].  However,  not  all  the  definitions  of  CIF  are 
recognized,  e.g.,  round  flash  and  some  user— defined  extensions 
are  not  known  to  the  extractor  and  errors  will  result. 
Nevertheless,  there  is  an  option  **-u"  in  running  the  package  to 
ignore  the  errors  resulting  from  user-defined  extensions.  It 
should  also  be  pointed  out  that  the  coordinates  of  the  input 
rectangles  should  be  integers;  in  other  words,  the  fraction 
numbers  after  the  decimal  points  in  the  coordinates  will  be 
chopped  off. 

(2)  NMOS  technology  :  The  extractor  assumes  NMOS  technology  and  the 
design  rules  defined  in  [38]. 

(3)  Manhattan  geometry  :  Since  most  layouts  do  not  contain  angular 
geometries  and  for  the  sake  of  simplicity  in  the  extraction  pro¬ 
cess,  all  the  geometries  on  the  layout  should  be  rectilinear 
with  respect  to  either  the  x  or  y  axis.  However,  the  original 
version  of  the  first  part  of  the  code  from  Mike  Graf  was 
intended  to  include  also  the  non-Manhattan  geometry.  Conse- 

*The  first  part  of  the  code  which  reads  in  CIF  file  and  sets  up 
the  4-d  binary  trees  for  the  rectangles  on  every  mask  is  due  to  Mike 
Graf  [40]. 
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quently,  there  is  soae  redundancy  in  the  definition  of  the  rec¬ 
tangle  record. 

(4)  Rectangle-based  data  structure  :  The  baaic  entity  of  the  data 
atruetnre  is  a  rectangle  on  a  aaak.  All  the  rectangles  on  the 
saae  mask  are  then  linked  together  to  fora  a  4-d  binary  search 
tree  with  right  thread.  [26.40] 

(5)  Transistor  configuration  :  It  is  assuaed  that  the  overlap  of 
diffusion  and  polysilicon  regions  constitutes  a  translator  chan¬ 
nel.  For  the  sake  of  sia^licity  the  extractor  only  recognizes 
basically  two  configurations  of  transistors,  i.e..  those  foraed 
by  one  diffusion  rectangle  and  one  polysilicon  rectangle  or 
those  by  tvo  diffusion  rectangles  and  one  polysilicon  rectangle. 
Specialized  transistor  configurations,  e.g.,  those  long,  serpen¬ 
tine  channel  regions  of  output  buffers,  are  not  recognizable. 
Warning  Messages  will  be  generated  for  these  unrecognizable 
overlapped  regions  of  diffusion  and  polysilicon  in  the  output 
log  file.  Functionally,  the  extractor  will  categorize  all  the 
transistors  into  three  different  types  of  transistors,  naaely, 
load,  drive  and  pass  transistors.  The  channel  lengths  and  the 
widths  of  the  transistors  are  evaluated  and  included  in  the  out¬ 
put  SPICE  file. 

(6)  Output  foraat  :  The  transistors,  resistors  and  capacitors  are 
reported  in  a  fora  eoapatible  with  the  SPICE  input  foraat.  Also 
the  node  nuabering  is  autcaatically  done  in  the  extractor.  For 
exaaple,  the  node  "0"  is  reserved  for  the  ground  node,  and  ”5" 
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is  for  the  power  node.  Every  center  node  of  the  T  equivalent 
leaped  eirenit  for  an  interconnect  has  a  node  nnaber  5xx  in 
order  to  distinguish  it  froa  other  node  nnabers  which  start  froa 
"6*  and  are  increaented  one  by  one. 

(7)  Paraaeter  file  :  All  the  process  paraaeters  are  stored  in  the 
"const ant. h"  file.  For  exaaple,  the  field  oxide  thickness,  the 
thickness  of  the  interconnects,  sheet  resistance  of  polysilicon, 
aetal  and  diffusion,  etc.,  all  should  be  given  in  this  file,  and 
the  extractor  will  fetch  the  information  froa  this  file  in  the 
process  of  extraction.  In  order  to  screen  off  small  resistances 
and  capaoitanoes  froa  short  interconnects,  two  paraaeters 
"CLIXir*  and  "ELIMIT"  can  be  set  in  this  file  snch  that  the 
interconnect  resistances  and  capaoitanoes  smaller  than  those 
values  will  not  be  reported.  Other  paraaeters  in  "constant .h" 
include  "VIXEGAP"  and  "BXANCHGAP",  which  are  used  in  the  compu¬ 
tation  of  coupling  capacitances,  and  will  be  discussed  in  a 
later  section. 

1-1-  P»ta  ataaiasi 

In  this  section,  the  data  structure  nsed  in  the  program  is  dis¬ 
cussed.  Since  the  prograa  is  coded  in  PASCAL,  the  "record"  defini¬ 
tions  of  "rect",  "transistor",  "cntlist",  "wire",  etc.,  are  instruc¬ 
tive  to  illustrate  the  data  structure. 

The  record  of  a  rectangle  in  the  program  is  defined  as  : 

e 

rect 


record 


key  :  array[0. .3]  of  real;  (•  four  keys  : 

key[0]  “  rain 
key[l]  ■  yaia 
key [2]  ■  -xaax 
key[3J  ■  -yin ax  •) 

seek  :  aasklevel;  (•  -7... 7*) 

Ison,  (•  left  son  in  tree  structure  •) 

xjou  :  rectptr;  (•  right  son  in  tree  structure  *) 

ndisc  :  integer;  (•  discr ininator  of  successor - 

only  valid  if  rson  is  a -thread  •  ) 
lapdieator  :  integer;  (*  an  indicator  •) 
end; 

rectptr  «^rect; 


The  paraaMter  "mask*  in  the  above  definition  indicates  the  aask 
level  on  which  this  rectangle  resides.  The  mask  levels  are  desig¬ 
nated  as  follows  : 

1  -  diffusion 

2  ■  ion  iaqplant 

3  “  polysilicon 

4  *  contact  out 

3  -  aMtal 

6  ■  buried  contact 

7  ■  overgl seeing 

Usually  we  are  concerned  with  only  four  mask  levels,  i.e.,  diffusion, 
polysilicon,  contact  cut  and  aetsl.  The  levels  of  buried  contact 
cut  and  overglassing  are  neglected  in  this  extractor.  Note  also 
that,  as  a  rule,  the  pointer  to  a  certain  data  type  is  naaed  by 
appending  "ptr"  to  the  naae  of  that  data  type.  Thus,  "rectptr"  is 
the  pointer  to  rect,  and  "rectptr  -/vrect"  is  PASCAL'S  way  of  specif¬ 
ication. 


•Only  those  specifications  which  are  utilized  in  this  extractor 
listed  here.  Others  concerning  non-Manhattan  geoaetries  and  the 
optiaisation  of  the  binary  trees  are  not  included,  though  they  are  in 
the  program. 


A  4-d  biliary  search  tree  with  right  thread  is  built  for  the  rec¬ 
tangles  on  each  stash  level.  Those  three  paraaeters  "Ison",  "rson" 
and  " ndisc "  are  used  to  link  the  rectangles  in  a  binary  search  tree. 
Each  rect  record  is  stored  as  a  node  in  the  tree,  and  Ison,  rson  are 
its  two  pointers  which  are  either  null  or  point  to  another  node  in 
the  tree.  (Or  each  pointer  can  be  considered  as  specifying  a  subtree 
to  this  node.)  The  discriminator  of  a  node  is  an  integer  between  0 
and  3  in  a  4-d  binary  tree.  Its  purpose  is  to  specify  the  relative 
orientation  of  the  node  to  its  sons  and  its  parents.  It  is  deter¬ 
mined  as  follows  :  The  root  of  the  tree  has  discriminator  0.  Its 
two  sons  have  discriminator  1.  and  so  on.  The  discriminator  of  level 
4  becomes  0  again.  The  cycle  repeats  between  0  and  3.  Generally 
speaking,  for  any  node  P  in  the  4-d  binary  tree,  let  j  be  its 
discriminator,  then  for  any  node  Q  in  the  left  branch  of  P,  it  is 
true  that  key^CQ)  <  keyj(P);  likewise,  for  any  node  R  in  the  right 
branch  of  P,  key^(R)  >  keyj(P).  In  case  keyj(Q)  ■  keyj(P),  then  the 
next  key,  (j+l)aod4,  is  compared,  until  the  order  is  decided.  If  all 
keys  are  equal,  then  the  two  recta  are  identical,  no  node  will  be 
inserted  into  the  tree.  In  aiding  the  search  process,  a  right  thread 
is  added  to  the  tree  [39].  A  negative  mask  value  in  the  rect  record 
indicates  that  the  right  son  of  this  node  is  a  right  thread  not  a 
true  son,  and  "ndisc”  contains  the  discriminator  of  this  successor 
which  is  linked  by  this  thread.  It  should  be  noted  that  the  parame¬ 
ter  "ndisc"  is  only  valid  when  the  right  son  of  the  node  is  a  right 
thread,  i.e.,  the  mask  is  a  negative  integer. 


In  Fif.  4.2  4a  example  of  soae  rectangles  on  a  mask  and  its 
corresponding  4-d  binary  search  trees  are  depicted.  The  dotted  lines 
in  the  tree  represent  the  right  threads  and  the  solid  lines  are 
"true*  links  of  the  trees.  Those  right  threads  should  always  point 
to  prey ions  nodes,  i.e.,  nodes  that  have  been  visited  before  the 
current  node  in  the  search  process.  A  "nil"  at  the  rightmost  end  of 
the  tree  signals  the  end  of  the  search.  Note  that  the  binary  search 
trees  are  different  depending  on  the  insertion  order  of  the  rectan¬ 
gles.  To  be  more  specific,  two  different  insertion  orders  and  their 
resulting  binary  search  trees  are  shown  in  Fig.  4.2.  These  two  4-d 
binary  search  trees  all  represent  the  rectangles  shown  in  Fig.  4.2, 
yet  they  have  different  structures.  Since  the  efficiency  of  the 
search  process  obviously  depends  on  the  structure  of  the  binary  tree, 
the  insertion  order  which  is  random  in  this  extractor  affects  the 
efficiency  of  the  search  process.  One  attempt  to  optimize  the  search 
process  is  to  try  to  balance  the  heights  of  those  binary  trees.  This 
was  included  in  the  original  version  by  Mike  Graf,  but  is  not  con¬ 
sidered  here.  Finally,  the  parameter  "lapdicator"  is  an  indicator 
for  the  subsequent  search  process  to  indicate  that  this  rectangle  has 
been  used  already. 

The  transistor  records  are  defined  as  : 
transistor  •  record 

ttype  :  nmostype;  (•  three  types  : 

load,  drive,  pass  *) 

width, 

length  :  real; 

drain, 

source  :  remrectptr;  (•  source  and  drain  rects  *) 
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gate  :  pdlaptr; 
dnode , 

•node , 

gnode  :  integer;  (•  node  numbers  *) 

end; 

tranptr  -/^transistor; 

Moat  of  the  parameters  defined  in  the  above  record  are  self- 
explanatory.  Those  diffusion  rectangles  residual  of  the  overlap 
regions  of  the  transistor  channels  are  designated  as  MremrectsM.  All 
of  these  remnant  rectangles  are  built  into  a  remnant  tree  as  before. 
"Gate'*  is  a  "pdlaptr"  which  is  a  pointer  which  points  to  a  record 
"pdlap"  where  the  keys  and  direction  of  the  overlap  region  are  kept. 
All  of  the  overlap  regions  of  diffusion  and  polysilicon  represented 
by  this  pdlap  are  also  built  into  a  4-d  binary  search  tree  with  right 
thread.  "Dnode",  "snode",  and  "gnode"  are  integer  variables  to  keep 
track  of  the  node  numbers  of  the  drain,  source  and  gate  of  the 
transistor.  They  are  initialized  to  negative  integers  first.  After 
extraction  is  complete  and  all  the  elements  are  interconnected,  these 
variables  should  become  positive  and  indicate  the  node  numbers  of  the 
transistors.  If  the  extraction  process  fails  to  connect  this 
transistor,  the  node  number  will  stay  negative,  and  this  is  very 
helpful  for  debugging. 

Another  useful  data  structure  in  the  program  is  a  linear  list 
which  connects  some  relevant  rectptrs,  e.g.,  the  constituting  rectan¬ 
gles  in  an  interconnect  can  be  represented  as  a  linked  list,  and  it 
is  defined  as  : 


cntlist  -  record 


this 

:  rectptr; 

next 

:  cntlistptr; 

(*  successor  in  linear  list  •) 

multipty 

:  integer; 

(•  multiplicity  •) 

cntdir 

:  direction; 

(*  direction  of  this  rect - 

used  in  coupling  cap  •) 

width 

:  real; 

(•  width  perpendicular  to 
cntdir  •) 

found 

end; 

cntlistptr  -/scntlist; 

:  Boolean; 

The  parameter  "multipty"  is  used  in  the  computation  of  interconnect 
parameters  indicating  the  multiplicity  of  this  cntlist  in  the  search 
process.  "Cntdir"  and  "width*  are  designed  for  nse  in  the  computa¬ 
tion  of  coupling  capacitances.  "Found"  is  a  toggle  switch  used  in 
the  program  to  indicate  the  status  of  this  rectangle.  It  should  be 
mentioned  that  the  creation  of  this  data  record  to  form  linked  lists 
of  rectangles  instead  of  tagging  more  fields  to  the  rect  record  aims 
at  minimizing  modifications  to  the  already  defined  rect  record  in  the 
original  version.  This  way  there  may  result  in  some  cumbersome  algo¬ 
rithms  and  redundant  procedures.  Therefore,  future  improvements  on 
the  program  may  abolish  this  data  record  and  carefully  replan 
another. 

In  order  to  identify  all  the  connected  rectangles  linked  between 
transistor  "ports"  which  constitute  an  interconnection  system,  a  data 
record  wire  is  set  up  to  represent  it.  A  transistor  port  is  defined 
to  be  the  gate,  source,  or  drain  of  a  transistor.  An  example  of  such 
a  wire  is  shown  in  Fig.  4.3  with  the  corresponding  equivalent  lumped 


circuit.  The  record  of  wire  is  defined  as  : 


wire  ■  record 

bound  :  irriy[0..3l  of  reel;  (*  rectangular  bounds  *) 
endport  :  ar r ay [0. .MXPORT]  of  portptr; 

portnua  :  integer;  (*  number  of  portptrs  in  endport  *) 

knot ray  :  array[0. .MXPORT]  of  portptr; 

knotcnt  :  array[0. .MXPORT]  of  cntlistptr; 

knot nun  ;  integer;  (*  no.  of  eleaents  in  knotray, 

knotcnt  •) 

branehray  :  array[0. .MXPORT]  of  branchptr; 
brancbnua  :  integer;  (•  no.  of  eleaents  in  branehray  *) 
wir ahead  :  cntlistptr;  (•  head  of  wire  list  •) 
end; 

wireptr  «^wire; 


"Bound"  is  the  outaost  four— corner  boundary  points  of  this  wire,  it 
is  defined  the  sane  way  as  the  keys  in  reet.  Its  aain  purpose  is  to 
screen  off  unnecessary  coaparison  between  wires  in  the  calculation  of 
coupling  capacitances.  When  two  wire  bounds  are  separated  by  aore 
than  the  parameter  "1IREGAP"  specified  in  the  "constant .h"  file,  no 
further  coaparisons  between  branches  inside  the  wires  are  needed. 

Since  a  wire  links  a  few  transistor  ports  together  at  the  ends 
of  the  wire,  a  data  structure  "port"  is  built  to  represent  those 
transistor  ports.  "Endport"  is  an  array  which  stores  those  pointers 
of  the  transistor  ports.  A  "knot"  is  a  rectangle  which  is  connected 
to  aore  than  two  other  rectangles  in  the  wire  list  as  shown  in 
Fig.  4.3.  In  order  to  eonfora  with  the  transistor  ports  at  the  boun¬ 
dary  of  the  wire,  knots  are  also  defined  as  ports.  "Enotray"  col¬ 
lects  those  port  pointers  pointed  to  the  knots  in  the  wire. 
"Knotcnt"  is  also  an  array  of  cntlistptrs  which  point  to  those  knot 
rectangles.  A  "branch"  is  the  linked  list  of  cntlistptrs  which  con¬ 
stitute  the  interconnection  between  two  ports  in  the  wire.  A  T- 


equivalent  luaped  circuit  is  specified  for  each  branch  in  the  wire. 


and  the  equivalent  resistance  and  capacitance  are  oonputed  according 
to  the  branch.  As  the  others,  "branchray"  collects  those  branches. 
For  each  array,  there  has  to  be  an  integer  *to  indicate  the  nuaber  of 


eleaents  in  the  array,  and  "portnua",  “knot nun"  and  "branchnun”  in 
the  definition  are  those  integers. 

±.±.  Extraction  Procedure 

Since  the  first  part  of  the  prograa  vhich  reads  in  CIF  input 
lines  and  sets  up  the  4-d  binary  search  trees  for  the  rectangles  on 
each  mask  level  vas  obtained  fraa  Mike  Graf  as  aentioned  earlier,  we 
would  like  to  concentrate  on  the  discussion  of  the  second  part  of  the 
prograa,  i.e.,  the  extraction  of  the  transistors  and  the  interconnect 
paraaeters.  The  prograa  is  coded  in  PASCAL  and  separated  into  nine 
subprograas  and  one  aain  prograa.  The  first  part  of  the  code  is  con¬ 
tained  in  subprograas  "exfirst.p"  and  "exsecond.p".  Note  that  "*.p" 
is  the  defaulted  format  on  the  VAX/UNIX*  system  to  indicate  a  PASCAL 
prograa. 


£*4.*i*  Find  the  diffusion  and  oolvsilicon  overlap 

First  the  transistor  channel  region  should  be  identified. 
According  to  the  NMOS  process  in  [38],  the  overlap  of  diffusion  and 
polysilicon  rectangles  constitutes  a  transistor.  Therefore,  a  search 
for  every  diffusion  rectangle  on  the  aask  through  the  polysilicon 
binary  tree  to  find  any  polysilicon  rectangle  which  overlaps  this 

*  VAX  is  a  trademark  of  Digital  Equipment  Corporation,  and  UNIX  is 
a  trademark  of  Bell  Laboratories. 
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diffusion  rectangle  should  be  done.  If  an  overlap  is  found,  a  record 
of  type  pdlap  is  set  up  in  nenory  and  a  pointer  to  it  is  generated. 
Next  all  those  pdlaptrs  are  stored  in  a  4-d  binary  search  tree  as  are 
those  of  rectangles.  The  insertion  algorithm  will  be  given  in  detail 
in  See.  4.6  and  it  is  similar  to  that  in  [26].  In  the  record  pdlap, 
the  coordinates  of  the  overlap  region,  the  pointers  to  those  two  rec¬ 
tangles  which  fora  this  overlap  region  are  kept. 

Note  that  in  the  process  of  finding  the  overlap  region,  all  the 
relative  positions  of  two  rectangles  on  a  plane  have  to  be  con¬ 
sidered.  They  are  :  overlap,  aeet  at  a  line,  meet  at  a  point,  and 
disconnection.  If  one  rectangle  is  totally  enclosed  inside  another, 
this  can  be  considered  as  a  special  case  of  overlap.  The  cases  in 
which  a  polysilioon  rectangle  and  a  diffusion  rectangle  meet  at  a 
line  or  at  a  point  are  illegal  for  the  design  rules  employed  here. 
Hence,  warning  messages  will  be  produced  in  the  log  file.  Last,  the 
above  xunctions  are  performed  by  the  subprogram  "exthird.p" . 

l.l.l.  se£.  amium 

After  the  overlap  regions  are  identified,  transistor  records 
should  be  established.  Since  the  chsnnel  region  of  a  transistor  can 
be  made  of  connecting  overlap  regions,  they  should  first  be  located. 
For  each  node  in  the  tree  of  overlap  pointers,  s  search  for  adjoining 
overlap  regions  in  this  tree  is  conducted.  All  the  found  connecting 
overlap  regiona  are  built  into  the  gate  of  a  transistor.  If  no  other 
adjoining  overlap  regions  are  found,  then  this  solitary  overlap 
region  also  constitutes  the  gate  of  a  transistor.  As  mentioned 
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before,  the  possible  configurations  of  s  transistor  are  limited  in 
this  extractor  for  the  sake  of  simplicity,  and  there  can  be  at  most 
tiro  connecting  overlap.  regions  to  form  a  tranaistor  channel. 

Next,  the  drain  and  the  source  regions  should  be  identified. 
They  are  those  remnant  regions  of  diffusion  rectangles  outside  of  the 
channel  regions  already  identified  as  gates.  The  above-mentioned 
limitation  is  mostly  intended  to  simplify  this  process.  Because  the 
reamant  regions  from  a  serpentine-shaped  channel  region  would  not  be 
rectangles,  but  rather  polygons,  they  are  not  easy  to  handle  within 
the  frame  of  the  present  data  structure.  Tith  the  above  limitations, 
the  remnants  are  still  rectangles,  and  they  can  be  separated  from  the 
channel  regions  and  identified  by  a  procedure  called  "split"  in  sub¬ 
program  "exfourth. p".  Then  respect  records  are  set  up  for  those 
remnant  regions,  and  all  of  them  are  also  linked  into  a  4-d  binary 
tree.  In  transistor  record,  the  pointers  to  the  source  and  the  drain 
"remrects"  are  also  kept  for  future  reference  in  the  extractor. 

Additionally,  the  type  of  the  transistor  should  be  defined. 
Since  the  "load"  transistor  is  assumed  to  be  depletion  type  and  has 
its  gate  connected  to  its  source,  the  special  characteristic  of  a 
contact  cut  overlapped  with  the  gate  region  is  used  as  a  criterion  to 
identify  it.  The  intersection  of  the  ion  implant  region  with  the 
gate  region  can  also  be  used  as  the  identifying  mechanism.  However, 
the  former  method  is  ea^loyed  here.  The  pointers  of  those  load 
tranaistors  are  kept  in  an  array  "loadtran"  with  an  integer  variable 
"loadnum"  to  indicate  its  number  of  elements.  The  other  transistors. 
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be  they  pass  transistors  or  drive  transistors,  are  kept  in  another 
array  "othertran" .  These  procedures  are  performed  in  subprogram 
"exfourth. p" . 

±.£•1.  Find  the  power  and  the  around 

It  is  assumed  that  all  the  ground  nodes  are  connected  by  a  net¬ 
work  of  rectangles  and  that  they  collapse  into  a  single  node  "0*  in 
the  circuit  simulation.  This  applies  to  the  power  supply  too,  and 
the  power  network  is  designated  as  node  "5".  In  order  to  string  both 
networks  in  two  linear  linked  lists,  coordinates  of  two  points  have 
to  be  fed  in  interactively.  Then  the  two  rectangles  which  encompass 
these  two  points  are  selected  as  the  "heads"  of  the  respective  lists. 
A  "depth-first"  search  [31]  is  then  performed  to  find  all  the  rectan¬ 
gles  in  the  lists.  The  procedures  "connect"  and  "setconnect"  perform 
this  function.  These  procedures  are  used  constantly  in  the  latter 
part  of  the  program  to  find  connecting  rectangles  for  "wires". 

The  procedure  "setconnect"  is  used  to  search  through  the  speci¬ 
fied  mask  level  given  the  pointer  of  a  cntlist  and  its  rectangle. 
The  data  record  cntlist,  which  was  discussed  in  the  previous  section, 
is  employed  and  its  pointer  is  named  "cntlistptr" .  Note  also  that 
the  parameter  lapdicator  in  the  reot  record  is  set  for  the  rectangle 
which  has  been  selected  so  that  it  will  not  be  used  again  in  the 
search  process.  The  "depth  first"  search  starts  with  the  head 
cntlistptr  and  its  rectptr  as  the  seed,  finds  all  the  rectangles  on 
the  given  mask  level  which  are  connected  to  if,  and  links  them  in  the 
list.  It  then  goes  to  next  cntlistptr  on  the  list,  uses  it  as  the 


toad*  and  continues  to  search  the  saae  tree  for  connecting  rectangles 
without  a  set  lapdictor  until  no  rectangles  can  be  found. 


Since  the  constituting  rectangles  of  an  interconnect  or  the 
power  network  may  not  necessarily  be  all  contained  on  one  mask  level, 
the  search  process  should  be  able  to  extend  to  other  aask  levels  as 
well.  Utilising  the  fact  that  those  interlevel  connections  can  only 
be  node  through  contact  cuts,  the  procedure  "connect"  then  performs 
this  extended  search.  It  will  call  upon  "setconnect"  and  search  not 
only  the  aask  level  of  the  seed  but  also  the  aask  level  of  the  con¬ 
tact  cuts.  The  pointers  of  the  found  contact  cuts  are  also  linked  in 
the  linear  list  in  the  "depth  first"  search.  Then,  whenever  a  con¬ 
tact  cut  is  encountered  as  a  seed  in  the  process,  all  three  levels 
(diffusion,  polysilicon,  and  aetal)  have  to  be  searched  by  the  pro¬ 
cedure  "setconnect".  This  way  all  the  connecting  rectangles,  even 
across  different  aask  levels,  can  be  found,  and  the  power  network  and 
the  ground  network  are  stored  in  two  lists  headed  by  "powerlist"  and 
"groundlist" .  The  above  procedures  sre  contained  in  the  subprogram 
"exf ifth.p". 

4.. 4,-4..  Set  un  power  and  sround  nodes 

First  the  rearects  of  "loadtrans"  sre  searched  to  see  if  they 
are  connected  with  any  rectangles  in  the  powerlist.  If  a  contact  is 
found,  then  this  rearect  is  designated  as  the  drain  of  the  transistor 
and  a  default  node  nuaber  5  is  set.  The  other  resurect  is  automati¬ 
cally  set  as  the  source  of  this  load  transistor.  Next  the  rearects 
of  "  othertrans"  are  searched  and  compared  with  the  rectangles  in  the 
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groundlist.  Any  remrect  connected  to  this  list  is  set  as  the  source 
of  the  transistor  and  a  default  node  nuaber  0  is  designated.  Another 
reareet  of  this  transistor  will  be  set  as ‘the  drain  of  the  transis¬ 
tor,  meanwhile,  the  transistor  will  be  naaed  as  "drive"  transistor. 
The  procedures  "se tpowernode"  and  "setdrivenode"  in  snbprograa 
"exsixth. p"  per fora  the  aboye  functions. 

Note  that  the  resistance  and  capacitance  of  the  power  and  the 
ground  networks  are  not  considered  in  this  extractor.  The  power  and 
the  ground  are  assuaed  to  be  at  constant  yoltages  and  are  designated 
to  be  node  5  and  node  0,  respectively.  Also  the  unset  nodes  will 
retain  their  initialized  negative  values  after  this  step. 

!•!•£•  Ei SA  »lres 

After  extracting  the  power  and  the  ground  network,  all  the 
reaaining  undecided  sources,  drains  and  gates  of  transistors  should 
be  linked  together  by  "wires".  Starting  with  the  unset  source  or 
drain  nodes  of  load  transistors,  and  then  the  unset  source,  drain  or 
gate  nodes  of  other  transistors,  we  proceed  to  find  all  the  wires. 

First,  using  an  unset  port  of  a  transistor  as  the  head  of  a 
list,  we  can  search  for  connecting  rectangles  by  the  procedure  "con¬ 
nect",  and  the  data  record  wire  is  built  in  the  process.  The  origi¬ 
nal  starting  node  of  the  transistor  is  the  first  element  in  the  array 
"endport"  in  the  wire  record.  In  the  process  of  depth  first  search, 
each  rectangle  in  the  list  is  also  coshered  with  other  unset  nodes  of 
other  transistors  in  arrays  "loadtran"  and  "othertran".  If  a  contact 
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is  found,  then  those  nodes  will  be  included  ss  ports  of  this  wire  in 
the  array  "endport".  Note  thst  this  process  only  searches  for  those 
rectangles  which  have  not  been  selected  in  the  previous  search  pro¬ 
cess.  (Lapdicator  in  the  reot  record  will  indicate  to  that  effect.) 
Thus,  the  search  process  will  end  when  all  the  connected  rectangles 
in  an  interconnection  are  found,  and  the  ending  ports  of  the  wire 
will  also  be  found  and  kept  in  the  array.  The  above  functions  are 
performed  by  the  procedure  "findwire"  in  the  subprogram 
"exeighth2.p". 


±.±.£.  Compute  interconnect  resistance  and  self  capacitance 
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From  the  data  record  wire,  the  equivalent  circuit  parameters  can 
be  computed,  i.e.,  resistance  and  capacitance.  As  shown  in  Fig.  4.2, 
the  interconnection  of  a  wire  may  be  multiple-connected,  i.e.,  there 
are  knots  at  the  intersection  of  the  branches.  In  order  to  identify 
those  knots,  the  linked  list  of  rectangles  has  to  be  searched  first. 
Any  rectangle  in  the  list  having  more  than  two  connections  to  other 
rectangles  is  considered  a  knot  and  they  are  stored  in  the  arrays 
"knotray"  and  "knotent".  A  branch  is  defined  as  a  linked  list  of 
rectangles  between  two  ports  remembering  that  knots  are  also  ports. 
For  each  branch  in  the  wire,  a  corresponding  T  equivalent  lwaped  cir¬ 
cuit  is  set  up  as  shown  in  Fig.  4.2.  Heuristically  the  two  resis¬ 
tances  are  assumed  equal  and  half  the  value  of  the  dc  resistance 
which  is  calculated  from  the  linked  list.  Furthermore,  the  capaci¬ 
tance  to  ground  is  the  self-capacitance  calculated  by  the  "cylindri¬ 
cal  approximation  formula"  given  in  Chapter  3.  Note  that  the  formula 
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is  «  two-dimensional  one,  and  the  calculated  value  from  the  formula 
has  to  be  multiplied  by  the  length  of  the  branch  for  the  total  capa¬ 
citance  value.  No  special  considerations  are  given  to  corners  and 
bends  in  the  branches. 

It  is  important  to  mention  that  the  rectangles  of  the  knots  are 
included  in  all  parameter  computations  of  the  branches  which  are 
connnected  to  them.  If  a  long  rectangle  is  designated  as  a  knot, 
then  the  calculated  parameters  of  the  branches  may  be  too  large.  In 
other  words,  this  is  a  pessimistic  estimation  of  the  parameters.  It 
is  hoped  that  this  situation  may  be  alleviated  by  careful  layout  or 
by  developing  further  procedures  to  recognize  it  and  split  the  rec¬ 
tangles  to  make  it  congruous  to  the  calculation.  Some  of  the  pro¬ 
cedures  in  the  subprograms  ”exseventh2.p”  and  wexeighth2.p"  perform 
the  above  tasks. 

Set  up  node  numbers 

Before  reporting  the  calculated  electrical  parameters,  we  have 
to  first  check  them  with  the  criteria  in  file  "constant.h*.  There 
are  a  few  cases  to  consider  :  (i)  If  both  the  capacitance  and  the 
resistance  values  are  larger  than  the  criteria,  ’'CLIMIT'  and 
"BLUUT" ,  all  components  of  the  T  equivalent  circuit  are  reported, 
(ii)  If  only  the  resistance  is  smaller  than  MRLIMITW,  then  the  node 
numbers  at  the  endports  of  this  branch  should  be  the  same,  and  the 
capacitance  should  be  inserted  between  this  node  and  the  ground. 
Furthermore,  if  the  capacitance  is  also  smaller  than  "CXIMIT" ,  just 
assign  the  two  connecting  endports  with  the  same  node  number  and 
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forget  the  T  equivalent  circuit,  (iii)  If  only  the  capacitance  is 
saaller  than  "CLIMIT" ,  no  center  node  is  needed  in  T  circuit,  and  the 
two  resistors  are  collapsed  into  one  resistor  connecting  the  two 
endports. 

Since  knots  are  considered  as  ports,  node  numbers  are  assigned 
to  knots  also.  Generally,  in  determining  the  node  numbers  at  the  end 
of  resistors,  first  check  the  connecting  ports.  If  node  numbers  have 
been  assigned  to  the  connecting  endports,  then  the  resistor  will  use 
those  node  numbers  at  its  ends.  Otherwise,  new  node  numbers  are  gen¬ 
erated.  The  center  nodes  where  capacitors  are  connected  to  ground 
are  numbered  as  5xx,  so  that  they  can  be  easily  recognized. 

£.4.4.  Compute  coupling  capacitances 

In  the  computation  mentioned  above,  only  the  self-capacitances 
of  the  interconnects  are  considered.  To  further  compute  the  coupling 
capacitances,  the  proximity  of  the  "wires"  has  to  be  checked.  First, 
the  bounds  of  the  wires  are  compared  pairwise.  If  they  are  discon¬ 
nected  by  a  margin  larger  than  "VIREGAP"  specified  in  the  file 
"constant. h",  then  no  further  comparison  will  be  made  to  avoid  possi¬ 
ble  confusions.  Otherwise,  the  branches  in  both  wires  have  to  be 
compared  pairwise.  In  this  check,  another  parameter  "BRANCH GAP*  in 
the  "constant. h"  file  is  used  to  test. the  separation  of  the  branches. 
Only  when  the  separation  between  two  interconnects  is  smaller  than 
this  parameter  is  the  coupling  capacitance  between  them  computed  by 
the  approximate  formula  in  Chapter  3.  Since  there  are  quite  a  few 
different  scenarios  of  two  close  interconnects,  e.g.,  they  can  be 


114 


parallel  on  the  same  level  or  on  different  levels,  or  a  metal  line 
can  cross  over  a  polysilicon  line  perpendicularly,  etc.,  every  case 
of  two  close  branches  has  to  be  considered  individually.  Also  these 
scenarios  may  be  dependent  upon  the  techno 'ogy  employed  and  the 
parameters  of  the  fabrication  processes.  In  the  present  version  of 
the  extractor,  silicon  NMOS  technology  is  assumed  and  only  parallel 
lines  on  the  same  level  and  the  perpendicular  cross-over  are  con¬ 
sidered  in  the  coupling  calculation.  In  Fig.  4.4,  an  example  of  two 
parallel  interconnects  is  shown.  Note  that  the  capacitance  computed 
by  Eq.  (3.12)  is  capacitance  per-unit-length;  thus,  it  must  be  mul¬ 
tiplied  by  the  total  length  in  question.  In  the  case  shown  in 
Fig.  4.4,  the  length  used  in  the  computation  is  lx+la.  This 
apparently  is  less  than  the  actual  length  because  of  the  omission  of 
the  turn.  The  T  equivalent  circuits  for  these  two  interconnects  are 
also  shown  in  Fig.  4.4.  One  of  the  advantages  of  using  the  T  circuit 
is  that  the  coupling  capacitances  can  be  easily  inserted  between  the 
center  nodes  of  the  two  T  circuits  of  the  adjacent  interconnects. 
The  procedures  "sortwire”  and  "comparewire”  in  the  subprogram 
"exninth2.p*  perform  the  above  functions.  It  should  also  be  men¬ 
tioned  that  the  self-capacitances  of  an  interconnect  will  be  reduced 
slightly  due  to  another  interconnect  in  the  vicinity  as  discussed  in 
Sec.  3.3.1.  However,  this  effect  is  not  taken  into  account  in  the 
computation  of  the  self-capacitances  in  this  extractor. 
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WIRE  #1 


h 


WIRE  #2 


4,.4..£.  Report 


Basically,  The  extractor  generates  two  output  files.  One  is  the 
log  file  and  the  other  is  the  SPICE  file.  The  log  file  records  warn¬ 
ing  and  error  messages  and  all  the  steps  taken  in  the  extraction  pro¬ 
cess  along  with  the  time  and  date  for  reference.  The  original  CIF 
input  file  and  some  data  structures  like  wires  and  transistors  are 
also  included  for  debugging.  It  is  defaulted  to  "list. out"  or 
another  name  given  by  the  user.  A  detailed  explanation  of  the  usage 
and  commands  of  the  extractor  can  be  found  in  the  procedure  header  of 
"getcmdlnargs"'  in  the  subprogram  "exsecond.p" . 

Another  file  is  the  SPICE  file  which  includes  all  the  parameters 
of  the  transistors  and  the  interconnect  resistances  and  capacitances. 
The  format  on  this  file  is  compatible  with  the  SPICE  input  format. 
The  node  numbers  are  set  in  the  extraction  process  and  cannot  be  set 
externally.  This  may  need  to  change  in  the  future  version  so  that 
the  user  can  set  the  crucial  node  number  to  his  choice  to  improve  the 
readability  of  the  output.  This  file  is  defaulted  to  the  user's  ter¬ 
minal,  if  not  otherwise  redirected  or  renamed.  A  few  examples  will 
be  discussed  in  the  next  section  and  some  SPICE  files  are  shown 
there. 

£.£.  Examples  and  Results 

In  this  section,  a  few  examples  to  illustrate  the  network 
extractor  are  presented.  The  circuit  layouts  are  drawn  manually 
employing  the  design  rules  comformable  to  those  in  [38] .  After 
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supplemented  with  some  model  cards  and  input  output  specifications, 
the  extracted  outputs  of  the  circuits  are  fed  into  SPICE,  and  the 
simulated  circuit  responses  are  compared. 

Example  I  :  Static  RAM  cell 

In  Fig.  4.5,  the  layout  of  a  six-transistor  static  RAM  cell  is 
shown.  After  extraction  by  the  network  extractor,  the  output  SPICE 
file  is  given  in  Table  4.1.  "C0001*,  "rOOO"  and  "rrOOO"  are  the  three 
components  in  a  T  equivalent  circuit  for  an  interconnect.  "ROOO"  and 
"rrOOO"  have  the  same  valne  of  resistance.  The  center  node  of  the  T 
circuit  where  the  capacitors  are  connected  to  the  ground  is  numbered 
from  500.  The  capacitance  "c003"  is  the  only  coupling  capacitance 
between  the  interconnects  of  nodes  7,10  and  nodes  6,11.  That  is,  the 
coupling  capacitance  of  the  metal  and  polysilicon  cross-over  at  the 
center  of  the  layout. 

A  schematic  circuit  diagram  of  the  circuit  is  also  shown  in 
Fig.  4.6.  Note  that  the  negative  node  numbers  on  "motOOO"  and 
"mot002"  are  due  to  the  dangling  nodes  connected  to  the  bit  line  and 
the  word  line  on  the  pass  transistors  at  both  ends.  Those  negative 
nodes  are  helpful  in  locating  not  only  the  input  and  output  nodes, 
but  also  the  erroneous  disconnected  nodes.  It  should  also  be  men¬ 
tioned  that,  in  the  extraction  process,  the  parameters  "RLIMIT"  and 
"CLIMIT"  in  the  file  "constant. h"  are  specified  as  zero,  so  all  of 
the  resistances  and  capacitances  of  the  interconnects  are  computed 
and  listed  in  Table  4.1. 
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TABLE  4.1  SPICE  Output  Listing  for  SRAM  Coll. 


cOOO 

500  0 

O.OllOpf 

rOOO 

8  500 

2.9C+02 

rrOOO 

500  9 

2.9o+02 

cOOl 

501  0 

0.0041pf 

rOOl 

10  501 

l.le+02 

rrOOl 

501  7 

l.le+02 

c002 

502  0 

0.0077pf 

r002 

11  502 

7. 5o+01 

rr002 

502  6 

7.5o+01 

c003 

501  502 

0.0020pf 

mldOOO 

5  6 

6  0  load 

1- 

9.00ua 

2.67um 

mid 001 

5  7 

7  0  load 

1- 

9 . OOum 

w*» 

2.67mm 

mot 000 

7  8 

-1  0  pass 

1- 

2 .OOum 

2. OOum 

motOOl 

6  10 

0  0  dr  ive 

1- 

2. OOum 

wm 

8. OOum 

mot002 

-2  9 

6  0  pots 

1- 

2. OOum 

2. OOum 

mot 003 

7  11 

00  dr ive 

1- 

2. OOum 

w- 

8. OOum 

1.1.1.  Xi  :  Full  Adder  Circuit  I 

The  scheme tic  circuit  diagram  of  "Full  Adder  Circuit  I"  ia  shown 
in  Fig.  4.7.  Vith  the  same  specif icatious  in  the  " constant. h"  file 
as  in  the  previous  example,  we  can  extract  the  SPICE  file  from  the 
extractor.  Combined  with  some  model  cards,  input  and  output  specifi¬ 
cations.  and  a  few  modifications  to  the  circuit,  a  few  simulations* 
of  the  circuit  are  performed.  In  order  to  distinguish  the  effects  of 
the  interconnect  parameters,  three  cases  are  considered,  namely  :  (i) 
the  circuit  without  taking  into  account  the  interconnect  parameters; 
(ii)  the  circuit  with  only  the  self-capacitance  and  the  resistance  of 
the  interconnects;  (iii)  the  circuit  with  all  interconnect  parame¬ 
ters  including  the  coupling  capacitances.  These  three  simulated 
responses  of  "Full  Adder  Circuit  I"  are  shown  in  Fig.  4.8  to 
Fig.  4.10.  The  full  SPICE  deck  for  case  (iii)  is  given  in 
Appendix  E.  It  should  be  pointed  out  that  all  eight  possible  combi¬ 
nations  of  the  three  inputs  to  the  full  adder  are  covered  in  the 
sismlationa,  and  only  the  first  input  waveform  is  shown  in  the 
graphs.  The  other  two  input  waveforms  are  not  plotted  for  the  sake 
of  neatness  of  the  graph. 

In  comparing  Fig.  4.9  with  Fig.  4.8,  it  is  obvious  that  the  rise 
and  the  fall  times  in  Fig.  4.9  are  larger  than  those  in  Fig.  4.8. 
This  is  naturally  due  to  the  SC  time  delay  contributed  by  the  inter¬ 
connects  in  the  circuits.  It  should  be  mentioned  that  the  time 


*  These  simulations  are  actually  done  by  SLATE  [48,49]  which  has  a 
similar  format  as  SPICE  and  does  a  better  job  at  transient  analysis. 
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delays  shown  in  Fig.  4.8  may  be  attributed  to  the  capacitive  loads  of 
0.1  pF  at  both  output  ports.  Furthermore,  these  loads  are  used 
throughout  all  the  circuits  discussed  iu  these  two  examples  of  full 
adder  circuits. 

The  circuit  response  of  the  third  case  which  also  includes  the 
coupling  capacitances  are  shown  in  Fig.  4.10.  They  differ  little 
from  the  responses  shown  in  Fig.  4.9.  On  close  examination,  there 
may  be  just  a  little  more  delay.  Overall,  it  can  be  said  that  the 
coupling  capacitances  for  this  circuit  do  not  influence  the  circuit 
responses  much.  However,  the  inclusion  of  the  self-capacitances  and 
the  resistances  has  more  visible  effects  on  the  responses,  and  they 
should  be  considered  in  the  circuit  simulation. 

Example  III  :  Full  Adder  Circuit  II 

Another  full  adder  circuit  is  used  as  an  example.  Its  schematic 
is  shown  in  Fig.  4.11.  There  is  a  basic  design  difference  from  the 
previous  full  adder  circuit.  In  this  case,  more  pass  transistors  are 
used.  Noticeably,  they  are  used  at  the  output  ports;  thus,  the  out¬ 
put  responses  of  this  circuit  cannot  reach  the  full  5  volts  of  VDD. 
At  most,  the  logic  "l"  for  the  outputs  can  reach  5  volts  minus  the 
threshold  voltages  of  the  pass  transistors  which  are  set  at  1  volt. 
This  point  is  manifested  in  the  following  graphs.  Those  three  cases 
mentioned  in  the  previous  example  are  also  considered  here.  In 
Fig.  4.12,  the  responses  of  "Full  Adder  Circuit  IIM  without  the 
interconnect  parameters  are  shown.  Figure  4.13  is  for  the  case  of 
taking  the  self-capacitances  and  the  resistances  into  account. 


f  s  s 
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Finally,  in  Fig.  4.14  the  responses  with  all  the  interconnect  param¬ 
eters  included  are  considered.  The  "constant. h"  file,  the  model 
cards,  and  the  load  capacitances  are  chosen  the  same  as  in 
Example  II;  the  listing  of  the  SPICE  input  deck  is  given  in 
Appendix  F. 

In  comparing  the  responses*  for  these  three  cases,  it  is  obvious 
that  the  coupling  capacitances  in  this  circuit  have  greatly  affected 
the  output  responses.  Especially,  the  waveforms  in  Fig.  4.14  differ 
drastically  from  those  in  Fig.  4.12  and  Fig.  4.13.  Not  only  are  the 
rise-  and  the  fall-time  delays  of  the  responses  larger  than  those  in 
Fig.  4.13  and  Fig.  4.12,  but  also  there  is  a  notch  developed  in  the 
sum  waveform.  This  may  result  in  unwanted  errors  or  hazards  for  the 
circuit.  Hence,  the  circuit  or  the  layout  of  the  circuit  should  be 
examined  more,  or  redesigned,  to  correct  this  possible  flaw.  In 
design  verification,  it  is  our  aim  to  discover  and  eliminate  those 
possible  errors.  On  the  other  hand,  the  responses  in  Fig.  4.13  with 
only  the  self-capacitances  considered  differ  from  those  in  Fig.  4.12 
by  larger  rise-  and  fall-time  delays  as  the  case  in  the  previous 
example . 

To  sum  up,  it  can  be  concluded  that  the  coupling  capacitances  of 
the  interconnects  are  important  in  the  circuit  simulation  in  this 
example.  However,  in  the  previous  example,  they  had  an  insignificant 
effect  on  the  circuit  responses.  Consequently,  coupling  capacitance 
effects  depend  heavily  on  the  mask  layout  and  the  circuit  itself. 


*  These  simulations  are  also  done  by  SLATE. 
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Their  effects  on  the  circnit  have  to  be  investigated  individually. 
As  to  the  self-capacitances  and  the  resistances  of  the  interconnects, 
they  sorely  contribute  to  the  time  delay  of  the  circuit  and  should  be 
considered. 

Computational  Complexity  of  the  Extractor 

Because  of  the  increased  circuit  complexity  and  myriads  of  data 
involved  in  VLSI  circuits,  it  is  essential  that  the  CAD  tools  consist 
of  computationally  efficient  algorithms.  The  insertion  algorithm  of 
the  4-d  binary  search  tree  with  right  thread,  which  is  the  basic  data 
structure  of  the  extractor,  is  given  partly  in  [26]  and  listed  as 
follows  : 

Given  a  node  P  to  be  inserted  into  the  tree. 

11.  [Check  for  empty  tree]  If  ROOT  ■  nil,  then  ROOT  <-  P, 

RSON(P)  <-  nil,  LSON(P)  <-  nil,  set  RSON(P)  as  a  thread, 

return;  otherwise,  Q  <-  ROOT. 

12.  [Comparison]  If  K± (P)  -  K^Q),  0£  i  £  3,  and 

i  is  the  discriminator,  then  return;  otherwise,  set  SON(Q) . 
(RSON  or  LSON  depending  on  the  comparison  of  Key(P)  and 
Key(Q) ) .  If  SON(Q)  -  nil,  then  goto  14. 

13.  [Move  down  the  tree]  Set  Q  <-  SON(Q) ,  goto  12. 

14.  [Insert  new  node  in  tree]  Set  SON(Q)  <-  P,  RSON  (P)  <-  nil, 

LSON(P)  <-  nil,  return. 


A  detailed  discussion  of  the  complexity  of  the  above  algorithm  was 
given  in  [26],  and  it  was  concluded  that  typical  insertions  and 
record  look-ups  in  a  4-d  tree  will  examine  approximately  1.3861og,n 
nodes,  where  n  is  the  size  of  the  tree.  In  the  extractor,  this 


insertion  algorithm  is  called  to  set  np  a  4-d  binary  search  tree  for 
rectangles  on  each  mask.  It  is  also  set  np  for  the  diffusion  and 
polysilicon  overlap  rectangles. 

Next  in  search  of  the  diffusion  and  polysilicon  overlap  regions, 
a  search  algorithm  is  needed.  For  example,  given  a  polysilicon  rec¬ 
tangle  P  with  its  four  keys,  ve  have  to  find  all  the  diffusion  rec¬ 
tangles  Q  which  are  overlapped  with  P.  In  other  words,  the  keys  of  P 
and  Q  have  to  satisfy  all  the  following  conditions  simultaneously  : 

*.<Q)  i  -*.( P) 

Ki«W  i 

*s«M  1  -Ie(P) 

*,<Q)  i  -Ks(P) 

This  search  is  classified  as  "Intersection  Queries"  in  [26,30].  The 
search  algorithm  can  be  described  as  follows  : 

Given  a  rectangle  P,  try  to  find  all  the  overlapping  rectangles  Q' s 
in  the  tree  rooted  by  SOOT. 

11.  [Initialization]  Q  <-  SOOT. 

12.  [Move  to  leftmost  corner]  While  LSON(Q)  t  nil  do 

Q  <-  LSON(Q)  and  INCREMENT* i) ,  i  is  the  discriminator. 

13.  [Check  intersection]  If  I.(Q)  £  -  K, j+2)mod4<P> » 

0  <  j  £  3 ,  then  report  Q,  continue  14. 

14.  [Go  to  predecessor]  If  K.(Q)  >  -  ] 

SSON(Q)  is  a  thread,  then  set  Q  <- 
If  Q  ■  nil,  then  return;  otherwise  goto  12. 


&ii££ctil§6s(Q)  !r 


The  above  algorithm  is  similar  to  that  described  in  [27];  basically. 


it  is  tn  "inorder"  traverse  of  the  binary  search  tree.  The  main 
advantage  that  this  data  structure  offers  in  saving  traverse  time  is 
due  to  the  check  on  step  14.  If  the  condition  E^(Q)  >  “*(  i+2)nod4^ 
is  satisfied,  then  the  right  subtree  of  Q  ueed  not  be  traversed,  and 
the  search  can  be  directly  positioned  to  the  predecessor  of  Q.  Note 
that  only  one  condition  is  checked  in  this  step,  and  i  is  the 
discriminator.  Instead,  all  four  keys  have  to  be  examined  in  step 
13.  Moreover,  the  procedure  to  find  the  predecessor  can  be  accom¬ 
plished  rather  easily  with  the  help  of  the  right  threads  woven  in  the 
tree  instead  of  employing  a  stack  register  to  store  the  nodes  on  the 
way  down  the  tree  as  suggested  in  [27].  Hence,  this  is  another  sav¬ 
ing  in  the  search  time.  However,  due  to  this  step  in  the  search 
algorithm  and  the  random  nature  of  the  tree,  it  is  difficult  to 
evaluate  the  computational  complexity  for  this  algorithm.  Some 
empirical  data  had  been  collected  to  test  the  complexity  as  discussed 
in  [27] .  It  may  be  inferred  that  the  search  algorithm  of  the  4-d 
binary  tree  is  rather  efficient. 

Other  possible  expensive  algorithms  in  the  extractor  may  be 
those  for  establishing  the  interconnect  data  structure  wire  and  find¬ 
ing  the  coupling  capacitances.  Because  the  transistors  are  stored  in 
linear  arrays,  a  search  through  the  whole  arrays  has  to  be  done  for 
each  interconnect  rectangle  to  check  if  it  is  connected  to  the 
transistor  ports.  Also,  the  algorithm  of  finding  the  coupling  capa¬ 
citances  is  done  by  first  comparing  the  wires  pairwise,  though  "wire- 
bounds"  are  used  first  to  screen  off  far-away  candidates;  then,  all 


the  branches  inside  the  vires  are  compared  pairvise.  Both  of  these 
algorithms  are  approximately  n*  processes,  and  they  can  be  computa¬ 
tionally  intensive.  Thus',  they  are  in  want  of  some  improvements  in 
the  future  versions. 

In  order  to  experimentally  estimate  the  computational  complexity 
of  the  extractor,  a  fev  cases  of  a  chain  of  ''Full  Adder  Circuit  II”, 
presented  in  the  previous  section,  are  tested.  The  number  of 
transistors  in  the  circuit  is  used  as  a  measure  of  the  complexity  of 
the  circuit.  The  results  are  given  in  Table  4.1.  An  empirical  rela¬ 
tion  between  the  number  of  transistors  (N)  and  CPU  time  (T)  can  be 
assumed  as  :  T  *  oN0.  From  the  data  in  Table  4.1,  the  exponent  £  can 
be  estimated  to  be  in  the  range  of  1.3  to  1.7.  Therefore,  the  compu¬ 
tational  complexity  of  the  extractor  is  proportional  to  0(Nlal)  to 
0(N**  ),  vhere  N  is  the  total  number  of  transistors.  This  is 


TABLE  4.1  CPU  Time  Consumption  of  the  Extractor. 


no.  of  transistor  (N) 

CPU  Time  (T)  in  sec 

21* 

13.3 

42 

33 

63 

58.1 

84 

84 

105 

122.4 

168 

264 

336 

851.2 

*  One  "Full  Adder  Circuit  II"  contains  21  transistors. 
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comparable  to  some  other  extractors,  e.g.,  the  extractor  discussed  in 
[28]  has  complexity  of  OCN1**). 

In  summary,  from  the  examples  shown  in  this  chapter,  it  can  be 
concluded  that  the  inclusion  of  interconnect  parameters  in  the  cir¬ 
cuit  simulation  is  important.  One  can  discover  possible  timing 
errors  in  the  design  or  the  layout.  In  order  to  handle  more  complex 
circuitry,  the  algorithms  employed  in  the  extractor  need  to  be  effi¬ 
cient,  and  hierarchy  should  be  built  into  the  program. 


CONCLUSION 


The  complexity  of  VLSI  circuits  has  necessitated  the  development 
of  sure  comprehensive  CAD  systems.  One  of  the  problems  addressed  in 
this  thesis  is  concerned  vith  the  interconnection  circuit.  It  is 
essential  to  take  into  account  the  electrical  parameters  of  the 
interconnects  in  the  more  complicated  VLSI  circuits  as  they  gradually 
become  the  dominating  factors.  Furthermore,  it  should  also  be  "com - 
putationally  feasible"  to  extract  this  information  from  the  layout  in 
the  design  verification.  Hence,  there  are  basically  two  parts  in 
this  thesis.  The  first  part  is  concerned  vith  the  accurate  computa¬ 
tion  of  the  interconnect  capacitances  vith  numerical  methods,  more 
specifically,  the  integral  method  vas  employed  here.  Then  simple 
formulas  for  the  interconnect  capacitances  are  developed  in  the 
thesis  vith  a  viev  to  easily  incorporate  them  into  an  extractor. 
The  second  part  of  the  thesis  deals  vith  the  extraction  of  intercon¬ 
nect  parameters,  and  a  netvork  extractor  vas  developed  for  it. 

In  Chapter  2  a  detailed  formulation  of  the  integral  method  vas 
given.  It  encompasses  the  derivation  of  the  Green’s  functions  in 
homogeneous,  tvo  layers  and  three  layers  of  media;  the  even-odd  mode 
analysis  for  multi-conductor;  and  the  use  of  the  "method  of  moments" 
to  solve  the  integral  equation.  Pulse-type  functions  vere  used  for 
both  the  basis  functions  and  the  testing  functions  for  the  sake  of 
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closed-form  evaluation  of  kernel  matrix  elements.  This  method  is 
integrated  into  tiro  FORTRAN  programs  "CAP2D"  and  "CAP3D".  All  the 
self  capacitances  with  respect  to  the  ground  of  every  interconnect  in 
question  and  all  the  coupling  capacitances  between  interconnects  can 
be  evaluated  up  to  three  layers  of  media  and  at  most  ten  intercon¬ 
nects.  For  the  two-dimensional  case,  slant  side  walls  of  the  inter¬ 
connects  were  also  considered.  As  the  examples  in  Chapter  2  and 
Appendix  B  show,  a  variety  of  interconnect  capacitances  can  be 
evaluated  by  these  programs,  and  they  are  reasonably  efficient  and 
versatile.  However,  in  three  layers  of  media,  the  Green's  function 
has  become  quite  complicated  and  the  computation  may  be  intensive. 
If  detailed  capacitance  information  is  needed  for  irregularly  shaped 
boundaries,  the  integral  method  may  not  be  suitable. 

As  for  the  si^le  formulas,  the  "cylindrical  approximation  for¬ 
mula"  was  derived  in  Chapter  3  for  the  self-capacitance.  It  is  based 
on  the  known  formula  of  the  capacitance  of  a  cylinder  above  a  ground 
plane.  As  shown  in  Chapter  3,  it  compares  favorably  with  other  sim¬ 
ple  formulas  and  is  reasonably  accurate  with  respect  to  the  numerical 
results.  Then  a  "least  square  fit"  method  was  employed  to  find  sim¬ 
ple  formulas  for  the  coupling  capacitances  between  two  parallel 
interconnects  or  two  interconnects  on  different  levels  which  cross 
over  each  other.  These  prove  to  be  useful  in  the  network  extractor. 
As  the  interconnect  structure  becomes  more  complicated,  simpler  for¬ 
mulas  may  be  necessary  to  include  all  the  capacitances  in  the  extrac¬ 


tion 
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Finally,  a  detailed  description  of  the  network  extractor  was 
given  in  Chapter  4.  It  is  a  rectangle-based,  Manhattan-geometry- 
only,  NMOS  extractor.  It  reads  in  CIF  input  of  the  layout  informa¬ 
tion,  then  pnts  ont  circuit  components  in  SPICE  compatible  form.  The 
SPICE  output  file  includes  transistors  and  their  parameters,  e.g., 
channel  width  and  channel  length  and  the  resistance  and  capacitance 
of  the  T-equivalent  lumped  circuit  of  the  interconnect.  A  few  exam¬ 
ples  were  given  in  Chapter  4,  along  with  some  circuit  simulations. 
From  these  examples,  it  is  apparent  that  the  interconnect  parameters 
certainly  play  an  important  role  in  the  circuit  simulation.  Although 
the  original  design  did  not  emphasize  the  efficiency  of  the  extrac¬ 
tor,  it  was  found  out  by  experiment  that  the  computational  complexity 
of  the  extractor  is  approximately  proportional  to  0(N**’),  where  N  is 
the  number  of  transistors  in  the  circuit.  However,  due  to  simplifi¬ 
cation  of  the  transistor  extraction,  only  a  limited  repertoire  of 
transistor  configurations  is  recognized  in  the  extractor.  Also,  some 
redundancy  in  the  data  structure  and  possible  inefficiencies  in  the 
algorithms  may  need  further  improvements.  In  all,  it  fulfills  the 
need  to  have  the  interconnect  parameters  extracted  and  utilizes  the 
simple  formulas  derived  in  Chapter  3.  Besides  increasing  the  reper¬ 
toire  of  transistor  configurations  recognizable  to  the  extractor  and 
improving  some  of  the  algorithms  and  data  structures  in  the  program, 
it  is  also  necessary  to  link  the  extractor  hierarchically  with  other 
CAD  tools  to  handle  larger  circuits. 
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THE  VALIDITY  OF  THE  LUMPED  CIRCUIT  MODEL 


Introduction 

The  lumped  circuit  model  for  the  interconnect  is  employed  in 
this  thesis.  A  discussion  is  presented  here  to  establish  the  vali¬ 
dity  of  this  approach  through  the  transmission  line  theory. 

Formulation 

Consider  a  transmission  line  circuit  and  its  T-equivaleut  cir¬ 
cuit  shovn  in  Fig.  A.l  and  Fig.  A. 2.  The  length  of  the  line  is  ”1" 
snd  **I  "  and  MZj”  are  the  current  source  and  the  line  load,  respec¬ 
tively.  The  characteristic  impedance  and  the  propagation  constant 

are 

Z,  -  [(R+j«L)/(G+j«C)]1/l 
r  -  t(R+j«L) .(G+j»C)]l/* 

vhere  R,  L,  G,  C  are  the  resistance,  inductance,  conductance  and 
capaoitance  per  unit  length  of  the  line.  From  the  transmission  line 
theory,  the  input  impedance  and  the  voltage  response  at  the  load 
end  of  the  line  Vj  can  be  expressed  as  [46] 


ZjCosMyD  +  Z0sinh(  yD 


(A.l) 


Z 

s 


Z. 

ZjjinJi(Yl)  +  Z0cosh(Yl) 


Wi 

Z^sinh(Yl)  +  Z#cosh(Yl) 


(A. 2) 


These  two  parameters  define  the  characteristics  of  the  transmission 
line,  and  they  will  be  compared  with  those  of  the  "equivalent"  lumped 
circuit  model. 

Next  consider  the  T-eqnivalent  lumped  circuit  in  Fig.  A. 2.  All 
the  parameters  are  defined  the  same  as  those  in  Fig.  A.l.  Then  the 
input  impedance  Z'  and  the  voltage  response  of  this  lumped  circuit 
can  be  easily  obtained  as 


Z' 

s 


ZjU+O.Sy*!*)  +  Zo(y1+0.25y*  1* ) 

Zo(1+0.5yV)  +  ZjCyD 


(A. 3) 


Wi 


Z,(1+0.5y  1  )  +  Z.(T1) 


(A. 4) 


In  order  to  show  that  the  lumped  circuit  in  Fig.  A. 2  is  an 
"equivalent"  circuit  to  the  transmission  line  circuit  in  Fig.  A.l, 
the  input  impedance  and  the  voltage  response  should  be  correspond¬ 
ingly  equal.  In  other  words,  Eq.  (A. 3)  should  be  equal  to  to 


Eq.  (A.l)  and  Eq.  (A. 4)  to  Eq.  (A. 2) .  Note  that  if  the  hyperbolic 
functions,  sinh  and  cosh,  in  Eq.  (A.l)  and  Eq.  (A. 2)  are  expanded  in 
a  Taylor  series,  the  first-order  approximations  are  : 

sinh(yl)  -  yl  +  . . . . 

cosh(yl)  ~  1.0  +  0.5y*l*  + - , 

Furthermore,  if  the  above  approximations  are  substituted  into 
Eq.  (A.l)  and  Eq.  (A. 2) ,  Eq.  (A. 3)  and  Eq.  (A. 4)  are  obtained  except 
for  the  term  0.25y  1  in  Eq.  (A. 3) .  Therefore,  it  may  be  concluded 
that  the  equivalent  lumped  circuit  is  valid,  as  long  as  the  parameter 
yl  is  small  enough. 

A.  3_.  Discussion 

A.i-i,.  Lossless  transmission  line  First,  let  us  consider  a  loss¬ 
less  transmission  line.  In  this  case,  the  resistance  and  the  conduc¬ 
tance  of  the  line  are  relatively  small,  i.e.,  R  <<  u»L,  and  G  <<  u>C . 
Then  the  propagation  constant  becomes 

2n 

r  -  jp  -  j - 

X 

vhere  X  is  the  wavelength  of  the  propagating  wave  on  the  transmission 
line.  In  order  to  satisfy  the  criterion  arrived  earlier,  i.e.,  the 
factor  yl  -  j2nl/X  must  be  small,  the  length  of  the  line  "1"  should 
be  small  compared  with  the  wavelength  on  the  line  X.  This  is  the 
usual  criterion  for  the  validity  of  the  lumped  circuit  model. 
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A. 3,. 2.  RC  lossy  line 

Next,  consider  a  lossy  line  with  negligible  inductance  and  con¬ 
ductance,  i.e.,  &  >>  wL,  and  G  <<  uC.  In  this  case  the  propagation 
constant  becomes 


r 


<l+j)-< 


uRC 

2 


1/» 


To  keep  yl  small  as  discussed  before,  the  factor  (uRC)1/*!  should  be 

a  i 

kept  small  too.  In  order  words,  u  RC1  should  be  small.  The  term 
RC1*,  as  a  whole,  can  be  considered  as  the  "time  delay'*  of  the 
transmission  line,  and  the  factor  1*  is  due  to  the  per  unit  length 
basis  of  the  resistance  and  the  capacitance.  If  T  is  the  period  of 
the  propagating  wave,  and  u  »  2tt/T,  then  u»RCl  *  2ji(RC1*/T)  .  The 
criterion  for  the  validity  of  the  lumped  circuit,  i.e.,  small  yl,  can 
be  interpreted  in  this  case  as  implying  that  the  delay  of  the  line 
should  be  small  compared  with  the  period  of  the  wave.  For  a  typical 
example  [47J ,  a  doped  polysilicon  line  with  30Q/square  and  1  cm  long 
has  a  time  delay  in  the  order  of  10  sec.  This  corresponds  to  a  pro¬ 
pagating  wave  on  the  line  of  frequency  about  10  MHz.  This  is  usually 
unacceptable  for  the  typical  circuit  employing  silicon  technology. 
Therefore,  long  lines  usually  lie  on  the  metal  level.  For  the  case 
of  large  R  or  C  of  the  line,  multiple  sections  of  a  lumped  circuit 
can  be  used  to  approximate  the  transmission  line.  In  summary,  for  RC 
lossy  lines,  the  criterion  for  the  validity  of  the  lumped  circuit  is 
that  the  RC  time  delay  of  the  line  be  small  relative  to  the  period  of 
the  corresponding  frequency  of  the  line. 
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APPENDIX  B 


CAP  USER' S  MANUAL 


B.i.  Introduction 

"CAP2D*  and  "CAP3D1*  are  two  FORTRAN  programs  for  two-dimensional 
and  three-dimensional  capacitance  computations.  A  detailed  descrip¬ 
tion  of  their  usage  and  output  format  is  given;  some  examples  are 
also  given  at  the  end. 

B.2..  Description 

The  numerical  computation  of  capacitance  of  these  two  programs 
is  based  on  the  integral  method  which  is  essentially  the  Green's 
function  approach  combined  with  the  moment  method.  Pulse-type  basis 
functions  over  subdivisions  on  the  conductor's  surface  are  assumed, 
and  the  testing  functions  are  chosen  the  same  as  the  basis  functions. 
The  subdivisions  can  be  chosen  to  be  variable-length  or  constant- 
length.  If  variable-length  subdivisions  are  chosen,  by  default,  they 
are  divided  according  to  the  roots  of  the  Chebyshev  polynomials. 
Since  subdivisions  with  variable  length  usually  yield  more  accurate 
results  for  the  same  number  of  subdivisions,  the  constant  case  is 
basically  included  for  comparison  purposes. 

Both  programs  have  the  capability  to  handle  up  to  10  traces  (or 
conductors)  in  at  most  3  layers  of  different  media.  In  the  two- 
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dimensional  case,  "CAP2D"  is  also  capable  of  taking  into  account  the 
pitch  angle  of  the  slant  side  walls,  assuming  that  the  traces  are 
left  and  right  symmetric.  However,  "CAP3D"  for  the  three-dimensional 
case  deals  with  vertical  side  walls  only.  In  the  case  of  multiple 
layers,  i.e.,  two  layers  or  three  layers  of  media,  an  infinite  summa¬ 
tion  is  involved  in  the  Green's  function  and  a  cut-off  criterion 
should  be  chosen. 

B.2.J..  Usage 

The  programs  are  designed  to  be  run  interactively.  Alterna¬ 
tively,  the  UNIX  system  also  provides  a  way  to  run  it  through  an 
input  file  and  redirect  ("<”)  command.  The  following  parameters  will 
be  proaipted  in  running  the  programs. 

(1)  "number  of  traces"  and  "number  of  layers"  they  are  self- 
explanatory. 

(2)  "want  kernel  matrix  ?"  and  "want  grids  ?"  :  If  "y"  is  chosen, 
the  kernel  matrix  elements  and  the  coordinates  of  all  the  subdi¬ 
visions  will  be  listed  in  the  output  file.  These  are  basically 
intended  for  debugging  purposes.  Usually  "n"  should  be  the 
answer  to  both  requests. 

(3)  "wx.nx”  :  These  are  the  horizontal  or  x-direction  width  of  the 
trace  and  its  corresponding  number  of  subdivisions.  Note  that 
all  the  widths  and  lengths  are  stored  in  double  precision  vari¬ 
ables,  so  "wx*  should  be  in  double  precision  form.  This  also 
applies  to  all  the  other  widths  in  the  program. 
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(4)  "wy,ny" ,"wz,nz"  :  These  parameters  are  similar  to  the  previous 
ones,  except  they  are  in  vertical  or  y  direction.  "wz.nz"  only 
appears  in  the  three-dimensional  case. 

(5)  "ang"  :  This  is  the  pitch  angle  for  the  slant  side  vails.  It 

only  occurs  in  the  two-dimensional  case.  The  angle  should  be  in 
degrees,  and  within  the  range  of  0°  to  90°. 

(6)  "subdiv  method"  :  This  concerns  the  choice  of  constant-  or 
variable-length  subdivisions.  A  "0"  means  constant  subdivisions 
and  any  other  integer  will  result  in  defanlt  Chebyshev  variable 
subdivisions. 

(7)  "orientation"  :  This  only  occurs  in  running  "CAP3D" .  It  speci¬ 
fies  the  direction  of  the  trace.  For  example,  a  trace  which  is 
parallel  to  the  x-axis  has  an  orientation  1.  By  the  same  token, 
the  orientation  is  2  if  it  is  parallel  to  the  y-axis.  However, 
orientation  3  is  prohibited.  It  should  be  pointed  out  that  the 
charges  on  the  cross  sections  perpendicular  to  the  orientation 
axis  of  the  trace  are  neglected,  i.e.,  charges  at  each  end  of 
the  line  are  not  computed. 

(8)  "lolx.loly"  :  In  "CAP2D",  the  coordinates  of  the  lower  left 
corner  of  the  trace  are  needed. 

(9)  "cx,cy,cz"  :  On  the  other  hand,  the  coordinates  of  the  center  of 
the  trace  are  used  in  "CAP3D". 

Note  that  the  above  parameters  from  (3)  to  (9)  have  to  be  fur¬ 
nished  for  every  trace  in  question. 
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(10}  "enl,en2"  :  These  are  for  the  dielectric  constants  of  the  media. 
Because  the  outmost  layer  in  the  three-layer  case  is  assumed  to 
be  free  space,  at  most  two  different  dielectric  constants  need 
to  be  specified.  It  is  important  to  note  that  the  permittivity 
of  the  free  space  (8.854x10  14  F/cm)  is  not  built  into  the  pro¬ 
grams.  Thus,  usually  the  output  capacitance  values  are  normal¬ 
ized  by  this  factor.  If  capacitance  values  in  Farad  or  F/cm, 
depending  on  three-  or  two-dimensional  problem,  are  desired, 
then  the  permittivity  of  the  free  space  should  be  included  in 
the  dielectric  constants. 

(11)  "a.b.h"  :  These  are  the  heights  of  the  dielectric  layers.  They 
should  be  double  precision  numbers.  MaM  and  "b"  are  for  the 
case  of  three  layers  of  media,  and  it  is  assumed  that  b  >  a. 
WH"  is  the  height  of  the  lower  dielectric  layer  in  the  two-layer 

case. 

(12)  "errbnd"  :  This  is  the  cut-off  criterion  for  the  infinite  summa¬ 
tion  in  the  two— layer  or  three-layer  Green's  functions,  l.d-3 
is  recoeuaended  in  the  program.  7et,  bear  in  mind  that  a  small 
"errbnd”  will  render  more  computations,  consume  more  CPC  time, 
and  result  in  more  accurate  answers.  This  trade-off  between 
accuracy  and  computation  is  also  observed  in  the  case  of  the 


choice  of  the  number  of  subdivisions  for  the  traces 
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B.2.2,.  Output 

The  result  of  the  capacitance  computation  will  be,  by  default, 
stacked  in  files  "fort. 7"  for  "CAP2D"  and  "fort. 9"  for  "CAP3D" .  If 
these  files  already  exist  in  the  directory,  the  current  output  will 
be  appended  to  the  corresponding  files.  The  report  includes  all  the 
pertinent  input  parameters  for  the  traces  and  the  media.  Further- 
sure,  the  following  parameters  are  also  included  : 

"Number  of  kernels"  indicat  s  the  number  of  elements  in  the 
upper  triangular  part  of  the  kernel  matrix  including  diagonals.  If 
the  total  number  of  subdivisions  is  n,  then  the  "number  of  kernels" 
will  be  n*(n+l)/2.  Note  also  that  the  kernel  matrix  will  be  listed 
in  "packed  symmetric"  form  if  it  is  requested.  In  essence,  the  upper 
triangular  and  the  diagonal  elements  are  listed  linearly.  For  exam¬ 
ple  ,  the  (i.j)th  element  (i  <  j)  in  the  kernel  matrix  will  be  the 
[j*( j”D /2+i]th  element  in  the  linear  list. 

"Iterm"  indicates  the  maxiaram  number  of  terms  taken  in  the  com¬ 
putation  of  multi-layer  kernels.  It  is  dependent  on  the  cut-off  cri¬ 
terion  "errbnd"  as  discussed  previously.  And  it  is  0  for  the 
single-layer  case.  Since  the  do  loop  index  for  the  three-layer  case 
is  set  at  12,  it  has  not  neccessarily  converged  if  "iterm"  equals  13. 
If  divergence  is  encountered  in  the  summation  process,  i.e.,  the 
absolute  value  of  the  term  in  the  infinite  summation  is  larger  than 
its  predecessor,  "no.  of  divgnt"  will  be  incremented,  and  the  summa¬ 
tion  will  be  terminated.  Since  this  divergence  usually  occurs  at 
terms  with  small  magnitude  (sl.d-2),  the  current  sum  can  still  be  a 


reasonable  value  for  the  kernel  element,  and  the  computation  of  other 
kernel  elementa  will  continue. 

The  parameters  "rcond"  and  "info"  are  indicators  from  the  "lin- 
pack"  programs  "dppco"  and  "dppal"  which  are  utilized  to  solve  the 
positive  definite  kernel  matrix.  If  "info"  is  other  than  0,  the 
solution  is  erroneous.  It  indicates  the  order  of  the  principal  sub¬ 
matrix  which  is  not  positive  definite.  Usually  the  subdivision  at 
fault  can  be  detested  this  way.  "Rcond"  is  the  reciprocal  of  the 
condition  number  of  the  kernel  matrix.  Only  when  it  is  in  the  order 
of  the  machine  epsilon  (about  l.d-14  for  double  precision  on  VAX 
11/780) ,  then  is  the  linear  system  ill-behaved  and  the  solution  may 
be  erroneoua.  Usually  this  number  is  within  0.1  and  l.d-5. 

The  solution  vector  included  in  the  report  represents  the  even¬ 
mode  charges  on  the  subdivisions,  i.e.,  the  charge  distribution  when 
all  the  traces  are  at  a  constant  positive  potential  with  respect  to 
the  ground  plane.  It  delineates  the  charge  distribution  in  the  sys¬ 
tem  and  is  helpful  in  gaining  a  physical  understanding  of  the  sys¬ 
tem.  Also  it  serves  as  a  good  gauge  for  debugging. 

All  the  self  capacitances,  C  ,  and  the  coupling  capacitances, 
C  ,  n  j*  m  are  reported.  Note  that  the  accuracy  of  the  computed 
capacitances  depends  on  the  cutoff  criterion  "errbnd",  the  number  of 
subdivisions  and  the  subdivision  method.  Also  it  has  been  shown  in 
the  literature  that  the  integral  method  will  always  result  in  a  lower 
bound  for  the  true  capacitance  value,  and  the  true  value  is  also 
variationally  stationary  for  the  first  order  perturbation.  Thus, 
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sore  subdivisions  and  a  tighter  "errbnd"  will  yield  larger  capaci¬ 
tance  values.  This  characteristic  is  also  very  helpful  for  debug¬ 
ging,  and  for  the  estimation  of  the  true  value.  Generally  in  order 
to  ensure  the  accuracy  of  the  capacitances,  the  same  case  can  be  run 
with  successively  more  subdivisions  on  the  dimensions  of  the  traces. 

Finally,  the  time  usage  shows  the  consumption  of  CPU  time  for 
various  computations.  "Kernel  set  up  time"  is  the  time  to  compute 
the  kernel  elements.  It  usually  takes  up  the  lion’s  share  of  the 
total  user  time.  Since  the  "number  of  kernels"  is  the  total  number 
of  elements  computed,  it  is  directly  responsible  for  the  "kernel 
setup  time".  "Kernel  solve  time"  is  the  time  spent  by  "dppco"  and 
"dppsl"  to  solve  the  corresponding  positive  definite  kernel  matrix. 
"Capacitance  time"  is  the  time  spent  to  obtain  the  various  self  and 
coupling  capacitances  from  the  solution  vectors  (i.e.,  charge  distri¬ 
butions)  of  even-mode  and  odd-mode  excitations.  This  time  is  usually 
the  smallest  and  is  negligible  for  single  traces. 

£.1.  F_Ue? 

The  FORTRAN  source  files  are  "cap3d.f",  "subp3d0.f", 
"subp3dl.f",  "subp3d2.f",  "subp3d3.f"  for  the  three-dimensional  case. 
The  main  program  is  in  "cap3d.f",  all  the  others  are  subprograms. 
"subp3dl.f"  contains  the  subprograms  for  the  homogeneous  case, 
"subp3d2.f"  for  the  two-layer  case  and  so  on.  The  "make"  file  for 
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"cap3d.ex"  is  as  follows  : 


# 

#  make  file  for  cap3d.f 

# 

FFLAGS  ■  -g  -u 

cap3d.ex  :  cap3d.o  subp3d0.o  subp3dl.o  subp3d2.o  subp3d3.o 
coaptiae .o 

f77  -g  cap3d.o  subp3d0.o  sabp3dl.o  subp3d2.o  subp3d3.o 
coaptiae. o  -o  cap3d.ex  -llinpk 


"Cap2d.ex"  has  a  siailar  "make"  file  and  "source"  file  structure, 
with  names  changed  to  2d,  accordingly.  The  object  code  "coaptiae. o" 
is  the  "object"  file  of  the  two  C  procedures  which  are  utilized  for 
the  tiae  usage  calculation,  "-llinpk"  indicates  that  "linpack"  is  to 
be  loaded,  and  this  flag  has  to  be  the  last  flag  in  the  command. 
Finally,  the  output  will  be  in  "fort. 7"  for  "cap2d.ex"  and  "fort. 9" 
for  "cap3d.ex". 


The  calculation  of  capacitances  in  three-layer  cases  is  very 
computationally  intensive,  due  to  the  infinite  summation  and  the  dou¬ 
ble  suamation  inside.  A  few  hundreds  of  seconds  of  CPU  tiae  on  the 
UNIX  systea  is  not  uncommon  for  a  single  trace  (see  example  IV).  So 
run  this  case  cautiously  and  prudently. 


In  the  following  are  soae  examples  of  interconnect  systems.  The 
configurations  of  interconnects  are  presented  in  graphs,  and  the 
interactive  sessions  of  running  the  programs  and  the  results  from  the 
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run  are  also  included.  Hopefully,  these  examples  will  illustrate  the 
usage  of  both  programs  in  more  detail. 

£•£•!.<  Example  J. 


Following  is  the  interactive  session  of  running  the 
program  wcap3d"  for  the  traces  in  Fig.  B.l. 

% 

%  cap3d.ex 

number  of  traces  ntrc  ■  2 

number  of  layers  nlay  =  1 

want  kernel  matrix  Cy'  or  'n')  ?  ' n* 
want  grids  ('y*  or  'n* )  ?  'n* 
trace  no.  1  : 
orientation  (x*l,y»2)  »  1 
subdivision  (constmO,  Cheby“other)  -  1 
x  width  and  no.  of  div.  :  wx,nx  *  6.5dO  4 

y  width  and  no.  of  div.  wy,ay  ■  1.3d0  3 

z  width  and  no.  of  div.  wz,nz  *  1.2d0  3 

position  of  center  :  cx.cy.cz  ■  O.dO  O.dO  l.ldO 

trace  no.  2  : 
orientation  (x»l,y»2)  »  2 
subdivision  (const*0«  Cheby*other)  =  1 
x  width  and  no.  of  div.  wx.nx  *  1.8d0  3 

y  width  and  no.  of  div.  :  wy.ny  »  9.d0  4 

z  width  and  no.  of  div.  wz.nz  ”  1.7d0  3 

position  of  center  :  cx.cy.cz  >  O.dO  O.dO  3.35d0 

dielectric  const.  :  enl  ”  4.d0 
STOP  in  tort. 9  statement  executed 

END  of  the  interactive  session. 


The  output  from  the  above  session  is  as  follows  : 


*•••*  three-dimensional,  1-layer,  2-trace  •*•*•  Apr  1  83  00:17:29 

subdivisions  :  (x,  y  and  z) 

.44  2.81  2.81  .44  .19  .92  .19 

.18  .85  .18 


trace  no.  1  with  wx  =  6.50  nx  ■  4  wy 
wz  ■  1.20  nz  ■  3  center  at  :  .000 

orientation  :  1  subdiv.  method  :  1 


1.30  ay  *  3 

.000  1.100 


subdivisions  :  (x,  y  and  z) 

.26  1.27  .26  .60  3.90 


3.90 


.60 


4 


trace  no.  2  with  wx  ®  1.80  ax  *  3  wy  ■  9.00  ny  * 
wx  *  1.70  az  -  3  ceater  at  :  .000  .000  3.350 

orieatatioa  :  2  sabdiv.  method  :  1 

1-layer  with  diel.  coast,  eal  «  4.00 


total  ao.  of  sabsectioaa  : 


iterm 

■ 

0 

no.  of 

divgnt 

areas 

• 

e 

.08 

.54 

.54 

.08 

.40 

.54 

.08 

.08 

.54 

.54 

.08 

.54 

.54 

.08 

.08 

2.39 

.37 

.08 

.49 

.49 

.37  2 

.39 

2.39 

.37 

.08 

.16  1 

.03 

4.96 

1.03 

1.03 

.16 

.77 

.16 

1.03 

4.96 

.77 

.16 

.15 

.97 

.97 

.15 

.97 

.97 

.15 

.15 

4.68 

.72 

.15 

.97 

.97 

ao. 

0 

of  kernels  : 

4656 

2.59 

2.59 

.40 

.08 

.54 

.08 

.40 

2.59 

2.59 

.40 

.49 

.49 

.08 

.37 

2.39 

.08 

.08 

.49 

.49 

.08 

.49 

.49 

.08 

.16 

.77 

4.96 

1.03 

.16 

.77 

.16 

1.03 

1.03 

4.96 

1.03 

.16 

.15 

.72 

4.68 

4.68 

.72 

.97 

.15 

.97 

.15 

.72 

4.68 

rcond  -  .21049453 e-01  iafo  «  0 


solntion(even)  : 
1.447  6.351 

6.351 

20.83 

5.809 

1.447 

.7007 

.6746 

.6746 

.6321 

2.876 

.7007 

1.210 

4.889 

4.889 

6.229 

3.374 

.6668 

1.210 

4.889 

4.889 

6.229 

3.374 

.6668 

1.277 

5.118 

1.277 

2.299 

5.510 

2.299 

1.016 

4.028 

1.016 

2.639 

6.598 

2.639 

1.196 

2.281 

2.281 

6.039 

4.189 

.9724 

1.196 

2.281 

2.281 

6.039 

4.189 

.9724 

1.447 

5.809 

20.83 

6.351 

6.351 

1.447 

.7007 

2.876 

.6321 

.6746 

.6746 

.7007 

1.210 

3.374 

6.229 

.8890 

.8890 

.6668 

1.210 

3.374 

6.229 

.8890 

.8890 

.6668 

2.299 

5.510 

2.299 

1.277 

5.118 

1.277 

2.639 

6.598 

2.639 

1.016 

4.028 

1.016 

1.196 

4.189 

6.039 

2.594 

2.594 

.9724 

1.196 

4.189 

6.039 

2.594 

2.594 

.9724 

ao.  of  traces  :  2 

cap(  1.  1)  ■  166.0239059  cap(  2,  2)  ■  140.5242337 

cap(  1.  2)  -  51.89136663 


TIME  OSAGE  :  kernel  setup 
kernel  solve 
capacitance 
total  aser 


time  :  346.467 
time  :  5.93333 
time  :  .283333 
time  :  352.683 
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end  of  report 


Apr  1  83  00:27:43 


B.£.2.  Er  enrol  a  2_ 


The  interactive  session  and  outpnt  of  running  "cap2d”  for  the 

traces  in  Fig.  B.2  : 

cap2d.ex 
number  of  traces  ntrc  *=  4 

number  of  layers  nlay  “  2 

want  kernel  matrix  ('y*  or  'n')  ?  *n' 
want  grids  ('y'  or  *n')  ?  #n’ 
trace  no.  1  : 

horiz.  width  and  no.  of  div.  :  wx.nx  =  1.5d0  3 

verti.  width  and  no.  of  div.  :  wy.ny  “  .44d0  3 
position  of  low  left  corner:  lolx.loly  311  2.d0  ,38d0 
slant  angle  of  the  trace  :  ang  (in  deg)  -  7.0dl 
division  method  (consta0  Cheby=other)  *  1 
trace  no.  2  : 

horiz.  width  and  no.  of  div.  :  wx.nx  =  1.5d0  3 

verti.  width  and  no.  of  div.  wy.ny  *  .44d0  3 

position  of  low  left  cornet:  lolx.loly  ■  4.5d0  .38d0 
slant  angle  of  the  trace  :  ang  (in  deg)  »  7.dl 
division  method  (const~0  Cheby«other)  »  1 
trace  no.  3  : 

horiz.  width  and  no.  of  div.  wx.nx  3  2.2d0  3 

verti.  width  and  no.  of  div.  wy.ny  =  1.7d0  3 

position  of  low  left  corner:  lolx.loly  3  O.dO  2.37d0 
slant  angle  of  the  trace  :  ang  (in  deg)  *  7.5dl 
division  method  (const~0  Cheby™other)  ■  1 
trace  no.  4  : 

horiz.  width  and  no.  of  div.  :  wx.nx  =  2.2d0  3 

verti.  width  and  no.  of  div.  :  wy.ny  ■  1.7d0  3 

position  of  low  left  corner:  lolx.loly  ■  3 .2d0  2.37d0 
slant  angle  of  the  trace  :  ang  (in  deg)  *  7.5dl 
division  method  (const*0  Cheby*other)  ■  1 
diel.  const,  of  lower  layer  :  enl  *  4.d0 
diel.  const,  of  higher  layer:  en2  ■  7.5d0 
interface  position  :  h  *  1.07d0 
relative  error  bound  :  errbnd  (l.d-3)  *  l.d-3 

STOP  in  tort. 7  statement  executed 

END  of  the  interactive  session. 


two-dimensional,  2-layer,  4-trace  ***•*••  Aprl5  83  01:56:11 


trace  no.  1  with  wx  *  1.50 
position  :  2.00  .38 

subsections  : 


nx  *  3  wy  =  .44  uy=3 

slant  angle  (deg)  :  70.00 

.834  .173  .220  1.061 


173 


.220 
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.069  .331 

.069 

trace  no.  2  with 

wx  =  1.50  nx 

position  :  4.50 

• 

U) 

00 

subsections  : 

.173 

.069  .331 

.069 

trace  no.  3  with 

wx  ”  2.20  nx 

position  :  .00 

2.37 

subsections  : 

.189 

.258  1.244 

.258 

trace  no.  4  with 

wx  -  2.20  nx 

position  :  3.20 

2.37 

subsections  : 

.189 

.258  1.244 

.258 

=  3  wy  =  .44  ny  =  3 

slant  angle  (deg)  :  70.00 

.834  .173  ,220  1.061  .220 


=3  wy  »  1.70  ny  *  3 

slant  angle  (deg)  :  75.00 

.911  .189  .322  1.556  .322 


“3  wy  =  1.70  ny  =»  3 

slant  angle  (deg)  :  75.00 

.911  .189  .322  1.556  .322 


2-layer  with  brdy  h  =  1.07  rel  errbnd  »  .100e-02 

diel.  const,  enl(low)  »  4.00  en2(top)  *  7.50 

total  no.  of  subsections  :  48  no.  of  kernels  :  1176 

iterm  *  7  no.  of  diwgnt  :  0 


widths  : 

22  1.06 

.22 

33  .07 

.22 

07  .07 

.33 

26  1.24 

.26 

91  .19 

.26 

.17 

.83 

.17 

1.06 

.22 

.17 

.07 

.32 

1.56 

.26 

1.24 

.26 

1.24 

.26 

.26 

.07 

.33 

.07 

.83 

.17 

.07 

.32 

.19 

.91 

.32 

1.56 

.32 

1.24 

.26 

.07 

.33 

.19 

.19 


rcond  ■  .24604396e-0l  info  ■  0 

solution( even)  : 


2.9S2 

11.14 

2.911 

.2795 

1.387 

1.317 

.2139 

1.248 

2.910 

11.14 

3  .008 

.1394 

1.249 

.9512 

.1137 

1.441 

1.759 

2.804 

.1911 

.6801 

1.494 

2.748 

.7878 

.9962e- 

-01 

.1260 

.7004 

.6998 

.3601 

•7113e-01 

.3019 

.3588 

.8920 

no.  of  traces  :  4 

eap(  1,  1) 

-  23.02923288 

cap(  2,  2) 

* 

cap(  3,  3) 

-  12.85232106 

cap(  4,  4) 

ss 

cap(  1,  2) 

-  1.766451058 

cap(  l,  3) 

= 

cap(  1,  4) 

-  4.5225 

07157 

cap(  2,  3) 

SB 

cap(  2.  4) 

*  7.549320560 

cap (  3,  4) 

3 

USAGE  :  kernel  setup 

tine  :  198.850 

kernel  solve 

t  ime  :  1 

.15000 

capacitance 

t  ime  :  . 

233333 

.3718 

.9462 

.7297 

1.565 

1.291 

.2684 

1.254 

2.231 


24.11008633 
f  386332101 
4./19475070 
.2622865412 
15.43729017 


.1257 

.1012 

.5101 

.3493 

.3642 

.3656 

.6398 

.7506 


#  * 


total  user 


time  :  200.233 


end  of  report 


Aprl5  83  02:03:17 


B.5,.3,.  Example  3. 

The  interactive  session  and  ontpnt  of  running  Mcap3dM  for  the 
traces  in  Fig.  B.3  : 


cap3d.ex 
number  of  traces  ntrc  ”  2 

number  of  layers  nlay  =*  2 

want  kernel  matrix  (*y'  or  'n')  ?  'n* 
want  grids  ('y'  or  'n')  ?  *n* 
trace  no.  1  : 
orientation  (x**l,y»»2)  *  1 
subdivision  (const~0,  Cheby**other)  =  1 
x  width  and  no.  of  div.  wx.nx  =  S.dO  3 
y  width  and  no.  of  div.  :  wy.ny  *  2.d0  3 

z  width  and  no.  of  div.  :  wz,n z  ■  O.dO  0 

position  of  center  :  cx.cy.cz  =*  3  .ldO  3.d0  1.2d0 
trace  no.  2  : 
orientation  (x-l,y*2)  ■  1 
subdivision  (const*0,  Cheby*other)  «  1 
x  width  and  no.  of  div.  wx.nx  ■  5.d0  3 
y  width  and  no.  of  div.  :  wy.ny  ■  3.d0  3 

z  width  and  no.  of  div.  :  wz.nz  *  O.dO  0 

position  of  center  :  cx.cy.cz  *  3.1d0  6.5d0  3 ,2d0 
diel.  const,  of  lower  layer  :  enl  =  4.d0 
diel.  const,  of  higher  layer:  en 2  *  7.5dO 
interface  position  :  h  *  2.d0 
relative  error  bound  :  errbnd  (l.d-3)  *  l.d-3 

STOP  in  fort. 9  statement  executed 

END  of  interactive  session. 


•••••••  three-dimensional,  2-layer,  2-trace  •••*•••  AprlS  83  01:38:41 

subdivisions  :  (x,  y  and  z) 

.73  3,54  .73  .29  1.41  .29 

trace  no.  1  with  wx  ■  5.00  nx  *  3  wy  ■  2.00  ny  =»  3 

wz  »  .00  nz  ■  0  center  at  :  3.100  3.000  1.200 

orientation  :  1  subdiv.  method  :  1 


subdivisions  :  (x,  y  and  z) 
3  3.5 


161 


l 


*N 


trace  no.  2  with  wx  ■  5.00  nx  *  3  wy  *  3.00  ny  •  3 

wz  -  .00  nz  ■  0  center  at  :  3.100  6.500  3.200 

orientation  :  1  subdiv.  method  :  1 

2-layer  with  brdy  h  ■  2.00  rel  errbnd  *  .100e-02 

diel.  conat.  enl(low)  -  4.00  en2(top)  «  7.50 


total  no.  of  subsections  :  18  no.  of  kernels  :  171 

item  *  6  no.  of  divgnt  :  0 


areas  : 
.21  1.04 

1.55  .32 


.21  1.04 

1.55  7.50 


5.00  1.04 

1.55  .32 


.21  1.04 

1.55  .32 


.21  .32 


rcond  *  .53305083 e-01  info  *  0 

solution(even)  : 


4.092 

14.04 

4.092 

9.938 

27.23 

9.938 

3.553 

11.40 

3.553 

6.053 

16.18 

6.053 

17.37 

35.94 

17.37 

7.045 

20.92 

7.045 

£ 


no.  of  traces  :  2 

eap(  1,  1)  -  87.83536715  cap(  2,  2)  -  133.9861517 

cap<  1,  2)  -  15.45701657 


TIME  USAGE  :  kernel  setup  time 
kernel  solve  time 
capacitance  time 
total  user  time 

••••••••••  end  of  report  •* 


:  65.4833 
:  .183333 
:  .166667e— 01 
:  65.6833 

•••••**  Aprl5  83  01:42:41 


1.4.1.  fatfpip  1 

The  interactive  session  and  the  output  for  the  trace  in 
Fig.  B.4  : 

% 

%  cap2d.ex 

number  of  traces  ntrc  ■  1 

number  of  layers  nlay  =  3 

want  kernel  matrix  ('y*  or  'n* )  7  ’y' 
want  grids  ('y'  or  'n* )  7  'n' 
trace  no.  1  : 

horiz.  width  and  no.  of  div.  wx.nx  ”  1.7d0  3 
verti.  width  and  no.  of  div.  wy.ny  *  1.2d0  3 
position  of  low  left  corner:  lolx.loly  -  0,d0  .5d0 
slant  angle  of  the  trace  :  ang  (in  deg)  **  7.5dl 


^  ^  m 
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division  method  (const«0  Cheby*other)  *  1 
diel.  const,  of  lower  lsyer  :  enl  *  4.d0 
dial,  const,  of  higher  lsyer:  en2  -  7.5dO 
interface  positions  :  a,b  (s<b)  »  0.5d0  2.5d0 
relative  error  bonnd  :  errbnd  (l.d-3)  *  l.d-2 
STOP  in  tort. 7  statement  executed 

END  of  interactive  session. 


***••*•  two-dimensional ,  3-layer,  1-trace  ••**••*  Aprl5  83  00:02:36 

trace  no.  1  with  wx  *  1.70  nx  *  3  wy  *  1.20  ny  *  3 

position  :  .00  .30  slant  angle  (deg)  :  75.00 

subsections  :  .155  .747  .155  .249  1.202  .249 

.182  .878  .182 

3-layer  with  bdry  heights  a  *  .50  b  •  2.50 

diel  const  enl(low)  «  4.00  en2(top)  *  7.50 

relative  errbnd  »  .100e-01 

total  no.  of  subsections  :  12  no.  of  kernels  :  78 

e 


iterm  * 

12  no. 

of  divgnt  : 

0 

widths  : 

.25  1.20 

.25 

.15  .75 

.15  .18 

.88 

18  .18 

.88  .18 

packed  symm  matrix  ker  : 

.  84743 e-01 

.20391e-01 

•40871e-01 

.3 7890 e-02 

.20391e-01 

.84743 e-01 

.18020e— 01 

.16603 e-01 

.11365e— 01 

,99889e-01 

.15007 e-01 

.17498e-01 

.15007 e-01 

.47521e— 01 

,65681e-01 

.11365e-01 

•16603e-01 

,18020e-01 

.27945e-01 

.47521 e-01 

•99889e— 01 

.  62266e-01 

,17036e-01 

.37810e-02 

.193 81 e-01 

.15633e-01 

•11525e-01 

.883 93 e-01 

.30659e-01 

.186 03 e-01 

•8l513e-02 

.35820e— 01 

.26691 e-01 

•18518e-01 

.35191e— 01 

.55031e-01 

.19202 e-01 

.16595e-01 

.10561e-01 

.70566e-01 

.3 97 81 e-01 

,24433e-01 

.20815e-01 

.40914e— 01 

.95082e-01 

.37810e-02 

.17036e— 01 

.62266 e-01 

.11525e-01 

•15633e— 01 

.193 81 e-01 

.33764e-02 

.773 81 e-02 

.10605 e-01 

.883 93 e-01 

.81513e-02 

.18603 e-01 

.30659e-01 

.18518e— 01 

.26691e-01 

.35820e-01 

. 77381 e-02 

.123 57 e-01 

.16755e-01 

.3  5191e— 01 

.55031e— 01 

,1056le-01 

,16595e-01 

.19202 e-01 

.25121e-0l 

.40466e-01 

.68090 e-01 

,10605e-01 

.16755e-01 

.223l2e-01 

,20815e-0l 

.40914e-01 

,95082e-01 

rcond  - 

.61052279e-01  info 

-  0 

solution(even)  : 

2.958 

11.40 

2.958 

.9070 

2.174 

.9881 

3.017 

5.438 

1.200 

3.017 

5.430 

1.156 

no.  of  traces  :  1 

cap(  1,  1)  -  40.64073429 


COMPUTATION  OF  TWO-DIMENSIONAL  KERNEL  MATRIX  ELEMENTS 


C.i.  Introduction 

In  reference  to  Eq.  (2.14)  the  kernel  matrix  element  A^  is  a 
doable  integral  containing  the  Green's  function,  the  basis  function 
and  the  testing  function.  In  order  to  have  a  closed  form  for  the 
evaluation  of  these  elements,  the  basis  functions  are  chosen  as 
pulse-type  functions,  and  the  testing  functions  are  the  same  as  the 
basis  functions  for  Galerkin's  method.  Two  types  of  interconnect 
side  vails  are  considered  in  the  tvo-dimensional  case,  i.e.,  vertical 
and  slant  side  vails.  Of  course,  the  vertical  case  can  be  considered 
as  a  special  case  of  the  slant  side  vail  vith  a  pitch  of  90°.  For 
the  sake  of  comparison,  both  cases  are  considered.  Additionally,  let 
us  define  a  special  convention  to  simplify  the  expressions  of  those 
closed  forms. 

[f(x,y)l  a*'  4  ”  f(a,,y)  -  f(ax,y)  +  f(a4,y)  -  f(a,,y)  (C.l) 

[f(x#y)l  ■*  —  f(a»  ,Pi)  —  f(ax,px) 

1  y-Pi 

+  (C.2) 

vhere  f(x,y)  can  be  any  real-valued  function,  and  ax,  a* ,  a,.  a« ,  0X, 
Px  are  real  constants. 
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C.2.  Vertical  Side  Wall 

For  the  two-dimensional  and  vertical  side  wall  case,  the  basis 
functions  which  divide  the  perimeter  of  the  cross-section  of  an 
interconnect  into  subdivisions  can  generally  be  categorized  into  two 
groups  :  horizontal  pulses  and  vertical  pulses.  They  are  defined  as 


p$(x.y>  -  [  1/SJ'  for  xi  1  1  *  xl-  7 
0,  otherwise 

and  Sj  m  <Xj“Xj) 


pj<,.y>  -  [ 1/Sj'  ,or  i  7 1  1 '  *5 

0,  otherwise 
and  S.  -  (y]-yj)* 


(C.3) 


(C.4) 


Since  the  "cores"  of  the  two-dimensional  Green's  functions  in 
Eqs.  (2.15,2.16,2.21,2.22)  are  logarithmic,  the  following  integrals 
result  from  incorporating  the  pulses  in  Eq.  (C.3)  and  Eq.  (C.4). 
Assume  that  2  *j.  *i  2  *7-  7j  2  7j*  7j  2  7~i  A,  K  are  real  con- 

o 

stants  which  can  also  be  zeros;  moreover,  z^  is  the  abscissa  of  the 
vertical  pulse  P*,  and  yj  is  the  ordinate  of  the  horizontal  pulse  P^. 


(1)  I 


zz 


J  ^dz  j*  dz'  ln[(z-z')#  +  A*] 

X7  Xi 


*  Normalizing  the  pulse  aaplitude  with  the  pulse  width  or  area. 
Si,  is  essentially  the  same  as  using  the  charge  Q  as  the  "state  vari¬ 
able"  as  described  in  [10].  The  advantage  is  a  numerically  more 
stable  kernel  matriz. 
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(2)  I. 


(3)  I 


(4)  I 


J7Jdy  J7*dy'  ln[A*  +  (y+y'+K)*] 

7j  yi 

J^dy  J^dy'  la[A*  +  (y-y*  +K)*1 

y7  n 

j*y^dy  J7^dy'  ln[A*  +  (ly-y'l+K)*] 

7I  n 

f  1  «•  f  t  *  —  *  ^  \  ^  X  /  ^  Xwl  -X 1 


(5)  Ixyl  -  JXjdx  J7  dy'  ln[(x-xj)‘  +  (yj+y'+K)*] 


Xj  yi 


(6)  -  f^dx  J7^dy'  ln[ (x-x® )  *  +  (yJ-y'+K)*] 


Xj  yi 

+  + 


(7)  I 


J  ^dx  Jyidy*  ln[(x-xj)*  +  ( lyj-y' l+K) *] 

xI  n 


Because  of  tlie  similarity  of  the  above  integrals,  the  closed- 
form  expressions  for  them  can  be  given  in  terms  of  the  following 
functions  : 


F(x.y)  «  -y-(x*-y*)  ln(x*+y*) - 1—  x*  +  2xytan-1-y- 


(C.5) 


H(x,y)  ■  xyln(x*+y*)  -  4xy  +  2y*tan-1-p- 


(C.6). 


G(a,p,a)  »  xtan-1~-rix. 


?2« 


0.3l(a*+x*)tan_1-|-  +ax]  ^ ,  0>a>O  or  02P>a 

xma  (C.7) 

A 

0.5t(a*+x,)tan“1-i-  +»x]  - 2-a*sign(a),  a<O<0 


x*a 


where  sign(a) 


I  a  I 


Note  that  F(x,y)  is  an  even  function  in  both  variables  x,y,  while 
H(x,y)  is  an  odd  function  on  both  variables.  The  special  considera¬ 
tion  given  to  G(a,0,a)  as  described  above  is  due  to  the  singularity 
or  branch  cut  at  x-0. 

After  sone  algebraic  manipulation,  the  integrals,  1^,  i  r 
etc.,  can  be  expressed  in  terms  of  Eqs.  (C.5,C.6,C.7)  as  follows  : 


(1) 


I 

xx 


(2)  I 


.  rF(.  A)i  t1**P*- 
711  1  VyjTl-i:. 


(3)  I 


772  L  J7-7j-y]+K,  y]-yt+K 


r  l  iyt-7iI+K'  iyT-yll+K 

L  Vlyj-y^l+K,  ly~-yt|+K 

yj>y]>yi>yi  or  yt>yl>yj>yj 

(it)  [*<’■»]*’-'’  7tT  -  KtJtJ). 

1  Vy Ti-7j+i  J  1 


1  2 . 


(Eln(I  +A")-2E+2Atan"’1-^— )  ,  y?>yt>y->y- 


(4)  Iy73  -  ( iii)  interchange  yj<_>yi.  yj<“>yt 

in  thr  above  expxession,  if  yj>yi>y"i^yj 


r  i  yT“yI+K'  yi“ yj+K 

[F<y,A)]  7i  *  ll  -2(yW7). 
L  VyryJ+K,  7i-yj+K 


(Kln(r+A')-2K+2Atan~1-|-)  ,  yt>yt>y7>y- 


(▼)  interchange  yj<->yji»  yj<~>yi 

in  the  above  expression,  if  yj>yi>yj>yi 


(5)  I 


[[H(x.y)]  *  +  2G(xt-x^,Xj-Xj,y) ]  7j  ?1 

x-xj-xi  Vyj+yi+K 


(6)  I. 


xt-*! 


yj“y!+K 


-  [«*•»>  Wj  * 

J  ^  J  1 


xt-x^ 


sign(y-T-y*)[[H(x,y)]  j  J 
J  x-x7-Xj 


(7)  I. 


yj>yl  °x  yj<yi 


[[H(x,y)J  ^  o  +  2G(Xj-x^,xJ-x^,y) ]  7*  7j 

x-xJ-xA  Jy*5  K»  K 

yi>y]>y'i 


£•2*  Slant  Side  Wall 


If  the  pitch  angles  of  the  side  walls  of  the  interconnects  are 
taken  into  account,  assnaing  that  the  interconnects  are  left  and 
right  syaaetric,  the  pnlse-type  basis  functions,  be  they  on  the  hor¬ 
izontal  planes  or  on  the  slant  side  walls,  generally  can  be  expressed 
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parametrically  with  the  help  of  polar  coordinates  as 


oV/  ,  r  1/A  { *  for  x-x.+tcosrii.  y»y .+tsinq  . ,  0  £  t  £  A. 
P,(x.y)  “  I  J  JJJJ  J 


0,  otherwise 


(C.7) 


Note  that  (x^.y^)  is  the  coordinate  of  one  of  the  end  points  of  the 

subdivision,  and  is  the  polar  angle  of  the  subdivision.  This 

angle  should  be  limited  within  0°  and  360*.  A ^  is  the  width  of  this 

subdivision,  and  as  in  the  previous  case,  it  is  used  to  normalize  the 

o 

pulse  amplitude.  For  horizontal  subdivisions,  the  angle  *  o  or 
e 

180  .  Hence,  all  the  subdivisions  on  the  two-dimensional  cross- 
section  of  an  interconnect  can  be  represented  by  Eq.  (C.7) . 

As  in  the  previous  section,  after  incorporating  the  "cores"  of 
the  two-dimensional  Green's  functions  with  the  pulse  in  Eq.  (C.7), 
the  following  integrals  are  needed  for  the  evaluation  of  the  kernel 
matrix  elements. 


(1)  Ix  ■  P^dt  dt'  ln[  (x,+tcosq  ,-x^-t' cosq^) 

0  0  j  j 

+  (yj+tsin^j-y^t'  sinr^+k,)  ] 


Ix  *  j  ^dt  J  dt*  ln[(x .+tcosq .-x --t' cosq .) 

J0  J0  J  J  % 

+  (yj+tsinUj+y^+t'  sinq^+k#)  ] 


Ij  *  P^dt  J  *dt'  ln(  (x  .+tcosq  .-x  .-t '  cosq  • ) 

J0  J0  J  J 

+  (  l3rj+tsinqj-yi-t'sinqil+ko)1] 
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where  k0  i$  a  real  constant ,  and  Aj ,  A ^  are  the  lengths  of  the  jth 
and  the  ith  subdivisions  in  question. 

In  order  to  encompass  the  evaluation  of  all  three  of  the  above 
integrals,  let  us  first  consider  a  more  general  integral,  assuming  a, 
a,  b  and  k,  are  all  real  constants. 


I 


Aj  fa+bt  4 

dt  J  dt'  ln[(x  .+tcosTi --x,-t' cosn.) 

°  a  J  J 

+  (yj+tsinUj-yj- t' sinnj+k0)  ] 


■Aj  -a+bt-ki 
Jdt  I  dt'  ln( t 

0  a-kj 


*+<) 


(C.8a) 


where  kx  «  at  +  yx 

9 

k*  *  <*ct  +  r* 


(C.8b) 


tnd  as  *  “  ®c  **  cos^j-i^) 

y x  *  (xj-xi)siaT\i  -  (yj-yi+koJcosnj 
y »  M  (xj-xi)cosr\i  +  (yj-yi+k0)sinni 


(C.8c) 


Furthermore,  two  more  functions  are  defined  to  represent  the  integral 
in  Eq.  (C.8a).  They  are  : 


_  rAj  _■»  (b-o  )t+(a-r») 

PsfAj.a.b)  -  J  dt  (Qst+yi)  tan  - - 


(C. 9a) 


F*(A.,a,b)  *  J  dt  [  (b_a  t+(a— Yj)  ]  ln{[(b-a  )t+(a-yj)]  +(a  t+Y*)  } 
J  0 

(C.9b) 


Finally,  the  integral  I  in  Eq.  (C.8a)  can  be  vritcen  in  terms  of  the 


above  F*  and  F,  as 


I  “  2F1(Aj,a.b)  +  F,(Aj  ,a,b)  -  2F1(Aj,a,0)  -  F^Aj.a.O) 


-  2(a-a)Aj  -  bAj 


(C.10) 


After  some  lengthy  derivation,  the  integrals  in  Eq.  (C.9)  can  be 
expressed  in  closed  forms  as  : 


(i)  (-2-a"Aj  +  rxAj)taa' 


-1  (b“®c)Aj+(a-y,) 
Vj+Tx 


-,  if  r*  =  o 


_  a  t  _i  (b— a  )A  -  +  ( a-y ,) 

F^A-.a.b)  -  (ii)  (  -----  A  ■  +  TjAjtan  1 - — — J - 

J  2  J  J  asAj+Ti 


2T7  Aj 


T t  YtA;  +2y,A.+y, 

■- — J-(b-a  )  in  - J - J - 

2rT  r. 


(C.lla) 


1  ,  asT»  T«Y.(b-a  ) 

+  - [~* -  +  - # 

2  Tt 


,  -i  r.+T7Ai  r* 

(tan  A J - tan  x - )  ,  if  r.  a  o 

T«  ft  * 


(  i)  (»-T,)A.  Inly, | ,  if  r,  =  0 


(C.llb) 
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(b-ac) 


Fj(Aj,a,b)  *  (ii)  [ — 2  1 (b“°c) +U-*Y»)Aj]  lalrfl  +  2  ' 


r,  *  r«  a  r*  j-^j-hra/r? 

I(x  +  — r~)  ln(x  +  — 3-)-x  P  +  c  < a — JT a ) 

1  *1  **  j*-y./y* 


r» 

r7 


+y»/yt 


r  a  Y  *  T  •  _i  ay  t  -iA  i 

-(b-a  )  ]  |  xln(x  +  — r— )-2x+2 - tan  1 - 

c  L  v  r*  J*-r,/y, 


if  rT  jfe  0 


where  r«  -  (b-ac)Yi  -  (a-Y,)as 
Y7  -  a*  +  <b-ac>* 

Y.  *»  (b-ac)  (a-r,)  +  asYx 

rf  *  Y*  +  (*~Y»)*  (C.llc) 


Finally,  the  integrals  described  at  the  beginning  of  this  sec¬ 
tion,  Ix,  IJf  I,,  can  be  expressed  in  terms  of  the  more  general 
integral  I  which  in  torn  is  an  algebraic  combination  of  closed  forms 
Fx  and  F2 .  Assume  that  I  is  a  tonction  of  its  variables  and  can  be 
written  as 

I  -  I(Aj,a,a,b,xj,yi,T,.,xj,yj,nj,k0) 

Then  Ix,  lx  can  be  written  in  terms  of  I  as  follows. 
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"  l(Aj  ,0,A,»0,x,»y,,Ti,»*j  »7j  j'Hj  *k0) 

^2)  Xj  ■  I(A  .  ,OiA.  iO,x.  ,-y,  ,i|  •  ,x  .  ,y  .  ,ti  •  >kg) 

J  i  *  *  *  J  J  J 

Assuming  sirnij  2  0  and  sinn,  2  0,  I,  is  expressed  as 
(i)  7^  >  y,  +  Aisinni 


*»  **  I(Aj»°»Ai.O,xi,yi,iii,Xj,yJ.,Hj,ko) 


(ii)  1  yj  +  AjSinnj 


I(Aj  •0,A^,0,Xj,y^,H£,Xj,yj,Tjj  ,-k0) 


(iii)  y^AjSinn,  >  y^A^sin^  and  7j  2  y£  **d  sinn^O 

’  ‘  <Ai'°'  Jiniij 

_  T,A  .  _l£2i  sinTh 

VV  .inn, 

(  iv)  yi+Aisinx\i  >  yj+Aj.inr\j  and  jj  2  y,  »ad  sinn,  *  0 
X,  ■  I(Aj  »0,A,,0,x,,y,,Hj[»*j»yj»nj»A*) 


(v)  yj+AjSinnj  2  yi+Aisinniand  y,  2  7j  and  sinnj  t  0 

t  y^T^  sian, 

1  ■  IIA‘-0'  .i»n. 

.  T/*  .  7i~7\  sinTU 

*'V  .1- 


«  A  fin 


(vi)  yj+AjSinHj  2.  y^+A^im^and  yi  2  yj  aad  sim^  -  0 

I,  -  I(A^,0,Aj,0,X£,y£4Ti£,Xj,yj<Hj,ko) 


(vii)  yj+A^aini^  >  yj+Ajiinrij  >  yA  >  yj  and  sinrjj  t  0  and  sinijj  £  0 


I,  «  I< 


y^-y^A^sinn^  y^y*  *iar|  i 


.  *  sinn.  '7J 'V  k0) 


sinr\^  sinrij 


JC7\ 


Y  i~7  \  *ia^i 


»inn.  +AJ  sini^  ’AJ*  sini^  '  sinn  ' x  i  * y  i ' 71  i ' x  j  'yj  »k*) 


+  KA..O.A.-A.- 

J  A  J 


4£^.0.«i.5ri.,i.Xj.yj.V-k.) 


(viii)  y j+Aj *inHj  >  yi+Ai$inni  >  y^  >  yi  and  sini^  t  0  and  sini^  +  0 

T  7C7\  *imii  yryi  *iniii 

*j  M  I( . ■1"1  +A  . — - ,0 , . ^  .  , . **  ,x4  »y4  ,x- ,y  ■  ,n  •  ,k0) 

sim\4  1  sm^.  smi^  sini^  11  1  J  J  J 


yty, 


ain^lj  y^-yj  sinri 


sini^  smi) 


” — *“»xi.yi.ni»x:  >Yi  ,n4 ,-k0) 
sinn.  1  1  1  J  J  J 


,  sinnj 

+  KA.-A. - J 

J  *  sinTj. 


— : - ■1—,0>Ai,0,xi,yi,ni»x  .  ,y .  ,11 .  ,k0] 

s  mn .  1  1  1  1  J  J  J 


With,  the  help  of  the  above  closed-fora  expressions,  the  evaluation  of 
the  kernel  matrix  elements  easily  can  be  obtained. 


EVALUATION  OF  THREE-DIMENSIONAL  KERNEL  MATE II  ELEMENTS 


D.l,.  Introduction 

Following  the  sue  strategy  outlined  in  Appendix  C  for  the  two- 
dimensional  case  and  employing  the  same  convention  in  Eq.  (C.1,2),  we 
can  evaluate  the  three-dimensional  kernel  matrix  elements  in  closed 
form  as  before.  However,  only  the  vertical  side  wall  case  is  con¬ 
sidered  here. 

D.2.  Formulation 

As  before,  in  reference  to  Eq.  (2.14),  the  pulse  type  basis 
functions  on  the  surface  of  the  interconnects  are  chosen  as 


Pj (x.y.z)  -  [ 


1/Sj ,  for  xj  <  x  i  xj , 
0,  otherwise 


yj  <  j  i  z 


(D.la) 


pju.7.,>  -  [  1/sr  for  <  r  i  rj.  *J  <  *  i  *  -  xj 


otherwise 


(D.lb) 


P?(x.y,z)  * 


£  1/Sj,  for  zj  <  z  £  ij,  xj  <  x  £  x j , 


otherwise 


(D.lc) 
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where  is  the  area  of  the  subdivision,  e.g.,  *  ( x j -x J )( y j- -y J )  in 

Eq.  (0.1a).  Similarly  S.  can  be  obtained  for  Eq.  (D.lb)  and 

J 

Eq.  (D.lc) .  The  reason  for  using  as  the  normalization  factor  of 
the  pulse  amplitude  is  to  have  a  more  stable  kernel  matrix  as  dis¬ 
cussed  in  Appendix  C.  Note  that  in  Eq.  (D.la)  xT,  xf,  f~i,  y+,  are 

J  J  J  J 

the  bounds  of  the  jth  rectangular  subdivision,  and  Zj  is  the  ordinate 
of  this  "square"  on  the  z-axis.  The  pulses  described  in  Eq.  (D.l) 
are  on  the  planes  parallel  to  the  xy,  yz,  zx  planes,  respectively,  as 
the  superscripts  indicate. 

Since  the  "cores"  of  the  three-dimensional  Green's  functions  in 
Eqs.  (2.23,2.24.2.29)  are  of  the  form 

f(  [(x-x')a+(y-y')*+(z±z'+k#)l]"1/*  ) 


some  of  the  typical  integrals  that  result  from  incorporating  them 
with  the  pulses  in  Eq.  (D.l)  are  listed  in  the  following  : 


(1) 


J^dx  |Xldx'  J^dy  J^dy*  [  (x-x' )  *+(y-y' )  *+(  Zj-z^+k,)  '] 

x7  xi  y7  y7 


(2)  I, 


(3)  I, 


[  (x-x' ) *+(yj-y*) *+( z+z'+k0) *]-1/2 
[(x-x' )  *+(y-y^)  *+(Zj-z'+kl>) 


ffimrmmjgmr 
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it  Z+  ^  Z  + 

(4)  I4  ■  J  dz  j"  dz'  j~^dx  J  dz'  [  (x-xp  *+(7j-y)  *+(  z-z*  +k#)  * ]-1/* 


Zj  zi  xj  *i 


where  ka  is  a  real  constant  which  can  also  be  zero. 


Next  two  closed-fora  functions  are  defined  such  that  those 
integrals  presented  above  can  be  expressed  in  teras  of  thea.  They 


P(o.p,$) 


1  a  a  a  a / a  2  a  a  ai/a 
-  — (a  +8  +«  )  '  +  V-(a  +P  +5  )  ' 

6  2 


°P  ,  _i  B  _i  « 

+  (asinh  . .  .  +  Bsinh  .  ■  ) 

j  a  a/a  H  a  .a  a/ a  7 

2  (<*+$)  (P  +5  )  ' 


5  B  -1  a 

——•(Bsinh  1  "  -  . . -  +  asinh  —  •  -  ■  ■  ■ ) 

a  *  x/s  a  a  a/ a  7 

2  (a  +5  )  7  (B  +5  ) 


-  gaptan 


5<«W>1/1 


(D.2) 


a  a  a  i/a  _i  5 

Q(o.p,5)  -  - (a  +B  +5  )  '  ~  5oP*iai  1 - : — TTTT 

3  (a  +0  )  ' 


,  B  « 

0.5[p(? - )  * inh  1 - — 

3  (B  +«  )  ' 


x  a  ,  P 

+a(? - )5iah  - ■  —  ;  t:,-I 

3  (a  +?  )  ' 
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+  0.551  P*taa-1 - - — --- —  +  o‘t 

p(a  +0  +5  )  ' 


5a 


53 


111  1/1 
a(a  +0  +5  )  ' 


1  t 

+  - 5  tan  ■ 

6 


afi 


(D.3) 


It  should  be  pointed  out  that  P(a,0,5)  is  an  even  function  on  all 
three  variables,  while  Q(a,0,5)  is  an  even  faction  on  the  last  vari¬ 
able  5  sad  an  odd  function  on  the  first  two  variables  a  and  0. 

Finally,  adopting  the  convention  in  Eqs.  (C.1,C.2),  we  can 
express  the  aforementioned  integrals  as  follows  : 

yt-yl  yT~yi 


(1) 

-  [p(a,0 , Zj— z^+k0) 

i  Tl 

xj"xi 

■■  x 

a-Xj-Xi 

*  xj”zi 

(2) 

li  m  [-P(a,0,yj-y^) ] 

+-X  + 

Xj  xi 

— 

xrxi 
—  + 

a-Xj-Xi, 

Xj-Xi. 

r+k. 


(3)  I,  -  [-Q(a.0.5)]  _ 


yj-yl  xH+k»  xj“xi  xrxi 

t-.T.  x--xt 
J  1 


a-y]-yi»  P»*j-zi+k#,  5-Xj-Xi 


r  1  xj_xi  yj-yi  zj+zi+ko  Sj+xi+ko 

(4)  I4  -  [Q(a,0,t)]  i9 


u*x 


e  - 


I++Z+H 

J  1 


j-*!*  0-7]-7i,  5** j+zi+k0 ,  zj+zt+k. 


There  are  some  other  possible  integrals  when  incorporating  the 
three-dimensional  Green's  functions  with  the  pulses  in  Eq.  (D.l) ,  but 
they  are  all  similar  to  those  just  presented.  Furthermore,  they  are 
all  expressible  in  terms  of  the  two  closed-form  functions  P  and  Q. 


Therefore,  the  kernel  matrix  elements  can  all  be  evaluated  by  the 
closed-form  expressions  described  above. 
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APPENDIX  E 


SPICE  INPUT  DECK  FOR  FULL  ADDER  I 


Shown  in  the  following  is  the  SPICE  input  deck  for  "Full 
Circuit  I"  with  ell  the  interconnect  paraaters. 


full 

adder 

circuit  I  (with 

rOOO 

15 

Id 

3 .de+01 

rOO1 

Id 

15 

2.0e+01 

r002 

15 

Id 

5.0e+01 

r003 

17 

18 

5.0e+01 

r004 

17 

19 

8. de+01 

r005 

18 

17 

7 .Oe+Ol 

r006 

17 

18 

1.0e+01 

cOOO 

500 

0 

0.0158pf 

r007 

21 

500 

1 .59+01 

rr007 

500 

22 

1.5e+01 

cOOl 

501 

0 

0.0158pf 

r008 

21 

501 

1.5e+01 

rr008 

501 

23 

1.5e+01 

c002 

502 

0 

0.0158pf 

r009 

21 

502 

1.5e+01 

rr009 

502 

24 

1.5e+01 

c003 

503 

0 

0.0158pf 

rOlO 

21 

503 

1.5e+01 

rrOlO 

503 

25 

1.5e+01 

c004 

504 

0 

0.012dpf 

rOll 

2d 

504 

1.5e+01 

rrOll 

504 

27 

1.5e+01 

cOOS 

505 

0 

0.012dpf 

r012 

2d 

505 

1.5e+01 

rr012 

505 

28 

1.5e+01 

cOOd 

50d 

0 

0.0112pf 

r013 

2d 

506 

2.1e-01 

rr013 

50d 

29 

2.1e-01 

c007 

507 

0 

0.0014pf 

r014 

29 

507 

1.5e+01 

rr014 

507 

30 

1.5e+01 

c008 

508 

0 

0 .012dpf 

r015 

31 

508 

1.5e+01 

rr015 

508 

32 

l.Se+01 

c009 

509 

0 

0.0144pf 

rOld 

31 

509 

5.4e*01 

Adder 


rr016 

509 

33 

5.4e+01 

cOlO 

510 

0 

0.0l26pf 

r017 

31 

510 

1.5e+01 

rr017 

510 

34 

1.5e+01 

cOll 

511 

0 

0.0ll6pf 

r018 

35 

511 

3 .0e+02 

rr018 

511 

18 

3 .0e+02 

c012 

512 

0 

0.0084pf 

r019 

36 

512 

2.7e+01 

rr019 

512 

16 

2.7e+01 

c013 

513 

0 

0.0087pf 

r020 

36  513 

8.2e+01 

rr020 

513 

37 

8.2o+01 

c014 

504 

508 

0.0l77pf 

c015 

504 

509 

0.0211pf 

c016 

504 

510 

0.0177pf 

c017 

505 

508 

0.0178pf 

c018 

505 

509 

0.0211p£ 

c019 

505 

510 

0.0l77pf 

c020 

506 

508 

0.0177pf 

c021 

506 

509 

0.0211pf 

c022 

506 

510 

0.0177pf 

c023 

507 

510 

O.OOOlpf 

c024 

504 

500 

0.0210pf 

c025 

504 

501 

0.0209pf 

c026 

504 

502 

0.0209pf 

c027 

504 

503 

0.0209pf 

c028 

505 

500 

0.0209pf 

c029 

505 

501 

0.0209pf 

c030 

505 

502 

0.0209pf 

c03l 

505 

503 

0.0210pf 

c032 

506 

500 

0.0209pf  • 

c033 

506 

501 

0.0209pf 

c034 

506 

502 

0.0209pf 

c035 

506 

503 

0.0209pf 

c036 

507 

502 

O.OOOlpf 

c037 

508 

513 

0.0034pf 

c038 

509 

513 

0.0036pf 

c039 

510 

513 

0.003 ipf 

mldOOO 

5 

16 

16  0  lotd 

1- 

12.00am 

V® 

2 .17am 

mldOOl 

5 

10 

10  0  load 

1- 

8.00am 

w® 

2 .25am 

mld002 

5 

18 

18  0  load 

1® 

17 .00am 

w® 

2.24am 

mid 00 3 

5 

11 

11  0  load 

1- 

10.00am 

V® 

2 .20am 

mot 000 

6 

25 

0  0  drive 

1- 

2.00am 

V® 

6 .00am 

motOOl 

7 

30 

0  0  dr  ive 

1- 

2.00am 

w® 

8.00am 

mot 002 

8 

23 

0  0  dr  ive 

1® 

2.00am 

V® 

8.00am 

mot003 

9 

24 

0  0  drive 

1* 

2.00am 

V® 

8.00am 

mot004 

16 

30 

90  dr ive 

1- 

2.00am 

w= 

8 . 00am 

motOOS 

16 

34 

8  0  drive 

1® 

2.00am 

w= 

8 . 00am 

mot006 

15 

34 

7  0  drive 

1® 

2.00am 

V® 

8 .00am 

oot007 

10 

16 

0  0  drive 

1® 

2.00am 

V® 

8.00am 
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mot008 

11 

35 

0 

0 

drive 

1= 

2 . OOum 

w* 

8. OOum 

mot009 

12 

27 

0 

0 

drive 

1- 

2 .OOum 

w*= 

8. OOum 

motOlO 

18 

32 

20 

0 

drive 

1- 

2. OOum 

v= 

6. OOum 

mo  tO  11 

17 

37 

12 

0 

drive 

1- 

2. OOum 

w*» 

8. OOum 

mo  tO  12 

18 

37 

14 

0 

drive 

1= 

2 .OOum 

w*» 

8. OOum 

mot013 

20 

28 

6 

0 

drive 

1- 

2 .OOum 

W“ 

6 .OOum 

mot014 

13 

33 

0 

0 

drive 

1- 

2. OOum 

w= 

8. OOum 

mot013 

19 

37 

13 

0 

drive 

1- 

2 .OOum 

8. OOum 

mot016 

14 

22 

0 

0 

drive 

1« 

2 .OOum 

W“ 

8. OOum 

*load  capacitances 

*  sum 

ell  11  0  .lpf 

•  carry 

cl2  10  0  .lpf 

.model  load  omos  vto*»-2.  uo»600  level*l  lambda=0.02  tox=l . e-7 
.model  pass  nmos  vto**l.  uo=»600  level*l  lambda=0.02  tox=l.e-7 
.model  drive  nmos  vto*l.  uo**600  level»l  lambda**) . 02  tox=l.e-7 
* input  s 

v3  32  0  pnlse(0  3  120n  O.ln  O.ln  120n  240n) 

v2  28  0  pulse(0  3  60n  O.ln  O.ln  60n  120n) 

vl  23  0  pnlse(0  3  30n  O.ln  O.ln  30n  60n) 

vdd  3  0  5 

.options  lvltimol  itl4-20 
.tran  3n  240n 

•print  tran  v<25)  v(28)  v(32)  v(ll)  v(10) 

.end 

Note  that  the  capacitors  numbered  from  nc014M  to  nc039”  are  coupling 
capacitances.  They  are  "commented  out*'  in  the  case  of  only  consider¬ 
ing  the  self-capacitances.  For  the  case  of  no  interconnect  parame¬ 
ters,  another  extraction  is  done  with  large  parameters  "RLIMIT”  and 
"CLIMIT*  specified  in  "constant. h" . 


183 


APPENDIX  F 


SPICE  INPUT  DECK  FOR  FULL  ADDER  II 


The  inpat  SPICE  deck  for  "Fall  Adder  Circuit  II"  is  given  in  the 
following. 


fall  adder  circuit  II  (with  capacitance) 


rOOO 

10 

7 

3.5e+01 

r002 

12 

13 

6 .9e+01 

r003 

12 

14 

1.0e+01 

cOOO 

500 

0 

0.0114pf 

r004 

14 

500 

5 .le+OO 

rr004 

500 

15 

5.1e+00 

cOOl 

501 

0 

0 ,0109pf 

r005 

14 

501 

5.  le+OO 

rr005 

501 

16 

5. le+OO 

c002 

502 

0 

0.0ll4pf 

r006 

17 

502 

4.9e+01 

rr006 

502 

18 

4.9e+01 

c003 

503 

0 

0.0335pf 

r007 

17 

503 

5.4e+01 

rr007 

503 

19 

5.4e+01 

c004 

504 

0 

0 .0114pf 

r008 

-17 

504 

8 . Se-02 

rr008 

504 

20 

8 . 6e-02 

r009 

20 

21 

1.0e+01 

cOOS 

505 

0 

0.0316pf 

rOlO 

22 

505 

5.2e+00 

rrOlO 

505 

8 

5.2e+00 

rOll 

23 

24 

5 .Oe+Ol 

r012 

23 

25 

l.le+02 

r013 

23 

26 

1.0e+02 

r014 

27 

28 

3 .Oe+Ol 

r015 

29 

28 

5.5e+01 

r016 

30 

31 

2.5e+01 

c006 

506 

0 

0.0427pf 

r017 

32 

506 

2 .0e+02 

rr017 

506 

33 

2 ,0e+02 

c007 

507 

0 

0.0339pf 

r018 

34 

507 

3.5e+02 

rr018 

507 

6 

3.5e+02 

c008 

508 

0 

0.0043pf 

r019 

35 

508 

4.5e+01 
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rr019 

508  31 

4.5e+01 

c009 

509 

0 

0.0144pf 

r020 

36 

509 

1 .5e+01 

rr020 

509  37 

1.5e+01 

cOlO 

510 

0 

0.0193pf 

r021 

36 

510 

1.5e+01 

rr02l 

510  38 

1.5e+01 

cOll 

511 

0 

0.0370pf 

r022 

36 

511 

3 .2o+02 

rr022 

511 

.  39 

3.2e+02 

c012 

512 

0 

0 ,0065pf 

r023 

40 

512 

6 ,8e+01 

rr023 

512 

41 

6.8e+01 

c013 

513 

0 

0.0289pf 

r024 

40 

513 

3 .0e+02 

rr024 

513 

42 

3 .0e+02 

c014 

514 

0 

O.OlOlpf 

r025 

40 

514 

l.le+02 

rr025 

514 

25 

1.1 e+02 

c013 

515 

0 

0.0296pf 

r026 

43 

515 

3.1e+02 

rr026 

515 

44 

3.1e+02 

c016 

516 

0 

0.0279pf 

r027 

45 

516 

9.08+01 

rr027 

516 

46 

9.08+01 

c017 

517 

0 

0.0179pf 

r028 

42 

517 

3 .Oe+Ol 

rr028 

517 

47 

3 .Oe+Ol 

c018 

518 

0 

0.0173pf 

r029 

44 

518 

1.8e+02 

rr029 

518 

28 

1.8e+02 

c019 

519 

0 

0.0065pf 

r030 

37 

519 

6.8e+01 

rr030 

519 

7 

6  .8e+01 

c020 

505 

503 

O.OOllpf 

c021 

505 

506 

0.0222pf 

c022 

505 

507 

0.0222pf 

c023 

505 

512 

0.0229pf 

c024 

505 

513 

0.0229pf 

c025 

505 

514 

0.0265pf 

c02  6 

503 

506 

0.0285pf 

c027 

503 

507 

0.0285pf 

c028 

502 

511 

0.0030pf 

c029 

502 

513 

0 ,0059pf 

c030 

503 

513 

0 ,0314pf 

c031 

504 

513 

0 ,0030pf 

c032 

516 

507 

0.0120pf 

c033 

506 

507 

0.0096pf 

c034 

510 

518 

0.0175pf 

c035 

511 

517 

0.0052pf 

c036 

517 

515 

0 ,0052pf 
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c037  500  515  0.0030pf 
c038  501  515  0.0030pf 


mldOOO 

5 

31 

31 

0 

load 

1  = 

14.00um 

w= 

5 .71um 

mldOOl 

5 

6 

6 

0 

load 

1= 

14.00um 

w* 

5 ,71um 

mld002 

5 

8 

8 

0 

load 

1= 

1 4.00um 

w* 

5 ,71um 

mld003 

5 

7 

7 

0 

load 

1- 

14.00um 

w= 

5 . 71um 

mld004 

5 

25 

25 

0 

load 

1= 

14.00um 

w~ 

5 ,71um 

mld005 

5 

28 

28 

0 

load 

1- 

12.00um 

V* 

6 .OOum 

motOOO 

6 

32 

0 

0 

drive 

1® 

4.00um 

w= 

8. OOum 

motOOl 

10 

34 

0 

0 

drive 

1® 

4.00um 

w= 

8 .OOum 

mot002 

7 

35 

0 

0 

drive 

1® 

4.00um 

w= 

8. OOum 

mot003 

8 

19 

0 

0 

drive 

1= 

4.00um 

W** 

8. OOum 

mot004 

21 

39 

16 

0 

pass 

1  = 

4.00um 

vr= 

4. OOum 

mot 005 

15 

41 

21 

0 

pass 

1  = 

4  .OOum 

▼=* 

4 .OOum 

mot006 

13 

43 

22 

0 

pass 

1« 

4.00um 

4. OOum 

mot007 

26 

45 

0 

0 

drive 

1- 

4. OOum 

w* 

8. OOum 

mot008 

24 

33 

0 

0 

drive 

1® 

4. OOum 

w=* 

8. OOum 

mo tOO 9 

27 

42 

0 

0 

drive 

1= 

4. OOum 

ir= 

12. OOum 

motOlO 

5 

38 

9 

0 

pass 

1® 

4. OOum 

w= 

4. OOum 

mo tOll 

18 

44 

9 

0 

drive 

1® 

4. OOum 

w* 

4. OOum 

mo  tO  12 

9 

47 

0 

0 

drive 

1® 

4. OOum 

4. OOum 

mot013 

29 

37 

0 

0 

drive 

i  — 
a 

4. OOum 

w= 

8. OOum 

mot014 

30 

46 

0 

0 

drive 

1® 

4 .OOum 

w= 

8. OOum 

•load  capacitances 

•  carry 

ell  9  0  .lpf 

•  STUB 

cl2  14  0  .lpf 

•  model  cards 

.model  load  nmos  vto=-3.  uo®600  level®l  lambda=0.02  tox=l.e-7 
.model  pass  maos  vto®l.  uo®600  level=l  lambda=0.02  tox=l.e-7 
.model  drive  nmos  vto*l.  uo®600  level=l  lambda=0.02  tox=l.e-7 
•inputs 

v3  19  0  pnlse(0  5  120n  O.ln  O.ln  120n  240n) 

▼2  46  0  pnlse(0  5  60n  O.ln  O.ln  60n  120n) 

vl  32  0  pulse(0  5  30n  O.ln  O.ln  30n  60n) 

vdd  5  0  5 
. tran  3n  240n 

.print  tran  v(32)  v(46)  v(19)  v(14)  v(9) 

.end 


Note  that  the  capacitances  from  c020  to  c038  are  coupling 
citances.  Some  of  the  interconnect  parameters  consist  of  resist 
only,  e.g.,  rOll  to  r016,  because  the  capacitances  of  diff 
interconnect  lines  are  neglected  in  the  consideration.  Other 
siderations  are  the  same  as  in  the  previous  case. 


capa- 

ance  s 

u  s  i  v-  tl 
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